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ABSTRACT 

The Cosmic Infrared Background (CIB) reflects the sum total of galactic luminosi- 
ties integrated over the entire age of the universe. Prom its measurement the red-shifted 
starlight and dust-absorbed and re-radiated starlight of the CIB can be used to deter- 
mine (or constrain) the rates of star formation and metal production as a function of 
time and deduce information about objects at epochs currently inaccessible to telescopic 
studies. This review discusses the state of current CIB measurements and the (mostly 
space-based) instruments with which these measurements have been made, the obstacles 
(the various foreground emissions) and the physics behind the CIB and its structure. 
Theoretical discussion of the CIB levels can now be normalized to the standard cos- 
mological model narrowing down theoretical uncertainties. We review the information 
behind and theoretical modeling of both the mean (isotropic) levels of the CIB and their 
fluctuations. The CIB is divided into three broad bands: near-IR, mid-IR and far-IR. 
For each of the bands we review the main contributors to the CIB flux and the epochs 
at which the bulk of the flux originates. We also discuss the data on the various quanti- 
ties relevant for correct interpretation of the CIB levels: the star- formation history, the 
present-day luminosity function measurements, resolving the various galaxy contribu- 
tors to the CIB, etc. The integrated light of all galaxies in the deepest near-IR galaxy 
counts to date fails to match the observed mean level of the CIB, probably indicating 
a significant high-redshift contribution to the CIB. Additionally, Population III stars 
should have left a strong and measurable signature via their contribution to the cosmic 
infrared background (CIB) anisotropies for a wide range of their formation scenarios, 
and measuring the excess CIB anisotropies coming from high z would provide direct 
information on the epoch of the first stars. 

Subject headings: cosmology: diffuse radiation - galaxies: clusters: general - galaxies: 
high-redshift - surveys 
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1. Introduction 

Diffuse backgrounds contain important information about tlie liistory of the early Universe, 
when discrete objects either did not exist or are not accessible to telescopic studies. The cosmic 
infrared background (CIB) arises from accumulated emissions from early galaxy populations span- 
ning a large range of redshifts. The earliest epoch for the production of this background occurred 
when star formation first began, and contributions to the CIB continued through the present epoch. 
The CIB is thus an integrated summary of the collective star forming events, star-burst activity 
and other luminous events in cosmic history to the present time. As photons move to observer they 
lose energy to cosmic expansion and any stellar emission from high redshift populations will now 
be seen in the infrared (mostly near-IR with A ^ a few fim , unless it comes from very cold stars at 
very high z). Emission from galactic dust will be shifted to still longer IR wavelengths. CIB thus 
probes the physics in he Universe between the present epoch and the last scattering surface and 
is complementary to its more famous cousin, the cosmic microwave background radiation (CMB) 
which probes mainly the physics at the last scattering. 

Considerable effort has now gone into studying the luminosity sources during the most recent 
history of the universe (z < 2 — 3), but the period from recombination to the redshift of the Hubble 
Deep Field (HDF) remains largely an unexplored era because of the difficulty of detecting many 
distant galaxies over large areas of the sky. Significant progress in finding high redshift galaxies has 
been achieved using the Lyman dropout technique (Stcidcl ct al 1996), but uncovering substantial 
populations of galaxies at z^5 is extremely difficult and the total sample is still very small. It is 
there that the CIB measurements can provide critical information about the early history of the 
Universe largely inaccessible to telescopic studies. At more recent epochs galaxy evolution is also 
constrained by the measurements of the visible part of the extragalactic background light or EBL 
(Bernstein et al 2002a,b; cf. Mattila 2003). 

Chronologically, the importance of the CIB and early predictions about its levels (Partridge & 
Peebles 1968) followed the discovery of the CMB which observationally established the Big Bang 
model for the origin and evolution of the Universe. Bccrtisc of the Earth's atmosphere is was clear 
from the start that the measurement of the CIB must be done from space and even there the Solar 
system and Galactic foregrounds presented a formidable challenge. It took a while for technology 
to reach the required sensitivity and, when the first rocket and space based CIB measurements 
were conducted, the results were inconclusive (Matsumoto, Akiba & Murakami 1988, Noda et al 
1992). Following its launch in 1983 the IRAS satellite was the first to conduct all-sky infrared 
measurements of point sources between 12 and 100 fim (Soifcr, Neugebaucr & Houck 1987). It 
revealed that galaxies were efficient infrared emitters, but its design was not optimized for diffuse 
background measurements. The COBE DIRBE instrument, launched in 1990, operated between 
1.25 and 240 jim and was the first to be devoted specifically to CIB measurements (Hauser &: 
Dwek 2001). It led to the first reliable measurements of and limits on the CIB over a wide range 
of infrared bands and literally began the observational CIB era. 
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The status of the CIB measurements in the immediate post-DIRBE era has been reviewed 
extensively by Hauser & Dwck (2001) and a shghtly earher review by Leinert et al (1998) has pro- 
vided a detailed discussion of the issues involved in measuring the extragalactic background light 
(EBL) levels. Since then several important developments happened, most notably the establish- 
ment, through CMB measurements, of the standard cosmological model for the overall evolution of 
the Universe, particularly the existence of the so-called dark energy, and of the standard paradigm 
(ACDM) for structure formation. CMB polarization measurements have identified the epoch of 
re-ionization, presumably due to first collapsed objects in the Universe. This now allows narrowing 
of CIB predictions and puts refined searches for the CIB on a much firmer basis. Also, in the 
CIB context, several new measurements have now been accomplished and several old ones are now 
more firmly confirmed. Scientific exploration is never static and the very recently (August 2003) 
launched NASA Spitzer satellite already led to several important findings for CIB and many more 
are expected to follow as Spitzer observational programs get implemented. New space instruments 
in the IR and sub-mm bands are planned by NASA and ESA which will shed new light on the 
early Universe and the interconnection between the CIB and early luminous systems. This thus 
seems like an opportune moment for a new attempt to review the CIB measurements and their 
interpretation. 

Indeed, following the WMAP and balloon experiments we now know that the Universe is flat 
and dominated by the vacuum energy (cosmological constant) and/or even more exotic quintessence 
field. The recent WMAP results (Bennett et al 2003, Spergel et al 2003) imply that the Universe 
had a large optical depth since last scattering (r ~ 0.2) indicating unexpectedly early star formation 
[z ~ 20). Direct deep studies of the universe at z >2-3 show that much of the luminosity of the 
universe at that time was probably involved in the process of galaxy formation, and the birth of the 
first generations of stars. Many systems, such as elliptical galaxies, already seem well-formed by 
the epoch of the HDF, and the metallicities for many systems are already near solar. Somewhere 
between z ~ 20 — 30 and z ~ 6 reside the systems made up of the first population of stars, the 
so-called zero-metallicity Population III, about which little is known observationally. All these 
galaxies should have left their imprint in the CIB and its structure. 

Observationally, the CIB is difficult to distinguish from the generally brighter foregrounds 
contributed by the local matter within the solar system, and the stars and ISM of the Galaxy. A 
number of investigations have attempted to extract the isotropic component (mean level) of the 
CIB from ground- and satellite-based data as described below. This has in nearly all instances been 
a complicated task due to a lack of detailed knowledge of the absolute brightness levels and spatial 
variations across the sky for the many foregrounds that overlay the CIB signal. It is thus very 
important to understand and estimate to high accuracy the various foreground emissions which 
need to subtracted or removed before uncovering the CIB. 

The outline of this review is as follows: in Sec. 2 we start with a summary of the units, 
the various relevant astronomical magnitude systems and the list of frequent abbreviations used 
throughout the review. Sec. 3 discusses theoretical basis for CIB and its structure and how these 
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are related to information on the early galaxies and stellar systems. Sec. 4 reviews the foregrounds 
that inhibit the CIB measurements and Sec. 5 reviews the status of the current CIB measurements. 
In Sec. 6 we discuss the topics related to CIB measurements and their interpretation such as the 
contribution from ordinary galaxies composed of Population I and II stars, the role of star formation, 
dust, the current status of the present day galaxy luminosity density measurements, and the deep 
galaxy counts from near-IR to far-IR. Sec. 7 is devoted to the signature of Population III in the 
CIB and its structure and the almost unique role CIB plays in uncovering the Population III era. 
We end up with a review of the cTirrcntly planned instruments and space missions that will or can 
have direct bearing on the CIB in Sec. 8 followed by a summary section. 



2. Miscellaneous: definitions, units, magnitudes, abbreviations, etc. 

We start with definitions. The surface brightness of the CIB per unit wavelength will be 
denoted as Ix, per unit frequency as /j/, and per logarithmic wavelength interval F = XIx = uly, 
and we call them all "flux." Throughout the paper B^, will denote the Planck black-body function 
per unit frequency. 

Fluxes of astronomical sources are often measured in narrow band filters. The flux units 
commonly used in astronomy for In are Jy=10~^® W/m^/Hz. The surface brightness of the CIB is 
usually given in units of either MJy/sr or nWrn^^sr""-*^ . The conversion between the two is: 

, nW 3000 MJy 

1 ^ = T7 ^ 1 

m^sr A(/imj sr 

The range of wavelengths used in this review is divided into the following groups: near-infrared 
(NIR) covers 1 //m < A< 5-10 /nm . Mid-infrared (MIR) is defined to lie in 5-10 //m < A< 50-100 
/Ltm range. We call the Far-infrared (FIR) the region corresponding to 50-100 /Ltm < 500 jim. , 
and beyond that will be sometimes referred to as sub-mm. These definitions are neither exact nor 
unique and, although the band ranges cover (largely) different physical processes, this division is 
used for convenience only. 

Throughout the review both Vega magnitudes and AB magnitudes are used. By definition, 
Vega's magnitudes are zero in all filters. However, because of uncertainties in the absolute fiux 
calibration of Vega, magnitudes of this star have been slightly corrected over time. The zero point 
of AB magnitude system is constant flux of 3631 Jy for apparent brightness (Oke & Gunn 1983). 



Table 1. Magnitude system conversion 



Filter 


J 


H 


K 


IRAC-1 


IRAC-2 


IRAC-3 


IRAC-4 


A ( fxm. ) 


1.25 


1.65 


2.2 


3.6 


4.5 


5.8 


8 


mAB-mVega 


0.90 


1.37 


1.84 


2.79 


3.26 


3.73 


4.40 
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Thus in AB magnitude system an object with I^, =constant (flat energy distribution) has the same 
magnitude in all bands, and all colors are zero. For the wavelengths used frequently in the paper 
the conversion is given in Table 1. The numbers were adopted in the J,H)K photometric bands 
from the 2MASS Explanatory Supplement ^ and in the four Spitzer IRAC channels from the Spitzer 
Observer's Manual ^. 

Table 2 lists the abbreviations and acronyms that will appear below. 

3. Theoretical preliminaries 

The near-IR CIB arises mainly from the stellar component of galaxies and probes evolution 
of stellar component of galaxies at early times. The mid- and far-IR CIB originates from dusty 
galaxies reprocessing stellar light and other energetic output. As will be discussed in Sec. 5, the 
recent mutually consistent detections of the near-IR CIB from the COBE/DIRBE and Japan's 
IRTS datasets (see Hauser & Dwek 2001 for review) indicated a surprisingly high amplitude of 
both the CIB fluctuations (Kashlinsky & Odenwald 2000a, Matsumoto ct al 2000,2002) and mean 
levels (Dwek k Arendt 1998, Gorjian k Wright 2000, Wright & Reese 2000, Cambresy et al 2001). 
In the mid-IR firm upper limits on the CIB have been found by various indirect methods, but 
because foreground emission is so high no direct detection has been possible. In the far-IR there 
are mutually consistent detections of the CIB from the DIRBE (Schlegel et al 1998, Hauser et al 
1998) and FIRAS (Puget et al 1996, Fixsen et al 1998) datasets. 

In this section we provide a general mathematical and observational basis for the underlying 

parameters and physics that determine the CIB and its structure. Specific cosmological and galaxy 
evolution models are mentioned only briefly (Sec. 3.4) and we attempt to make the discussion as 
general as possible. 

3.1. Mean level 

In the Friedman-Robertson- Walker Universe with flat geometry and the Robertson- Walker 
metric, ds^ = (?dt^ — (1 + z)^'^[dx^ + x^{dO'^ + s\v? OdcfP')], the comoving volume occupied by a unit 
solid angle in the redshift interval dz is dV/dz = (1 + z)~^dj;^{z)cdt/dz, where dL = x{z)/{l + z) is 
the luminosity distance. Thus the flux density in band A from each galaxy with absolute bolometric 
luminosity L at redshift z is ^^^a'^;^^^) /a( -g)- Here f\d\ is the fraction of the total light emitted 
in the wavelength interval [A; A -|- d\] and the extra factor of (1 -|- z) in the denominator accounts 
for the fact that the flux received in band A comes from a redshifted galaxy. The contribution to 



^http://www. ipac.caltech.edu/2mass/releases/allsky/doc/exp lsup.html 
^http://ssc.spitzer.caltech.edu/documeiits 
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Table 2. Summary of frequently used acronyms and abbreviations. 



2dF 


2 degree Field 


2MASS 


2 Micron All Sky Survey 


AGN 


Active Galactic Nucleus 


CDM 


Cold Dark Matter 


GIB 


Gosmic Infrared Background 


CMB 


Cosmic Microwave Background 


COBE 


COsmic Background Explorer 


DIRBE 


Diffuse InfraRed Background Experiment 


EBL 


Extragalactic Background Light 


ELAIS 


European Large Area ISO Survey 


ESA 


European Space Agency 


FIR 


Far IR 


FIRAS 


Far InfraRed Absolute Spectrometer 


FOV 


Field-Of-View 


FSM 


Faint Source Model 


HDF 


Hubble Deep Field 


IGM 


InterGalactic Medium 


IMF 


Initial Mass Function 


IN 


Instrument Noise 


IPD 


Interplanetary Dust 


IR 


InfraRed 


IRAG 


InfraRed Array Camera 


IRAS 


InfraRed Astronomical Satellite 


IRTS 


InfraRed Telescope in Space 


ISM 


Interstellar Matter 


ISO 


Infrared Space Observatory 


JWST 


James Webb Space Telescope 


MIPS 


Multiband Imaging Photometer System for Spitzer 


MIR 


Mid IR 


NASA 


National Aeronautics & Space Agency 


NGP/SGP 


North/South Galactic Pole 


NEP/SEP 


North/South Ecliptic Pole 


NIR 


Near IR 


NSF 


National Science Foundation 


PAH 


Polycyclic Aromatic Hydrocarbon 


RMS 


Root Mean Square 


SGUBA 


Sub-mm Common-User Bolometer Array 


SDSS 


Sloan Digital Sky Survey 


SED 


Spectral Energy Distribution 


SFR 


Star Formation Rate 


SNAP 


SuperNovae Acceleration Probe 


UDF 


Ultra Deep Field 


WMAP 


Wilkinson Microwave Anisotropy Probe 


ZL 


Zodiacal Light 
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the total CIB flux from the redshift interval dz is given by 

i 

where the sum is taken over all galaxy populations contributing flux in the observer rest-frame 
band at A, and fx characterizes the spectral energy distribution (SED) of galaxy population i. 
Here Rh = cHq^ and the (present-day) luminosity density is given by: 

-^-(0) = 12 [ ^oALv)L,dL, (3) 

where $o.i is the (present-day)luminosity function or the number density of galaxies of morpholog- 
ical type i in the dLy interval at frequency band u. 

It is illustrative to study the redshift dependence of the flux production rate, Eq. (2). At small 
redshifts the factor (1 + z)^'^dt/dz varies little with z, and the flux production rate, dF/dz, is 
governed by the comoving bolometric luminosity density C{z) and the SED of the galaxy emission 
fx- If the luminosity evolution at these redshifts is small, the flux production rate is governed by 
the SED shape. If /a (A) increases toward shorter wavelengths then dF/dz increases with z. For 
Xfx = const, and no luminosity evolution, the rate is roughly constant with small z. At sufficiently 
high redshifts, modest luminosity evolution in Eq. (2) would be ofl"sct by the factor (1 + z)~'^dt/dz, 
so that the flux production rate would be cut off at sufficiently large z. This factor is responsible 
for resolving Olbers' paradox even for a flat SED. 

There arc three broad parameters that determine the possible modes of galaxy evolution at 
higher redshifts: 1) how bright individual galaxies shine (luminosity evolution), 2) their spectral en- 
ergy distribution (SED) which determines how much of the luminosity was emitted at the rest frame 
of the galaxy (K-correction) , and 3) how numerous the galaxies were (number density evolution). 
Consequently, one can separate the z-dependcnce in Ci,{z) into terms due to K-corrcction (/Cjy), 
pTuc luminosity evolution (Si,) and pure number density evolution {Afu) (e.g. Yoshii &; Takahara 
1988), i.e. 

Then the total CIB flux emitted by evolving galaxy populations becomes: 

47r J {l + zf dz ^ ^ 



where dt/dz is given by eq. 20 with the expansion factor R= {I + z) 



-1 



The details of galaxy spectral energy distributions (SED) will be discussed later, but a simple 
analysis can already be made using eq. (2). The SED at rest-frame wavelengths A^IO /xm is 
dominated by stellar emission, with a peak at visible wavelengths and a decrease for A > 0.7 /xm 
. Consequently, assuming no-evolution would mean that most of the predicted J band CIB comes 
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from redshifts z ~ 0.3 — 1, which shifts the visible emission of normal stellar populations to ~ 1 
^m . In the M band at 5 /um , most of the predicted CIB comes from z > 1 — 2. Evolution will 
likely push these redshift ranges toward earlier times. At 10 /um ^A^200 /um , the emission is 
dominated by galactic dust and the situation is reversed, so f\ increases with wavelength roughly 
as A"^ with a ~ 1.5. Hence, the dusty star-burst galaxies observed by IRAS at low redshifts should 
make the dominant contribution to the 10 /xm CIB. In the far-IR, the K-correction is strongly 
negative (Sec. 6.3) and the measured CIB found can have large contributions from high redshifts. 

It is useful to make a simple estimate of the expected CIB flux. Measurements of the galaxy 
luminosity function in both optical (Loveday et al 1992; Blanton et al 2001) and near-IR (Gardner 
et al 1997; Kochanek et al 2003,; Cole et al 2003) indicate approximately that the density of 
bright galaxies is ~ 10~^/i^Mpc~^. Each of these galaxies emits in the rough neighborhood 
of ~ 10^"^ W. The flux-dimensional quantity composed of and the Hubble constant is 

F ~ -^^^L^cHq^ ~ 25 nWm~^sr~^ . This is a crude estimate, but it already shows that in order 
to measure reliably the mean levels of the CIB, foregrounds must be eliminated to well below ~ 10 
nWm~^sr~^ levels. Section 4 shows how difficult a task that is in the infrared bands. 

3.2. CIB anisotropies 

Because of the difficulty of accurately accounting for the contributions of bright foregrounds 
such as Galactic stars, interplanetary and interstellar dust, which must be subtracted from the 
observed sky background (Arendt et al, 1998, Kelsall et al, 1998), Kashlinsky, Mather, Odenwald 
&: Hauser (1996) have proposed to measure the structure of the CIB or its fluctuations spectrum. 
For a relatively conservative set of assumptions about clustering of distant galaxies, fluctuations in 
the brightness of the CIB have a distinct spectral and spatial signal, and these signals can be more 
readily discerned than the actual mean level of the CIB. The most common source of luminosity 
in the universe arises in galaxies, whose clustering properties at the present times are fairly well 
known and are consistent with the ACDM model predictions (Efstathiou, Sutherland & Maddox 
1990, Percival et al 2002, Tegmark et al 2004). The CIB, being produced by clustered matter, must 
have fluctuations that reflect the clustered nature of the underlying sources of luminosity. This 
signature will have an angular correlation function (or angular power spectrum) that distinguishes 
it from local sources of background emission such as zodiacal light emission, and foreground stars in 
the Milky Way. Moreover, distant contributions of CIB will have a different redshift, and therefore 
spectral color, than nearby galaxies and sources of local emission. From galaxy evolution and 
cluster evolution models it is possible to match the predicted slopes for the power spectrum of the 
CIB against the power spectrum of the data. 

On the largest angular scales (the dipole component), there would be an additional source of 
anisotropy due to our peculiar motion with respect to the inertial frame of the Universe. This 
component may be measurable over a certain range of wavelengths and is discussed after the more 
canonical source of CIB anisotropies, the galaxy clustering. 
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3.2.1. CIB anisotropies from galaxy clustering 

Whenever CIB studies encompass relatively small parts of the sky (angular scales 9 < 1 si) 
one can use Cartesian formulation of the Fourier analysis. The fluctuation in the CIB surface 
brightness can be defined as 5F{0) = F{6) — (F), where F = XIx, 6 is the two dimensional 
coordinate on the sky and [F) is the ensemble average. The two-dimensional Fourier transform is 
5F{e) = (27r)-2 / 5Fq exp{-iq ■ e)(fq. 

If the fluctuation field, dF{x), is a random variable, then it can be described by the moments 
of its probability distribution function. The first non-trivial moment is the projected 2-dimensional 
correlation function C{9) = {5F{x + 6)6F{x)). The 2-dimensional power spectrum is P2{q) = 
{\6Fq\'^), where the average is performed over all phases. The correlation function and the power 
spectrum are a pair of 2-dimensional Fourier transforms and for an isotropically distributed signal 
are related by 

Ci9) = ^j^ P2{q)Mq0)qdq, (6) 
/•oo 

P2{q) = 27r / C{9)Joiq9)9d9, (7) 
Jo 

where J„(x) is the n-th order cylindrical Bessel function. If the phases are random, then the 
distribution of the brightness is Gaussian and the correlation function (or its Fourier transform, 
the power spectrum) uniquely describes its statistics. In measurements with a finite beam, the 
intrinsic power spectrum is multiplied by the window function W of the instrument. Conversely, 
for the known beam window function, the power spectrum can be de-convolved by dividing the 
measured power spectrum by the beam window function. 

Another useful and related quantity is the mean square fluctuation within a finite beam of 
angular radius i}, or zero-lag correlation signal, which is related to the power spectrum by 

1 r°° 

C(0) = {{SFf)^ = — P2{q)WTH{q^)qdq 

~ ^q^P2{q)U^n/2^. (8) 

For a top-hat beam the window function is Wth = ['^Jiix)/x]'^ = 0.5 at a; 2± 7r/2 where x = q'd, 
and hence the values of correspond to fluctuations on angular scales of diameter 2± ir/q. 

At small angles < 1 sr, the CIB power spectrum is related to the CIB flux production rate, 
dF/dz, and the evolving 3-D power spectrum of galaxy clustering, Psik) via the Limber equation. 
In the power spectrum formulation it can be written as (e.g. Kashlinsky &; Odenwald 2000) : 



where dA{z) is the angular diameter distance and the integration is over the epoch of the sources 
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contributing to the CIB. Eq. 9 can be rewritten as: 



1 [ ( dl^^V Pz{qd-j^ \ 
cj \ dt d^ 



Mq) = Z I {^] ^"g' (10) 



This is equivalent to: 



where to is the time-length of the period over which the CIB is produced and 

is the fluctuation in number of sources within a volume k^'^ct^. ^ 

To within a factor of order unity, the square of the fractional fluctuation of the CIB on angular 
scale ~ -k/q is = {{61^)^) /I^ ~ I^^q^P2{q)/27r. The meaning of eq. 11 becomes obvious 
if we assume ci/j,/ /di=constant during the lifetime of the emitters to- In this case the fractional 
fluctuation due to clustering of early galaxies becomes: 



a2 



J A\qd^';z)dt (13) 



In other words, the fractional fluctuation on angular scale 7r/g in the CIB is given by the average 
value of the r.m.s. fluctuation from spatial clustering over a cylinder of length cto and diameter 

The Cartesian formulation is equivalent to the spherical sky representation used in the cosmic 
microwave background (CMB) studies on small scales. It was used in some CIB analyses (Haiman 
& Knox 2000, Cooray et al 2003). In that case one expands the flux into spherical harmonics: 
F{e, (/)) = ul^ = J2Zo T!ni=-i <^imYim{0, <P). The correlation function of the CIB, 0(9) = {6F{x) • 
SF{x + 9)), is then C{9) = Z ^CiPi(cos^) with Ci = {\ai\^) = {21 + l)-i ^L-/ and 
Vi denoting the Legendre polynomials (e.g. Peebles 1980). The angular power spectrum, P2{q)-, 
is the two-dimensional Fourier transform of C{9) and for Z S> 1 is related to the multipoles via 
Ci = P2{1 + \)- This follows because at Z S> 1 the Legendre polynomials can be approximated as 
Besscl functions, Vi{cos9) ~ ^o((' + 1/2)^)- The magnitude of the CIB fluctuation on scale tt/Z 
radian for large / is then ~ \/PCi/2t:. In the limit of small angles {I 3> 1) the values of C/'s are 
related to the power spectrum of galaxy clustering and the CIB flux production rate via: 

1 [ fdl^:\^ 1 I ^ + \ 



^The A(fc) defined by eq. 12 and used throughout the review should not be confused with another quantity 



\/k'-^P{k)/2n - sometimes encountered in the hterature under the same symbol. 
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In surveys with arcsec angular resolution it is possible to identify and remove galaxies brighter 
than some limiting magnitude, mn^. Because on average fainter galaxies are at higher z, by improv- 
ing sensitivity and angular resolution, one can isolate contributions to the CIB fluctuations from 
progressively earlier epochs (Kashlinsky et al 2002). The luminosity density, £ = J^^^"^'™) ^{L)dL, 
is then peaked at some particular redshift - at lower z the (brighter than apparent magnitude miim) 
galaxies are removed and at larger redshifts the luminosity density gets dominated by the bright end 
of the luminosity function with a sharp drop-off in the galaxy number density. The following toy 
model is useful in estimating the effect: in the visible to near-IR bands the present day galaxy lumi- 
nosity function is of the Schcchtcr (1976) form $ = <1>^L~-'^(L/L*)~° exp(— L/L^). Measurements of 
the galaxy luminosity function from B to K bands indicate that within the statistical uncertainties 
a ~ 1 (Loveday et al 1992, Gardner et al 1997), leading to C = $*L*(1 — exp[— L(miim)/-^*])- This 
then defines the redshift window in eqs. 10,4 which contributes most to the power spectrum of the 
CIB. Contribution from low z galaxies, for which L(miijn) < is £ ~ 'I>*L(mi;m) oc d!j^{z). In prac- 
tice the typical redshift at which most of the contribution arises can be estimated as L(miijn) L^,. 
If one can further remove galaxies lying in narrow bins. Am, around progressively fainter apparent 
magnitude mum, one can hope to isolate contributions to the CIB by epoch (Kashlinsky 1992, 
Kashlinsky et al 2002). 



3.2.2. Shot noise fluctuations from individual galaxies 

In addition to fluctuations from galaxy clustering there would also be a shot-noise component 
arising from discrete galaxies occasionally entering the beam. The relative amplitude of these shot 

— 1/2 

noise fluctuations will be ~ -^beam 'w^'^^re iVbeam is the average number of galaxies in the beam. 
This component is important in surveys with good angular resolution where A^beam $ a few. 

If galaxies are removed down to some magnitude m, the shot noise contribution from the 
remaining sources to the flux variance, C(0), in measurements with beam of the area Wbeam steradian 
is given by: 

^sn = r F\m)'^dm (15) 

Here F{m) = FolO~°-^"^ is the flux from galaxy of magnitude m and dNg^i/dm is the number of 
galaxies per steradian in the magnitude bin dm. 

Fourier amplitudes of the shot noise are scale-independent and equation 8 implies that the 
shot noise contribution to the power spectrum of the CIB would be given by: 




Because galaxy clustering has power spectrum that increases towards large scales, the shot noise 
component becomes progressively more important at smaller angular scales. 
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3. 2. 3. Cosmic variance for CIB anisotropics 



Any measurement of the angular power spectrum will be affected by the sample or cosmic 
variance in much the same way as the cosmic microwave background measurement (Abbot k, Wise 
1984). This results from the fact that in the best of situations we can only observe Air steradian 
leading to poor sampling of the long wavelength modes. If the power spectrum is determined from 
fraction /gky of the sky by sampling in concentric rings of width A.q in angular wavenumber space, 

— 1/2 

the relative uncertainty on P2{q) will be Nq ' , where Ng oc qAq is the number of ring elements 
in [q; q + Aq] interval. Therefore, the relative uncertainty from cosmic variance in the measured 
power spectrum on scale c^ir/O will be: 



In order to get reliable and independent measurements at a given scale, it is useful to have narrow 
band Aq/q ~ 0.05 — 0.1. Therefore, for reliable measurements on scales up to 6, one has to cover 
an area a few times larger. 



The dipole anisotropy of the CIB arises from our local motion with respect to the inertial 
frame of the Universe, rather than galaxy clustering. It carries important cosmological information 
and its amplitude and wavelength dependence can be predicted in a model-independent way and 
may one day be measurable. 

If all of the dipole anisotropy of the CMB is produced by pecuUar motion of the Sun and the 
Local Group with respect to the inertial frame provided by a distant observer (the last scattering 
surface in the case of the CMB or early epochs, high z, in the case of the CIB), the CIB should 
have dipole anisotropy of the corresponding amplitude and in the same direction. The amplitude 
of the dipole anisotropy can be characterized by the first term, Ci in expanding the sky in spherical 
harmonics. Since Ii./i'^ is an optical constant along the ray's trajectory, the motion of the terrestrial 
observer at speed Vpec with respect to the observed background will produce dipole fluctuation of 
the amplitude SIu/{Iu) = (3 — a^){vpec/c) cos 6 = Cicos9 (e.g. Peebles Sz Wilkinson 1968). Here 
9 is the angle between the line-of-sight and the direction of motion and ai, = dlio.Ii,/dlio.i' is the 
spectral index of the radiation. 

For CMB measurements in the Rayleigh-Jeans part of the CMB spectrum the index is acMB — 
2. Hence, the CIB dipole is related to that of the CMB via Ci^ciB = (3 — ctu) Ci.cmb/?cmb ^cib- 
The CMB dipole is known very accurately to be Ci,cmb = 3.346 ± 0.017 mK (Bennett et al 2003). 
Hence the CIB dipole amplitude is expected to be 
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3.2.4- CIB dipole component 



<^i^CiB,dipoie ^ 1-2 X 10-3 (3 - a^) FciB 



(18) 
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100 1000 

Fig. 1. — The insert shows the mean spectral index of the far-IR CIB according to eq. 24 (solid line) with its approximate 
uncertainty (dotted lines). Main figure shows the expected amplitude of the CIB dipole anisotropy normalized to the CMB 
dipole. The mean CIB dipole according to eq. 24 is shown with solid line and the shaded area denotes its approximate 
uncertainty. Crosses with errors show the dipole expected for the DIRBE measured levels of the CIB at 140 and 240 fim 
(Hauser et al 1998), which are also plotted with crosses in Fig. 9. Thick dashed line shows the latest measurements of the CMB 
dipole (Bennett et al 2003). 



The current measurements of the CIB are discussed at length in Sec. 5. In the far-IR they show 
that the energy spectrum of the CIB has a 'window' where is strongly negative (sec Fig. 9 and 
eq. 24). Fig. 1 shows the spectral index of the CIB derived from the FIRAS based measurements 
of the CIB, eq. 24, (Puget et al 1996, Fixsen et al 1998) and the resultant CIB dipole amplitude 
normalized to the observed CMB dipole. This component may be measurable at wavelengths where 
the CIB has a strongly negative ai, and a sufficiently high flux level. The measurements show that 
ai, is strongly negative for A < 200;um reaching the levels of a,^ ~ — 2 to —4 there. Around ~ 100 
fim the CIB spectrum has ~ —4 and the CIB dipole should be ~ 7 times more sensitive than 
the CMB dipole. The CIB flux at these wavelengths is ~ 10 - 20 nWm-^sr-i and the CIB dipole 
at this range of wavelengths should have a non-negligible amplitude of ~ 0.1 — 0.2 nWm~^sr~^ and 
may be detectable in some future measurements. At longer wavelengths the CIB dipole will be 
difficult to measure due to steep increase in the residual uncertainty from the CMB dipole which 
is also shown in the Figure. 

The CIB dipole anisotropy could be additionally enhanced because the over-density that pro- 
vides our peculiar acceleration also presumably has excess IR luminosity. Indeed, there are persis- 
tent claims of large bulks flows on scales of ~ 150 — 200/7.^"*^ Mpc, whose direction roughly coincides 
with that of the CMB dipole (sec Willick 2000 and references cited therein). Furthermore, the 
dipole of the distribution of rich (Abell) clusters does not converge out to ~ 200^~^Mpc, while its 
direction roughly coincides with the CMB dipole (Scaramella et al. 1991). 
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The measurement of the CIB dipole anisotropy should be possible in the wavelength range 

100-300 /7,m and will be important to provide additional information on the peculiar motions 
in the local part of the Universe and will serve as additional, and perhaps ultimate, test of the 
cosmological nature of any CIB detection at that wavelength. 

3.3. Cosmological paradigms 

CIB levels and structure depend on the history of energy production in the post-recombination 
Universe. The energy production is driven mainly by nucleosynthesis which is related to the history 
of the baryonic component of the Universe. Ultimately, it is the gravity that drives baryon evolution 
and the latter is determined by the evolution of the density inhomogeneities, the nature of dark 
matter and the cosmological parameters. 

With the WMAP (Bennett et al 2003) measurements the structure of the last scattering sur- 
face has been mapped and a firm cosmological model has now emerged: the Universe is flat, 
dominated by the vacuum energy or an exotic quintessence field and is consistent with the in- 
flationary paradigm and a cold-dark-matter (CDM) model. In this paper, we adopt the ACDM 
model (Efstathiou et al 1990) with cosmological parameters from the WMAP and other observa- 
tions: fibaryon = 0.044, /i = 0.71, = 0.3, = 0.7,(78 = 0.84. Following the last scattering at 
z ~ 1000, the Universe entered an era known as the "Dark Ages", which ended with the star for- 
mation which produced Population III stars. The WMAP polarization results (Kogut et al 2003) 
show that the Universe had an optical depth since last scattering of r ^ 0.2, indicating an unex- 
pectedly early epoch of the first star formation (2;^ ~ 20). From the opposite direction in z, optical 
and IR telescopes are now making progress into understanding the luminosity history during the 
most recent epoch of the universe (z < 5), but the period from recombination to the redshift of 
the galaxies in the Hubble Deep Field (HDF) remains largely an unexplored era. CIB offers an 
alternative and powerful tool to probe those epochs. 

The evolution of dark halos, destined to convert baryons into stars, is fixed by the power 
spectrum of the (dark) matter inhomogeneities. The latter are believed to have been imprinted 
during inflationary era and COBE DMR and WMAP observations confirm that it started out with 
the scale invariant spectrum of the Harrison-Zeldovich slope. Within the CDM framework the 
later evolution of the density fluctuations is fairly well understood: during radiation-dominated 
era fluctuations inside the horizon remain frozen, whereas super-horizon modes grow self-similarly. 
When the Universe became matter dominated, all modes grew at equal rate. The epoch of the 
matter radiation equality thus determines the overall shape of the power spectrum with the horizon 
scale at that time (oc f^mattcr^^) being the only scale imprinted. The initial power spectrum is 
modified by the so-called transfer function and, in linear approximation, depends only on the 
cosmological parameters. Various approximations exist for its shape; we chose the approximation 
from Sugiyama (1995) with normalization to the COBE and WMAP CMB anisotropies for the 
numbers that follow. 



-17- 



The amplitude of matter density fluctuations on a given scale ir/k is ~ ^/k^P{k)/2'!r^, where 
P{k) is the power spectrum of density perturbations. Fourier modes evolution in the post-recombination 
Universe is determined by the cosmological parameters and the equation of state, p = wp. Cosmo- 
logical constant, or vacuum energy, would lead to w = —1; various quintessence models generally 
require —1 < w < —1/3 and the "normal" matter dominated Universe requires w = 0. The growth 
of fluctuations, * = 5{z)/S{0) is governed by the differential equation: 



^ + 2-* - 47rGpmatter* 







(19) 



where R = {1 + z) ^ is the expansion factor. This is coupled with the cosmic time - redshift 
relation, or Friedman equation, which for the flat Universe becomes: 



R^ = Hi 



. 1 



(20) 
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Fig. 2. — Left: cosmic time as function of redshift. Middle angular scale subtended by the physical scale l/i~^Kpc vs z; 
and Right the growth factor for linear density fluctuations. Solid lines in all pajiels correspond to = 1 with w = and sets 
of broken lines to flat Universe with Qm = 0.3 and equation of state P = wp(? with w = —1, —2/3, —1/2. 

Fig. 2 shows the various parameters out to z = 20, the redshift of the first star formation 
indicated by the WMAP polarization measurements at large angular scales, which will be useful 
throughout the review: left panel shows the cosmic time from the Big Bang to the given redshift 
z (WMAP analysis suggests the present age is ~ 14 Gyr), middle panel shows the angular scale 
subtended by physical scale lh~^Kpc, and the right panels shows the growth evolution of density 
fluctuations, '^{z). The middle panel also shows the necessity of the GIB studies: e.g. the candidate 
z ^ 7 galaxy (Kneib et al 2004) is estimated to have the total extent of ^ 1 Kpc and is likely to be 
among the largest systems at those epochs. Even if one resolves such compact systems, one would 
have to compromise on simultaneously observing a sufficiently large part of sky to gather robust 
statistics about the abundance and large-scale distribution of (very) high redshift galaxies. Thus 
the bulk of stellar material at high redshifts may be inaccessible to current and even future direct 
studies unless such work is complemented by the measurements of the GIB and its structure. 
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3.4. From cosmological paradigm to galaxies. 

Eq. 19 applies only in the linear regime, where the density contrast 5p = Sp/p « 1. As 
time goes on density fluctuations grow until they become non-linear, turn around, separate from 
the general expansion frame and collapse to form compact objects. For CDM models the typical 
density contrast increases toward smaller scales, so the small scale objects collapse earlier. At 
present epoch the RMS fluctuation in the counts of galaxies is close to unity at rg = 8/i~^Mpc 
(Davis & Peebles 1983) and the power spectrum of galaxy clustering has been accurately measured 
on scales up to ~ 100/i~^Mpc from the 2dF (Percival et al 2002) and SDSS (Tegmark et al 2004) 
surveys and found in good agreement with the ACDM model (Efstathiou et al 1990). The power 
spectrum should be normalized to reproduce at present {z = 0) the RMS fluctuations of erg = 0.84 
over a sphere of radius rg. On non- linear scales gravitational effects would modify the shape and 
amplitude of the mass power spectrum. We modeled its evolution with z using the Peacock & 
Dodds (1996) approximation. Fig. 3 shows the evolution of the mass power spectrum, or the RMS 
density contrast 5rms — s/f^P^i^)]^^^ ■, for ACDM model normalized to WMAP cosmological 
parameters as function of scale. Non-linear scales correspond to ■\/WP^{k)j2T^ ^1 and in the 
Harrison-Zeldovich regime ^/WP^{k)j2^ oc k'^. The total mass contained in a given scale is 
M ~ 1.25 X lO^^(r/l/i~^Mpc)^r2/i~^M0; the baryonic mass would be smaller by fibaryon/^^matter- 
Adopting spherical model for the evolution and collapse of density fluctuations would give that 
any mass that, in linear approximation, reached the density contrast of 5co\ = 1.68 could collapse. 
Therefore r] = 5rms/'^co1 gives the number of standard deviations a given object had to be in order 
to collapse at the given z. CDM and most inflation inspired models predict that the primordial 
density was Gaussian, so the probability that a given mass has collapsed at redshift z is given by 
Pm = erfc(r/(z)/V2). 

Once the fluctuation on a given scale turns around and collapses, gaseous processes will be 
critically important in determining how, when and which compact objects will form. Hoyle (1953) 
was the first to point out the importance of cooling in determining the galaxy masses, Rees &: 
Ostriker (1977) have set it in the cosmological context. Galaxy morphology (roughly speaking, 
ellipticals vs disk galaxies) may be related to either dissipative collapse and fragmentation with 
angular momentum determining the efficiency of fragmentation and star formation (Kashlinsky 
1982) or mergers (Toomre & Toomre 1972). 

A further complication is that the distribution of luminous systems may not trace that of the 
mass, or the so-called biasing. This is likely to be important at early times when stars formed in 
rare regions with fluctuations of the necessary amplitude. Biasing amplifies the 2-point correlation 
function in locations traced by the objects (Kaiser 1984). Various biasing schemes exist (e.g. 
Politzer k Wise 1986, Jensen & Szalay 1986, Kashlinsky 1987, 1991, 1998) which may (or may not) 
at least be good approximations to real situations. 

The era of metal-rich stellar populations composed of Population I and II stars, which we term 
'ordinary', was likely preceded by a (possibly brief) period of the first zero metallicity stars, the 
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Fig. 3. ACDM density field at z = 20, 10, 5 and (thin to thickest) using the Peacock & Dodds (1996) approximation for 

non-lineaj: evolution. Dotted lines corresponds to linear density field and solid lines include non-linear evolution. Filled diamonds 
with errors show the measurements of the present-day power spectrum of galaxy clustering from SDSS survey (Tegmark et al 
2004); assuming no biasing it would coincide with the power spectrum of mass fluctuations. The top axis shows the mass in 
h~^MQ contained within a comoving radius of 1/fe for Q = 0.3. 

so-called Population III. Theoretically this is a special case, much less rooted in observational data, 
but also very important for understanding the later formation and evolution of ordinary galaxy 
populations. These populations and their imprint on the CIB are discussed in a separate sec. 7. 

In order to interpret the measurements of the CIB, one needs to understand or model the details 
related to individual galaxy formation and evolution. Partridge &: Peebles (1967) were the first to 
discuss the levels of the CIB expected from galaxy formation in the Big Bang model. Tinsley (1976, 
1980) pioneered the studies of evolution of stellar populations, their SEDs and metallicity. Bond, 
Carr & Hogan (1980) have discussed a wide range of scenarios under which a measurable CIB can 
be produced, from contribution by primeval galaxies to those of decaying elementary particles; they 
also provided early estimates of the angular power spectrum of the expected CIB anisotropics. Later 
works continued the empirical approach or using a combined N-body and analytical machinery. The 
number of papers relevant to this brief theoretical section by far exceeds the limit of this review 
and in what follows we concentrate on only a few representative works to illustrate the tools used 
and the results obtained. 

• Backward and forward evolution: Lonsdale (1996) divided empirical modeling of CIB into 
the forward and backward evolution classes. In the backward evolution approach one uses the 
measurements of the CIB, galaxy luminosity function, their assumed or modeled SEDs, observed 
galaxy counts, etc to translate galaxy evolution back in time. The forward evolution approach 
takes the opposite route, starting with galaxies at some initial epoch and evolving their properties 
to the present time. Very often, both approaches are needed. 

Galaxy modeling generally consists of the following stages: 1) One specifies galaxy morphology 
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and morphological evolution with cosmic time; 2) One then needs to specify the stellar IMF, which 
may depend on cosmic time and vary among the various galaxy types; 3) Star formation rate must 
be assumed for each morphological type and z; 4) Stellar evolution tracks must be distributed 
according to the IMF and SED for each morphological type reconstructed at each cosmic time; 
5) Stellar emissions depend on the cosmic abundances and chemical evolution of the ISM must be 
accounted for and dust formation treatment is critically important at mid- to far-IR; 6) The various 
parameters described above may vary among galaxies of different mass or luminosity and this must 
be specified also; 7) For CIB fluctuations one must assume the model for the power spectrum of 
galaxy clustering, biasing and their time evolution. Additionally, these parameters depend on the 
total age of the Universe, the redshift of galaxy formation Zf, which may vary for different galaxies, 
etc. The models must be normalized to reproduce as many observational data as possible and are 
constrained by the data on galaxy morphology, metallicity Z (including at high z) , observations of 
galaxy counts at various wavelengths, colors, etc. Despite the many parameters (and uncertainties) 
involved in the above construction, one can arrive at meaningful and fairly accurate limits on galaxy 
formation and evolution from the CIB. 

Yoshii & Takahara (1988) used the stellar evolution models from Arimoto & Yoshii (1987) in 
order to compute the near-IR CIB (and the EBL) as function of the galaxy formation epoch and 
the deceleration parameter. The evolving galaxies were divided into five morphological types and 
their colors and counts were computed for the various cosmological and evolutionary parameters. 
They found that at visible bands most contribution to the EBL comes from late types (Sa to Sd), 
whereas at the NIR bands early types (E/SO) and late types give comparable contributions. Totani 
et al (1997) have extended the Yoshii &: Takahara modeling for various morphological types to 
model the observed evolution of the luminosity density from UV to 1 ^m (Lilly et al 1996). They 
noted that if the proportion of galaxy mixes remains constant from z ~ 1, then a flat Universe 
dominated by a cosmological constant is preferred by the data. 

Franceschini et al (1991) made an early detailed study with predictions of future extragalactic 
surveys from near- to far-IR. They modeled chemical evolution of early and late-type galaxies 
normalized to K galaxy counts and allowed for various evolutionary modes of AGNs. Mid- to far- 
IR emission from dust was modeled with data from IRAS galaxies allowing for realistic distribution 
of dust grain sizes and temperatures and included emission from PAHs. They presented estimates 
of confusion limits from both Galactic stars and high-z galaxies. 

Fall, Chariot & Pei (1996) presented an original way to relate the CIB to other observed 
parameters. They assumed that the dust at each redshift, having the same spatial distribution 
as stars, is traced by the neutral hydrogen. In turn, the neutral hydrogen column density can 
be normalized to reproduce the observed comoving density of HI in damped Ly-a systems out to 
z 2± 4. They assumed a Salpeter type IMF for all galaxies and solved for chemical evolution under 
various approximations. The resultant CIB at long wavelengths was in good agreement with the 
detection from FIRAS data (Puget et al 1996). In the near-IR the predicted CIB is produced 
by stars, and they discuss how warm dust can produce substantial levels of mid-IR CIB without 
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significantly affecting the far-IR part of the CIB spectrum, which was argued to come from the 
z^4: galaxies responsible for the damped Ly-a systems. They note that, for a wide range of the 
IMF, the agreement of the computed cmissivities with estimates from the Canada-France Rcdshift 
Survey or CFRS (Lilly et al 1996) requires that the initial density parameter of material that went 
into stars and dust at was between 10~^h~^ and 8 x 10~^/i~^. 

Jimenez &; Kashlinsky (1999) used synthetic models of stellar populations to analyze the con- 
tribution of normal galaxies to near-IR CIB and its fluctuations. The galaxies were assumed to 
form at fixed Zf and were divided into five morphological types from E to Sd/Irr. Stars in the 
disks were assumed to form with the Scalo IMF and in spheroids with the Salpeter IMF. Late type 
galaxies were modeled to follow the Schmidt law for star formation with the timescale which de- 
pended on the bulge-to-disk ratio. Ellipticals were assumed to form at the fixed high zj in a single 
burst of star formation and passively evolve. Because early type galaxies contribute a significant 
fraction of the near-IR CIB one has to be careful in normalizing to the observed properties of these 
populations; hence the early-type galaxies were normalized to the present-day fundamental plane 
of elliptical galaxies assumed to mimic the metallicity variations along their luminosity sequence. 
They found that, despite their simplicity, the models gave good fits to the observed near-IR galaxy 
counts, the CFRS measurements of the luminosity density evolution (Lilly et al 1996), and gave 
good matches to the observed metallicities of the damped Ly-a systems (Pettini et al 1997) and 
galaxy colors. This then allowed to refine the range of predictions for the CIB and its fluctuations 
in the near-IR. 

In the far-IR the bulk of the CIB comes from dust emission. Beichman &; Helou (1991) 
constructed a model of dust emission, which coupled with the observed IRAS galaxy luminosity 
function and evolutionary assumptions for the various components of galaxy emission, allowed 

(pre-COBE) theoretical estimates of the far-IR CIB. Dwek et al (1998) normalized the mid- to 
far-IR emissions from galaxies to the mean SED of IRAS galaxies and reproduce the far-IR CIB by 
assuming a simple evolution of thereby constructed galaxies. Haiman & Knox (2000) computed the 
angular spectrum of the far-IR CIB fluctuations assuming that the dust temperature is determined 
by the UV radiation, i.e. that both the rate of dust production and its temperature are determined 
by the (measured) SFR (Star Formation Rate). They get that for a wide range of models the 
far-IR CIB fluctuations should be ~ 10% of the mean CIB levels on degree and sub-dcgrcc scales 
consistent with more general arguments of Kashlinsky et al (1996) and Kashlinsky &: Odenwald 
(2000). Haiman and Loeb (1997) discussed the imprint in the CIB produced by smooth dust in IGM 
spread by remnants of the first stars. Knox et al (2001) have further expanded the analysis exploring 
the possibilities of measuring the FIR CIB anisotropies with the Planck space mission. Lagache 
et al (2003) considered a phenomenological model at far-IR to sub-mm wavelengths, constructing 
sample SEDs for starburst and normal galaxies and assuming only the luminosity function evolution 
with z. The model gave successful fits to the data on galaxy counts and rcdshift distribution, the 
far-IR CIB and they used it to give predictions for future observations with Herschel and Planck 
missions. 
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• Semianalytical galaxy formation modeling presents another track for theoretical insight and 

revolves around various numerical codes to generate dark matter evolution for a given hierarchical 
clustering model (usually CDM) and trace the merging history for present galaxy haloes. Prescrip- 
tions for galaxy formation inside the formed dark matter haloes are put in following gas dynamics 
and radiative processes and hydrodynamics with the choices constrained by data on e.g. Tully- 
Fisher (1977) relation, Lyman-a galaxies, galaxy colors, SFR history etc (e.g. Kauffman, White & 
Guiderdoni 1993, Baugh et al 1998, Sommerville, Primack &; Faber 1999). Guiderdoni et al (1998) 
used semi-analytic modeling, normalized to reproduce the observed far-IR/sub-mm part of the GIB, 
to study a plausible range of redshift distributions and faint galaxy counts at these wavelengths for 
the various evolutionary assumptions. 



4. Obstacles to measurement: confusion and foregrounds 
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Fig. 4. — IR emission from the various foregrounds in CIB measurements: zodiacal light at high Ecliptic latitudes (crosses, 
Kelsall et al 1998), Galactic cirrus at high Galactic latitudes (diamonds, Arendt et al 1998), atmospheric emission fluctuations 
at 1" after ~ 8 sec integration (filled circles, http://pegasus.phast.umass.edu/adams/airglowpage.html), and CMB (solid line) 

GIB measurements are very difficult because while the GIB signal is relatively weak, the various 
foregrounds conspire to be fairly bright as one moves to the IR range of wavelengths. Fig. 4 
summarizes the various foregrounds and the remainder of this section discusses them and their 
structure. The figure illustrates the difficulty of eliminating the foreground emission down to 
levels of ~ 10 nWm~^sr~^ in the infrared. While some of the foregrounds may be small at certain 
wavelengths, their sum-total always conspires to be well above the required level of ~ 10 nWm~^sr~^ 
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Foreground emission from the Galaxy and the solar system is the main problem in unveiling 

the expected CIB. At wavelengths less than 10 /im , the dominant foreground after removing 
the zodiacal light model is emission from stars in our Galaxy. At sub-mm wavelengths cosmic 
microwave background (CMB) emission dominates everything else. In ground-based measurements 
atmospheric emission is important. 

Intensity fluctuations of the Galactic foregrounds are perhaps the most difficult to distinguish 
from those of the CIB. Stellar emission may exhibit structure from binaries, clusters and associ- 
ations, and from large scale tidal streams ripped from past and present dwarf galaxy satellites of 
the Milky Way. At long IR wavelengths, stellar emission is minimized by virtue of being far out 
on the Rayleigh-Jeans tail of the stellar spectrum (apart from certain rare classes of dusty stars). 
At near-IR wavelengths stellar emission is important, but with sufficient sensitivity and angular 
resolution most Galactic stellar emission, and related structure, can be resolved and removed. 

4.1. Atmospheric emission 

For ground-based observations, the largest contribution to the sky background comes from the 
atmosphere itself, which at Kg band amounts to 295 MJy/sr. At wavelengths longer than ~ 2.5 
/xm the spectrum of the emission is characteristic of a black-body with a typical atmospheric 
temperature near 250 K. The emission at 1-2.5 fim range is dominated by many intrinsically 
narrow OH lines with some contribution from molecular hydrogen (at 1.27 iim ) and other 
species. The atmospheric seeing is typically ~ 1". On sub-arcminute angular scales fluctuations 
in the atmospheric emission have white noise spectrum both in the spatial and time domains 
and hence scale oc 9~^t~^^^. On larger angular scales atmospheric gradients become important 
(http://pegasus.phast.umass.edu/adams/airglowpage.html). Solid circles in Fig. 4 show typical 
atmospheric fluctuations at 1" during nighttime observations by 2MASS survey after 7.8 sec of 
integration. 

4.2. Galactic stars 

For low resolution experiments (e.g. DIRBE) Galactic stars are a major contributor to fore- 
ground emission at wavelengths ^ 3 fim . In narrow beam observations beam observations they can 
be excised out to fairly faint magnitude. In larger beam measurements, one can use their statistical 
properties for removal of their cumulative emission. For purposes of measuring the CIB, regions of 
the sky within 20° — 30° of the Galactic plane can be ignored. 

In the near-IR Galactic stars, have two useful properties: 1) at the Galactic poles star counts 
have a simple scaling with magnitude given by: 



(21) 
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and 2) outside the Galactic plane (Iboail > 20°) and away from the Galactic center (90° < ^Cai < 
270°) the Galaxy can be approximated as plane-parallel. In a plane-parallel Galaxy in which the 
radial structure can be neglected, the differential counts in the direction x = cosec|6| can be related 
to those at the Galactic pole by (Kashlinsky & Odcnwald 2000) 

dN 3 diV(m-51og^o x) 

-l—\x=X Pole- i^^j 

dm dm 

The right panel of Fig. 5 illustrates the degree of accuracy of the plane-parallel approximation 
for the Galaxy (eq. 22) using DIRBE data (Kashlinsky & Odenwald 2000). Using these approxi- 
mations, once the counts at the Galactic pole are measured, one can evaluate the expected number 
counts in any direction h and then compute the flux and fluctuations in the flux they produce via: 

1"°° dN r°° dN 

i^stars(<m)oc / lO-"-^™ — U^m; a^{x) <x lO'^-^"^ \,dm. (23) 
Jm dm Jm dm 

For B < 0.4 both are dominated by the brightest stars remaining in the field. 

The left panel in Fig. 5 shows compilation of stars counts at 2.2 /xm (K band) in the North 
Galactic Pole (NGP) region from various measurements. The data show that B ~ 0.3; the value of 
B = 0.6 would correspond to homogeneous distribution of sources or stars coming from well within 
the scale height of the Galactic disk. The NGP star counts were observed directly by Elias (1978). 
We show his data at K = 1, 2.5, 3.25 and 8 with N^^^ error bars and our binning of his data. 
Further NGP data were obtained by the 2MASS survey in band, almost identical to the DIRBE 
Band 2, and were kindly provided to us by Tom Jarrett (1998, private communication) who created 
it using the 2MASS point source catalog, now available to the public. The cumulative counts from 
these measurements were shown in Fig. 1 of Beichman (1997) out to Kg > 15, who found that 
they follow dN/dm oc loO-^"^ (cf. his Table 4). Actual 2MASS star counts from a region of 5 square 
degrees centered on the NGP are plotted in Fig. 5. The agreement between the DIRBE counts, 
the Elias (1978) and Jarrett (1998, private communication) data, and the B = 0.3 extrapolation is 
excellent over 15 magnitudes, or six decades in flux. South Galactic Pole counts from Minezaki et 
al.(1998) are also shown. At m^ < 1-5, the counts tend to the slope of B = 0.6 coming from stars 
much closer than the scale height; if B were less than 0.6, the integrated star brightness would 
diverge at the bright end. 

The star counts agree with model predictions. Both Beichman (1996) in his Fig. 1 and Minezaki 
et al. (1998) in their Fig. 1 show that the counts are fitted well by extensions of either the Bahcall 
&; Soneira (1980) or Wainscoat et al. (1992) models. An eyeball fit to their data gives B = 0.3—0.32 
at K = 11. The Wainscoat et al. model at K = 11 shown in Fig. 1 of Minezaki et al. (1998) gives 
\ogd,N/dm. 1.35, whereas continuation of the solid line in the left panel of Fig. 5 to if = 11 
gives logdN/dm = 1.3 if B = 0.3 and 1.4 ii B = 0.33. The agreement between the two slopes 
and normalizations is thus very good. Even the large-beam DIRBE instrument sees far beyond the 
scale height of the bright K band stars. 
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Fig. 5. — Left - Number counts measurements of Galactic stars in K band at the Galactic Poles. DIRBE based counts 
(crosses) are from Kashlinsky & Odenwald (2000) (at fainter magnitudes the DIRBE counts are affected by the confusion noise 
of the large DIRBE beam - sec Fig. 7 in Kashlinsky & Odenwald 2000). Filled triangles arc differential counts from Elias (1978) 
NGP measurements. Filled circles are differential NGP counts from 2MASS (Jarrett 1998). Open triangles are cumulative 
2MASS counts from Beichman (1998) multiplied by 0.3 In 10 to convert to differential counts for dN/dm oc 10""^'". Filled 
diamonds with error bars are South Galactic Pole counts from Fig. 1 of Minezaki et al. (1998). Solid line shows the B = 0.3 
fit. 

Right - DIRBE star counts at 2.2 ;um plotted in coordinates where a plane-parallel Galaxy would collapse on a single line; 
X = cosec|f)|. Data are from the ~ 20 o x20° DIRBE patches with |fe| > 20° and 90° < Z < 270°, with Poisson errors shown. 
For the DIRBE beam the confusion noise aifects counts at K > 5.5. Lines show the model Galaxy: solid is mean dN/dm 
and dashes are the ±l-sigma spread. 



4.3. Zodiacal emission 



Zodiacal emission from interplanetary dust (IPD) is the brightest foreground at most IR wave- 
lengths over most of the sky. There are some structures in this emission associated with particular 
asteroid families, comets and comet trails, and an earth-resonant ring, but these structures tend to 
be confined to low ecliptic latitudes or otherwise localized (Reach et al 1995). The main IPD cloud 
is inclined at ~ 2° with respect to the Ecliptic and is generally modeled with a smooth density 
distribution that also contains additional features of a circumsolar ring, a density enhancement in 
the Earth's wake and the dust bands at several AU (Kelsall et al 1998). Combining DIRBE and 
FIRAS observations. Fixsen &: Dwek (2002) find a sharp break in the dust distribution at radius 
of ~ 30 fim. and that the zodiacal energy spectrum beyond ~ 150 /nm can be fitted with a single 
black-body with emissivity and temperature of 240K. Observationally, intensity fluctuations 
of the main IPD cloud have been limited to < 0.2% at 25 fim (Abraham, Leinert, & Lemke 1997). 
Because the Earth is moving with respect to (orbiting within) the IPD cloud, the zodiacal light 
varies over time. Likewise, any zodiacal light fluctuations will not remain fixed in celestial coor- 
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dinates. Therefore repeated observations of a field on timescales of weeks to months should be 
able to distinguish and reject any zodiacal light fluctuations from the invariant Galactic and CIB 
fluctuations. 

4.4. Galactic cirrus 

IR emission from the ISM (cirrus) is intrinsically diffuse and cannot be resolved. Cirrus 
emission is known to extend to wavelengths as short as 3 nm. Statistically, the structure of the 
cirrus emission can be modeled with power-law distributions, P2{q) oc q"' (Gautier et al 1992), 
and has the power index n ~ —2, —2.5 (Wright 2000, Kashlinsky & Odenwald 2000, Ingalls et al 
2004). This power index lead to almost scale- independent fluctuations in cirrus emission which 
are typically ^ 10^^ of the mean flux level. Using the mean cirrus spectrum, measurements made 
in the far-IR can be scaled to 3.5 fim, providing estimates for the fluctuation contribution from 
cirrus. The extrapolation to shorter wavelengths is highly uncertain, because cirrus (diffuse ISM) 
emission has not been detected at these wavelengths, and the effects of extinction may become 
more significant than those of emission. 

4.5. Cosmic microwave background 

Ironically, the cosmic microwave background (CMB) has to be put in the category of "fore- 
grounds" when searching for the CIB. CMB is a relic from the Big Bang and has a strict black-body 
energy spectrum correspond to Tcmb = 2.725 it O.OOIK with very small, if any, deviations from 
the black-body spectrum (Mather ct al 1993; Fixscn ct al 1996: Fixscn & Mather 2002). CMB is 
also very homogeneous with the largest fractional spatial deviation being ~ 10~^ on dipole scales; 
the other angular moments are some two orders of magnitude smaller (Bennett et al 2003). Its 
brightness, shown with solid line in Fig. 4, overwhelms all other emissions at wavelengths longer 
than ~ 500 /xm and is even a non- negligible source of noise in many everyday electric appliances. 
Because its energy spectrum is known so accurately and its spatial structure is so small, the CMB 
can be subtracted down to levels <C 1 nWm~^sr~^ . 

5. Current CIB measurements 

This section discusses the current observational status of both the mean CIB levels and its 
structure (fluctuations) . 

Observationally, the CIB is difficult to distinguish from the generally brighter foregrounds 
contributed by the local matter within the Solar system, and the stars and ISM of the Galaxy. A 
number of investigations have attempted to extract the isotropic component (mean level) of the CIB 
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from ground- and satellite-based data. These analyses of the COBE data have revealed the CIB at 
far-IR wavelengths A > 100 //m (Hauser et al 1998), and probably at near-IR wavelengths from 1 
- 3 /im (Dwek &: Arendt 1998, Gorjian et al 2001, Cambresy et al 2001) with additional support 
from analysis of data from IRTS (Matsumoto et al 2000, 2003). However, it would probably be 
fair to say that none of the reported detections of the isotropic CIB are very robust (especially in 
the NIR), because all are dominated by the systematic uncertainties associated with the modeling 
and removal of the strong foreground emission of the zodiacal light and Galactic stars and ISM. 
Furthermore, the near-IR colors of the mean CIB do not differ greatly from those of the foregrounds, 
thus limiting the use of spectral information in distinguishing the true CIB from residual Galactic 
or solar system emission. 

In order to avoid the difficulty of exactly accounting for the contributions of these bright 
foregrounds in direct CIB measurements, (and the difficulty in detecting all the contributing sources 
individually), Kashlinsky, Mather, Odenwald and Hauser (1996) have proposed measuring the 
structure or anisotropy of the CIB via its angular power spectrum. They noted that for a relatively 
conservative set of assumptions about clustering of distant galaxies, fluctuations in the brightness 
of the CIB have a distinct spectral and spatial signal, and these signals can be perhaps more readily 
discerned than the actual mean level of the CIB. 

Figure 6 shows the filters and the IR wavelengths covered by the measurements from the 
instruments described below. 




Fig. 6. Filter responses for the various instruments discussed in the paper. The curves were normalized to a unity peak. 

DIRBE, FIRAS, 2MASS, IRTS, Spitzer (IRAC and MIPS) filters correspond to the colors shown in the upper left corner. The 
multiple ISO filters filters are not shown for clarity. 
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5.1. COBE DIRBE 

The primary mission of the DIRBE instrument on board COBE satehite was to find, or set very 
strict Hmits on, the CIB contribution from near- to far-IR. The DIRBE was described by Boggess 
et al. (1992) and is a 10-band, photometer system covering the wavelength range from 1.25 to 240 
/xm with an angular resolution of 0.7°. The 1.25. 2.2 and 3.5 fim DIRBE bands are similar, although 
not identical, to the ground based J, K and L bands. It was designed to achieve stray light rejection 
of less than 1 nWm^^sr^^ and a low absolute brightness calibration uncertainty. Accurate beam 
profile and beam response maps were obtained in fiight by observing multiple transits of bright 
point sources. All spectral channels simultaneously observed the same field-of-view on the sky. Its 
cryogenic stage lasted from November 1989 to September 1990 collecting 41 weeks of data. 

Absolute photometry was obtained through frequent (32 Hz) chopping between the sky and an 
internal zero flux surface maintained at temperatures below 2 K. System response was monitored 
every 20 minutes throughout the 10 month mission by observing an internal thermal reference 
source. The DIRBE was designed to achieve a zero-point to the photometric scale with an uncer- 
tainty below vl{v) = 10^^ Wm^^sr^^ at all wavelengths; tests conducted prior to launch and in 
flight indicate that this goal was indeed reached. 

Absolute calibration of the DIRBE photometric scale, and long-term photometric consistency, 
was established by monitoring a network of stable celestial sources. The primary calibrator for the 
J, K and L-bands was the star Sirius. Although DIRBE is a broad-band instrument, the measured 
in-band intensities are reported as spectral intensities at the nominal effective wavelengths of each 
band. Sources that have an SED significantly different from a flat would require color 

corrections. In the J, K and L-bands, these color corrections are of the order of a few percent or 
less over the temperature range from 1500 to 20,000 K. Long-term stability of the instrumental 
photometric scale at each wavelength was better than 1%. The stability of the band-to-band 
intensity ratios, or colors, was 1.4%. The absolute photometric uncertainty in the DIRBE J, K 
and L-bands was 4%. 

The calibrated DIRBE observations were binned into 20 arcmin pixels of approximately equal 
area for each of the forty one weeks and are available from (http://lambda/product/cobe/dirbe_overview.cfm). 
Each week covered about half the sky and four months of data gave one full sky coverage. 

In a real tour d'force the background and foreground results of the combined 41 week maps of 
DIRBE sky observations have been analyzed by the DIRBE team (Hauser et al 1998, Arendt et al 

1998, Kelsall et al 1998, Dwek et al 1998). From these maps, the models of contributions from the 
interplanetary dust (IPD) cloud and the Galaxy, both its interstellar dust and stellar components, 
were subtracted for each direction. 

The DIRBE IPD model is described in Kelsall et al (1998). It represents the zodiacal sky 
brightness as the integral along the line of sight of the dust emissivity times the 3-D dust density 
distribution function. The emissivity arises from both thermal emission and scattering; the former 
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dominates at A> 5/xm and the latter at shorter wavelengths. The thermal emission at each location 
assumes a single dust temperature for all IPD cloud components and is assumed to drop as a power 
law of the distance to the Sun. The density model included contributions from a smooth cloud, three 
pairs of asteroidal dust bands and a circumsolar dust ring. Parameters of the analytical functions 
that enter the model were determined by matching the observed temporal brightness variations in 
various directions on the sky. Once the optimal set of model parameters was determined, the IPD 
model was integrated along the line of sight at the mean time of each week of observations. The 
calculated IPD brightness map was then subtracted from each week of the DIRBE weekly maps. 
Because of the brightness of the zodiacal emission, the residual sky maps are dominated by the 
uncertainties in the IPD model at between 12 and 60 ^m . 

The Galactic emission was modeled by Arendt et al (1998) as made up of emission of bright 
and faint Galactic stars and discrete sources and the interstellar medium, or cirrus cloud emission. 
It was removed from the combined the maps after the zodiacal component removal. Bright(er) 
sources above a wavelength dependent threshold were removed from the maps at all bands. The 
contribution to the integrated light from faint sources was removed using a statistical source count 
Faint Source Model (FSM - Arendt et al 1998) which is based upon Wainscoat et al (1992). Prior 
to subtracting the FSM, bright stars with flux greater than 15 Jy (corresponding to K~ 4) have 
been subtracted directly. Stellar contribution is negligible at mid- and far-IR DIRBE bands. The 
model for the ISM emission was taken to be a product of the standard spatial template times 
a single spectral factor at each wavelength. At 100 /xm the spatial template was a map of HI 
emission constructed from three different 21 cm HI surveys (see Arendt et al 1998 and references 
cited therein); at other wavelengths the 100 /xm data were used. 

Hauser et al (1998) presented the limits on the CIB after removal of the zodiacal light model, 
the FSM and the Galactic cirrus contributions from the combined DIRBE maps. In the residual 
maps, they selected for the analysis "high-quality" regions that are located high Galactic and 
ecliptic latitudes, are free of foreground removal artifacts and cover at least 2% of the sky. The 
smallest and best region contains some 8,140 pixels in both northern and southern hemispheres 
allowing to test for isotropy of the residual signal. Their criteria for a CIB detection were that 
the residuals had to be in excess of 3-fj, a being both systematic and statistical uncertainty, and 
be isotropic in the high-quality region. They reported firm detections at 140 and 240 ^m and 
upper limits at shorter wavelengths. The 100 ^m residual they find at the level of 21.9 it 6.1 
nWm~^sr~^ also has high statistical significance, but is not isotropic. The dominant errors in the 
DIRBE CIB measurements are the systematic errors associated with the foreground subtraction. 
Dwek et al (1998) show in detail that the isotropic signal detected at 140 and 240 /xm is unlikely 
to arise from some Solar system and Galactic sources. Adopting the DIRBE detections at 140 and 
240 /im , Dwek et al (1998) claim a lower limit of 5nWm^^sr^^ at 100 ^m. assuming that the 
CIB emission arises from the coldest possible dust fitted to the 140, 240 jjcm detections. 

Figure 7 shows the map of the 240 /im residual emission from Hauser et al (1998); it is plotted 
here in Galactic coordinates (cf. Arendt et al 1998). The produced map corresponds to the original 
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240 /um map from which the zodiacal hght was removed along with a two-component ISM model, 
corresponding to a weighted subtraction of the DIRBE 140 and 100 /im maps. The bright sources 
near the Galactic plane are clearly visible as are the Small and Large Magellanic clouds. The 
remaining emission is isotropic and was identified with the CIB. The 140 /im map also passed the 
isotropy limits, although is significantly more anisotropic outside the high-quality region used for 
identifying the CIB. 



Fig. 7. — DIRBE map of the CIB at 240 shown in Galactic coordinates with the South Galactic Pole pointing toward 
the bottom of the figure. 

Chronologically, the DIRBE discovery of the far-IR CIB was preceded by a few months by an 
independent analysis of Schlegel et al (1998), who have combined the well calibrated DIRBE data 
that has relatively poor resolution (~ 0.5°) with the IRAS data, which have much better angular 
resolution (~ 5') but poor absolute calibration. The combined maps at 100 yum provides a very 
accurate to-date template of the diffuse Galactic emission and dust distribution. To construct the 
map, they have removed 1) zodiacal emission from the DIRBE maps, 2) striping artifacts from 
the IRAS data, 3) confirmed point sources, and then combined the maps preserving the DIRBE 
calibration and IRAS resolution. The ratio of the 100 /im to 240 /im emission gave the dust 
temperature leading to estimate of column density of the radiating dust. The 25 fim channel data 
was used to model the zodiacal light emission from the interplanetary dust (IPD) cloud. Fitting 
black-body functions from 12 to 60 /im they find variations in the IPD temperature of less than 
10 %. As the zodiacal light contributes ca. 300 nWm^^sr"^ at 100 fim the temperature variations 
of the cloud lead to absolute errors of ~ 30 nWm~^sr^^ if the 25 /im template is used linearly in 
subtracting from the longer wavelength maps. After accounting for the (non-linear) corrections for 
the IPD temperature variations in subtracting the zodiacal light emission, the CIB was determined 
as the zero offset in cross-correlation between the residual "de-zodied" emission with the Galactic 
HI column density from the Leiden-Dwingeloo 21 cm survey (Hartmann & Burton 1997). The 
analysis has been applied to the 100, 140 and 240 /im maps restricted to lie outside the ecliptic 
plane leading to null-detection at 100 /im and robust CIB values at the two longer bands for 
|/3|>20°. 
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Finkbeiner, Schlegel & Davis (1998) extended the analysis to the DIRBE data at 60 and 100 
/v,m . They apphed two different methods to remove the zodiacal light component: The first of 
these methods used the North-South asymmetry with a 1 year period of the zodiacal component and 
concentrated on the data within 5° of the ecliptic poles, where effects dependent on solar elongation 
cancel out. The method further assumed that the zodiacal light North- South asymmetry is the 
same in all bands. Cirrus was taken from the 100 //m map produced in the Schlegel et al (1998) 
analysis. They decomposed each of the North and South datascts into time-dependent (zodi) and 
tinic-indcpcndcnt (cirrus, CIB etc) components and obtain robust solutions for the latter. In the 
second method they remove zodiacal component by using its dependence on the ecliptic latitude in 
each of the DIRBE's 41 weeks of data with solar elongation 90°. A statistic was then constructed 
that ratios the flux in the various diagonal directions to that toward the ecliptic poles and whose 
advantage is that it has a known dependence on the ecliptic latitude (oc csc|/3|) with a negligible 
annual variation. The CIB flux was determined as a term at csc|/3| = 0. Both methods gave 
consistent and robust, if high compared to other measurements, levels of the CIB at 60 and 100 
/xm . 

Because of the difficulty in measuring the isotropic component of the CIB, Kashlinsky, Mather, 
Odcnwald & Hauser (1996) have proposed to probe the CIB structure instead. For a reasonable set 
of assumptions about the underlying distribution of emitters/galaxies, the fluctuations in the CIB 
can provide important information about the rate of production of the CIB and luminous activity 
in the early Universe. The method leads to particularly strong constraints at mid- to far-IR where 
the foregrounds are very smooth (Kashlinsky, Mather & Odcnwald 1996). The analysis, however, 
could not have been applied to the DIRBE's longest channels (140 and 240 fim ), which required 
the use of bolometers instead of the photoconductors, resulting in much larger noise than at shorter 
wavelengths. 

Kashlinsky & Odenwald (2000) have applied the method to the final DIRBE data release 
leading to detection of the near-IR CIB anisotropies. They started with the final DIRBE sky maps 
from which the zodiacal model from Kelsall et al (1998) was subtracted. The sky was divided into 
(384) patches of 32 pixels (~ 10°) on the side which were individually clipped of Galactic stars and 
other point sources by various procedures all of which led to the same final results. The RMS flux 
fluctuation, (5Frms = \/ C{0), or fluctuation on the DIRBE beam scale (~ 0.6°) was computed in 
each of the clipped patches. The residual fluctuation is the sum of contributions from the remaining 
Galactic stars, CIB, instrument noise, etc. When the data in the near-IR 1.25 to 5 //m were selected 
from high Galactic latitude (l^cail > 20°) and away from the Galactic center (90° < ^Cai < 270°), 
they showed a very tight correlation between cosec(|6Gai|) and C(0). That NIR fluctuation was 
explained in terms of the observed Galactic star counts in the plane parallel Galaxy, such as given 
by eqs. 21,22, shifted by a constant (isotropic) component. The latter presumably arises from CIB 
fluctuations as well as instrumental noise. Extrapolation of the correlation to cosec(|6Gal|) = 0, or 
zero contribution from Galactic stars, gave a robust value of the near-IR CIB anisotropies, which 
was demonstrated to be independent of the clipping methods, sky maps processing and other effects. 
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The instrument noise was evaluated in Kashlinsky et al (1996), Kashlinsky 8z Odenwald (2000) and 
Hauser et al (1998) and was shown to be significantly lower than the detected near-IR signal at 
wavelengths from 1.25 to 5 fim . At mid- and far-IR wavelengths the results of the analysis led to 
slightly tighter upper limits than in the previous analysis (Kashlinsky, Mather & Odenwald 1996) 
which used an intermediate zodiacal light model from the DIRBE team. 

Dwek & Arendt (1998) have presented another original method to isolate the isotropic com- 
ponent of the near-IR GIB by using prior information of the flux in one of the bands. They noted 
that the observed galaxy contribution to the K band CIB pretty much saturates at 2±9 nWm~^sr~^ 
. They started with the DIRBE maps at 1.25, 2.2, 3.5 and 4.9 fxm from which they subtracted 
the zodiacal light model from Kelsall et al (1998). Prom DIRBE near-IR maps the regions were 
selected with small errors in the subtraction of the foreground emission (Hauser et al 1998, Arendt 
et al 1998). Small cirrus emission at 3.5 and 4.9 //m was removed by scaling the average 100 //m 
intensity in DIRBE maps as in Arendt et al (1998). The 2.2 /Ltm DIRBE map minus the assumed 
CIB at that band was taken as a template of the Galactic stellar emission. That template was 
then cross-correlated with the maps at the other wavelengths to give a positive and statistically 
significant residuals (zero intercepts) STimmarizcd in Table 3. This technique was pushed farther 
by Arendt &: Dwek (2003) attaining a more isotropic removal of Galactic foregrounds, but at the 
expense of increased uncertainty in the absolute zero point. 

In order to reduce the uncertainty in the Galactic stars contribution to the NIR DIRBE maps, 
Gorjian et al (2000) have imaged a high Galactic latitude region of the sky of 2° on the side in K and 
L bands with arcsec resolution to K~ 9 and L ~ 8 magnitudes reducing the stellar contribution 
to confusion noise a factor of ~ 16 below that of DIRBE. They subtracted a modified zodiacal 
light model proposed by Wright (1997) which comes in two flavors: a "weak" zodi model which 
requires that the DIRBE high Galactic latitude sky at 25 /im be isotropic in addition to its non-ZL 
components remaining constant in time, and a "strong" zodi model which requires that the former 
component in addition is negligibly small. This reduced the residual 25 jim intensity by a factor 
of 7 compared to the DIRBE model (Kelsall et al 1998). They then selected 17 DIRBE pixels 
where overlap with star position is sufficient for robust statistics in subtracting the Galactic stellar 
component. The CIB was evaluated after subtracting the zodiacal components, the measured flux 
from the Galactic stars and the model flux from fainter and unobserved stars. Wright &: Reese 
(2000) used the above ZL model subtraction and the Arendt et al (1998) method for removing the 
ISM contribution to the residual emission. They found a constant flux offset in the flux histograms 
at 2.2 and 3.5 /xm compared to the Wainscoat et al (1992) model after testing the latter with 2MASS 
observations in several selected fields. The flux offset is identified as the CIB and is consistent with 
the Gorjian et al (2000) results. Wright (2000) used the 2MASS catalog stars in order to subtract 
Galactic star contribution to the DIRBE maps flux in selected regions near the Galactic poles. 
After the subtraction he finds a statistically significant flux excess at 2.2 /Ltm that is in agreement 
with the Gorjian et al (2000) and Wright & Reese (2000) results. 

Cambresy et al (2001) used the 1,400 deg^ of the 2MASS data from standard exposures (7.8 
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sec) in J, H, K^. They degraded the 2MASS maps (2") to 5' resolution. After subtracting the 

zodiacal light from DIRBE weekly maps using the zodiacal subtracted mission average maps, the 
residual intensity maps were averaged together. The subset of the 2MASS that was used to account 
for Galactic contribution to the DIRBE intensities came from high Galactic latitude (|6Gai| > 40°) 
and high ecliptic latitude {\Peci\ ^ 40°) regions of low cirrus emission for which a complete 2MASS 
coverage existed in each DIRBE pixel. After removing bright stars, which contaminate more than 
one DIRBE pixel, the final area for analysis was 1,040 dcg2. The intensities in the DIRBE sky were 
cross correlated with the observed fluxes from 2MASS stars. The zero-intercept of this correlation 
minus the (small) contribution from faint stars, which was computed from Galaxy modeling, is the 
GIB flux. 

Table 3 summarizes the results of the GIB measurements based on the DIRBE data. 

5.2. COBE FIRAS 

Far-Infrared Absolute Spectrometer (FIRAS) on board GOBE was designed for measuring the 
energy spectrum of the GMB and the far-IR GIB. It is a four-port Michelson interferometer. At 
each port (Left and Right) a dichroic filter split the beam into Low (30-660 GHz) and High (600- 
2880 GHz) frequency beams producing spectra with a resolution of 4.2 GHz or 16.9 GHz. The 
FIRAS data have angular resolution of 7 degrees, with pixel size of 2.6 degrees. FIRAS is described 
in Boggess et al (1992) and Mather et al (1993). The FIRAS Pass 4 data consist of spectra between 
104 ^m and 5000 /im in each of 6,063 out of the 6,144 pixels of the sky (81 pixels that have no 
data are omitted). They were calibrated using the method described in Fixsen et al. (1994, 1996). 

The Pass 4 FIRAS data have approximately half of the noise of the previous data releases 
(Pass 3). This is partly due to combining all FIRAS frequency bins and scan modes and partly 
due to improved understanding of the systematic errors. Particularly important was the reduction 
in the systematic errors since these have a complicated spectrum in the spatial domain. Since the 
fluctuations of the far-IR GIB may themselves look like noise with some a priori uncertain spatial 
power spectrum, finding a noise floor leaves one uncertain whether it is the fluctuations in the 
far-IR GIB or systematic effects in the instrument. The Pass 4 FIRAS data have many of the 
known systematics removed, which enables a direct and substantial lowering of the noise spectrum. 

Puget et al (1996) made a first cut at detecting the far-IR GIB from the Pass 3 FIRAS data. 
They used the particular spatial and spectral distribution of the different far-IR components of 
the sky emission in subtracting them from the Pass 3 FIRAS maps. First, DIRBE maps were 
degraded to the FIRAS resolution in order to map and subtract the zodiacal emission. In the 
DIRBE 25 jum channel the emission is dominated by the zodiacal light and they noted the average 
brightness ratio 5,y(100/xm)/Si,(25/Ltm) = 0.167 and used the FIRAS observations of the spectral 
dependence of the zodiacal light, Bi^^^odi (Reach et al 1995), to model the zodiacal light 

emission along each line-of-sight as 0.167(A/100/im)^^i?i,(100/Lim). The GMB along with its dipole 



Table 3. SUMMARY OF DIRBE MEASUREMENTS (nWm-^sr-^ ) 

Ref 1.25 //m 2.2 /im 3.5 //m 4.9 //m 12 /xm 25 /im 60 /im 100 /um 140 /um 240 /um 



1 32 ± 13 17 ± 4 

2 < 75 < 39 < 23 < 41 < 468 < 504 < 75 < 34 25.0 ± 6.9 13.6 ± 2.5 

3 26.9+2.3F2.2;^in - 9.9+0.3F2.2^m 23.3+0. lF2.2j^m 

±20.9 ±2.9 ±6.4 

28.1 24.6 
±1.8 ±7 ±2.5 ±8 

5 22.4 ±6 11.0 ±0.3 

^ 23.1 ±5.9 12.4 ±3.2 

^ 20.2 ±6.3 

8 54.0 ± 16.8 27.8 ± 6.7 



yc(o) 



15.5l?;^ 5.911;? 2.4t0| 2.ot°:f :£ 1- :£0.5 <0.7 



Note. — : Refs: ^ Schlegel, Finkbeiner & Davis (1998), ^ Hauser et al (1998), ^ Dwek & Arendt (1998) - ^2.2^111 is the 
flux at 2.2 fj,m in nWm~^sr~l , * Finkbeiner, Davis & Schlcgcl (2000) — * random and systematic errors shown, ^'^'^ Gorjian, 
Wright & Chary (2000), Wright & Rccsc (2000), Wright (2000), * Cambrcsy ct al (2001), Kashlinsky, Mather & Odenwald 
(1996), Kashlinsky & Odenwald (2000) - * 92% confidence level uncertainties are shown. 
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anisotropy was then subtracted from the raw FIRAS spectrum at each pixel. The Leiden-Dwingeloo 
HI survey (Hartmann & Burton 1997) degraded to the FIRAS resolution was used to subtract the 
Galactic interstellar dust emission which at far-IR correlates well with the neutral hydrogen emission 
(Boulanger et al 1996). At high Galactic latitude (l^cail > 20°) and low HI column density lines- 
of-sight the residual was positive, approximately isotropic, and homogeneous with amplitude of 
~ 3.4(A/400/xm)~^ nWm~^sr~^ in the 400-1000 /xm range. 

Fixsen et al (1998) used Pass 4 FIRAS data with significantly lower noise and systematic errors. 
From the FIRAS data they subtracted the DIRBE zodiacal model (Kelsall et al 1998) extrapolated 
to FIRAS bands and the low-frequency FIRAS CMB model with dipole anisotropy. In order to 
separate the Galactic emission they used three independent methods: 1) They noted that after 
subtracting CMB and zodiacal model, the residual emission has approximately identical spectral 
energy distribution across the sky. This spectral template was to model emission of the Galaxy 
for further subtraction. Using all the FIRAS channels and assuming a non-negative prior for the 
CIB, they find a statistically significant residual. 2) HI (21 cm) and CII (158 /um ) Galactic lines 
were used together with a quadratic fit to the Galactic HI map degraded to FIRAS resolution from 
the AT&T survey by Stark ct al (1992) to subtract Galactic IR emission. After the subtraction, 
the darkest 1/3 of the sky, where the HI data was sufficiently accurate, showed a statistically 
significant CIB residual at longer wavelengths. 3) Finally, they constructed Galactic templates 
from the DIRBE 140 /xm and 240 /xm maps with zodiacal emission subtracted together with 
the levels of the CIB detected by DIRBE at these bands. The template was degraded to FIRAS 
resolution. The CIB was obtained by extrapolating the the correlations with the DIRBE templates 
to the DIRBE measurements at 140 fim and 240 fim . All three methods gave consistent results. 

An attempt to study far-IR CIB fluctuations has been undertaken by Burigana & Popa (1998) 
on the FIRAS Pass 3 data between 200 and 1,000 /xm . From the data they subtracted CMB 
monopolc and dipole components, a two-temperature dust model and a zodiacal light model using 
the DIRBE 25 /xm channel and the v"^ frequency scaling. The analysis resulted in upper limits 
on the CIB fluctuations on scales greater than the FIRAS beam (~ 7°). 

5.3. IRTS 

Japan's InfraRed Telescope in Space (IRTS) was launched in March 1995. The Near InfraRed 
Spectrometer (NIRS) on board IRTS was designed to obtain spectrum of the diffuse background 

emissions (Noda et al 1994). It covers the wavelengths from 1.4 /xm to 4 /xm in 24 independent 
bands with spectral resolution of 0.13 /xm . IRTS observations covered ca. 7% of the sky in 30 
days (Murakami et al 1996). The resolution of IRTS is ~ 8', significantly better than DIRBE and 
allows better removal of faint stars. 

The results from the IRTS CIB analysis have been presented in Matsumoto (2000, 2003). 
A narrow strip at high Galactic latitude was chosen for the CIB analysis in order to reduce the 
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contribution from Galactic stars. Contribution of Galactic stars not resolved by the IRTS beam was 
estimated from the point source sky model of the Galaxy by Cohen (1997). They used the DIRBE 
physical model of the interplanetary dust (Kelsall ct al 1998) to remove the zodiacal component 
corresponding to the NIRS bands and the observed sky directions. There was a good correlation 
between the star-subtracted data and the modeled zodiacal component. The isotropic component 
was determined as the intercept of that correlation at zero modeled zodiacal brightness. The 
brightness of the isotropic emission found in this analysis was about 20 % of the sky brightness 
with an in-band energy of ~ 30 nWm~^sr~^ and was estimated to be a factor of ~ 2.5 greater than 
the integrated light of faint stars. 

In order to determine the spatial spectrum of the fluctuations, Matsumoto et al (2000, 2003) 
have added the NIRS short wavelength bands (with smaller noise) from 1.43 /um to 2.14 ^m 
. They then calculated the correlation function of the emission, C{6), which was converted into 
the power spectrum by Fourier transformation. The RMS fluctuations of the isotropic component 
was, after subtracting read-out noise and stellar fluctuations, consistent with the Kashlinsky & 
Odenwald (2000) results from the DIRBE sky analysis. The color of fluctuations was essentially 
the same as that of the isotropic component, consistent with both having the same origin. 

5.4. ISO 

Infrared Space Observatory (ISO) was operated by the European Space Agency between 
November 1995 and April 1998. Its infrared camera (ISOCAM) had a 32 x 32 pixel array and 
two IR channels with multiple filters: short wavelength channel operated between 2.5 and 5.2 //m 
and the long wavelength operated between 4 and 17 . The other instrument relevant to GIB 
studies was ISOPHOT which covered the wavelength range from 2.5 to 200 /xm but had much 
smaller arrays. The ISO telescope had no shutter making absolute calibration difficult for studies 
of the isotropic component of the GIB. 

There were several important findings concerning the mid- to far-IR GIB from the ISO surveys. 

IR number counts indicate that the IR-luminous galaxies evolved more rapidly than their 
optical counterparts and make a substantial contribution to the star formation at higher z (Elbaz et 
al 1999) which in turn makes them the dominant contributors to the GIB at mid-IR bands. Elbaz et 
al (2002) computed the contribution from distant galaxies observed by deep ISOGAM cxtragalactic 
surveys at the central wavelengths of 6.75 and 15 /im to the GIB. The 6.75 //m sample was strongly 
contaminated by Galactic stars, whereas at 15 /Ltm the stars are readily separated from galaxies 
using optical-MIR color-color plots. They estimate that in this way they accounted for the GIB 
from dust in the starburst galaxies out to z ~ 1.5 and ISO measurements at 15 and 170 //m reveal 
that the bulk of GIB at these wavelengths comes from galaxies at 1.2 (Dole et al 2001, Elbaz et 
al 2002). They find that galaxies above their completeness level of 50 ^Jy produce at least 2.4ib0.5 
nWm~^sr~^ contribution to the total GIB at 15 //m . (For comparison, the IRAS galaxies at 60 
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/xm with completeness limit of 0.5 mJy contribute only ~ 0.15 nWm ^sr ^ to the CIB at 60 fim .) 

Lagache Sz Puget (2000) have analyzed a sub-field from the deep extragalactic FIRBACK 
survey with the ISOPHOT instrument at 170 jum . The survey covered 4 deg^ and they selected 

for analysis a subficld of low cirrus emission covering 0.25 deg^. 24 extragalactic sources were 
detected and subtracted up to the confusion limit of 3a = 67 mJy. The pixel size of the maps 
was 9p = 89.4" implying the Nyquist frequency of ^Nyquist = (2^p)~^ = 0.3 arcmin"^. On scales 
corresponding to spatial frequencies below the Nyquist limit k < 0.3 arcmin"^ the results lead to 
upper limits on the CIB fluctuations at 170 /xm that correspond to about 0.5 nWm~^sr^^ at 0.5 
arcmin. 



5.5. 2MASS 

The 2-Micron All Sky Survey (2MASS) uses two 1.3-m Cassegrain telescopes, one at Mt. 
Hopkins in the Northern Hemisphere, and one at CTIO in the Southern Hemisphere. Each telescope 
is equipped with a three-channel camera, capable of observing the sky simultaneously at J (1.25 
lim ), H (1.65 iim ), and (2.17 /xm ) at a scale of 2'.' per pixel. As the telescopes are scanned 
in declination, individual 1.3-second sky frames are imaged on an overlapping grid by stepping 1/6 
the array. The frames are combined six at a time, to form the standard 2MASS Atlas images of 
size 512 X 1024 pixels with re-sampled l^pixels, and an effective integration time of 7.8 seconds 
per pixel. Hereafter, we will use the term 'image' to refer to the calibrated 2MASS Atlas images 
which have been co-added to an effective integration time of 7.8 seconds. The 2MASS photometric 
stability is <0.02 mag in all bands. A detailed description of the calibration process can be found 
in the 2MASS Explanatory Supplement (Cutri et al 2003). 

The standard 2MASS images are too shallow to be useful for CIB studies. But a limited 
number of standard stars had to be observed repeatedly each night, and for several months at a 
time, to establish the photometric zero-points for the data. The 2MASS standard stars were drawn 
from near infrared standard star catalogs (e.g. Persson et al. 1998, Casali and Hawarden 1992). 
Each calibration observation consisted of six independent scans of a calibration field. Each scan is a 
mosaic of 48 images, and the scans were made in alternating north-south directions, each displaced 
5 ''in RA from the previous one to minimize systematic pixel effects. Repeated observations of 
calibration fields during a night at a variety of elevation angles were also used to develop long-term 
atmospheric extinction statistics. By collecting the calibration scan data, which spanned nearly 6 
months of repeated observation, effective integration times exceeding one hour could be achieved 
for a small number of sky locations. 

Kashlinsky et al (2002) and Odenwald et al (2003) used the long exposure 2MASS data for 
CIB analysis and reported the first detection of small angular scale fluctuations in the near-IR CIB. 
They have selected data from the 2MASS Second Incremental Data Release between 1998 March 19 
and 1999 February 20, and used the 2MASS, on-line catalog of completed calibration fields in order 
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to identify fields that had the largest number of repeated measurements. One of these, 90565N, 
with a total exposure 1 1 hour was selected for further study. The field has Galactic and Eclip- 
tic coordinates (243°, 27°) and (21°, 35°) respectively, around the star P565-C in the constellation 
Hercules, 5° north of the bright star A Ophiuchus. The data analyzed by them consisted of 2080 
calibrated 8.6'xl5' frames covering an 8.6' x 1° swath oriented north-south, obtained during obser- 
vations at CTIO between Apr and Aug 1998. The images from individual exposures were co-added 
to produce a 8.6' x 1° field from which resolved sources (both stars and distant galaxies) and other 
artifacts were removed by various methods to ensure consistency and robustness of the final results. 
After assembling the final image with a total exposure >3,700 sees the data were calibrated and 
the field was divided into seven square patches 512" on the side. In each patch, individual sources 
(stars and galaxies) were removed by an iterative procedure where each pixel with surface bright- 
ness exceeding 3 standard deviations for the patch was excised along with 8 neighboring pixels. 
The 2MASS beam of ~ 2" subtends comoving scale ^30/i~^Kpc at z=l corresponding to a large 
galaxy so the sources near the the clipping threshold are unresolved. Because of the variation in 
background level and associated noise, the patches were clipped to different point source absolute 
flux levels. In band these varied from 18.5 to 19 Vega magnitudes. The clipping left >90% of 
the pixels in the patches, which provided a good basis for a Fourier analysis of the diffuse light. For 
the final images, the fiuctuation spectrum of the residual diffuse emission was computed; its slope 
in the co-added images is consistent with that of galaxy clustering. They estimated the contri- 
butions to the power spectrum from atmospheric fluctuations, remaining Galactic stars and cirrus 
emission, zodiacal light, instrument noise and extinction. These components had different slopes 
and negligible amplitudes compared to the detected signal, suggesting that the diffuse component 
fluctuations in the final image are dominated by the GIB. Because in this type of analysis fairly 
faint galaxies (K^IQ) have been removed, the resultant GIB provides a probe of still fainter sources 
which are presumably located at high redshifts, and allows to start isolating the contributions to 
the GIB from various cosmic epochs. 

Fig. 8 shows one of the patches in band after 7.8 sec integration (left panel - standard 
2MASS product) and after the integration used in the analysis (middle panel). The right panel 
shows the flnal 3-color image (J, H and K^) which was identifled with the map of the near-IR 
GIB from galaxies fainter than Kg — +19. An important step forward in this measurement is the 
removal of progressively fainter and more distant galaxies. This allows to begin to identify how 
much of the GIB comes from earlier times. Observations and theory both suggest that + 19 
galaxies are typically located aX z^l (Gowie et al 1996, Kashlinsky et al 2002). 



5.6. Results 

The results on the mean GIB levels from the measurements described above are shown in 
Fig. 9. Only detections arc shown; upper limits, in the absence of detections, are discussed 
below for each range of wavelengths. Hence the gap in the flgure at the MIR wavelengths where 
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Fig. 8. Patch 3 of 512" X 512" in size from the 90565 field used in the deep 2MASS data analysis of Kashlinsky, Odenwald 

et al (2002) and Odenwald, Kashlinsky et al (2003). (Left) Individual 2MASS images for 7.8-second integration in Ks band. 
(Middle) Corresponding deep-integration co-added Ks band image for an effective integration time of 3,900 seconds. (Right) 
3 color image of the patch after co-adding, clipping and de-striping in Fourier space. Red, green and blue correspond to J, H, 
Ks bands. The structure in the imago is identified with near-IR CIB from galaxies fainter than Ks ~ 19. 

zodiacal foreground is brightest. Total fluxes from the observed galaxy population are also shown 
for comparison; they are discussed at length in Sec. 6. 

Displaying measurements of CIB fluctuations is more complicated because of the extra spatial 
dependence and some of their measurements will be discussed for each wavelength range separately. 
Figure 10 shows the amplitude of CIB anisotropics at ~ 0.5° from the DIRBE data analysis of 
Kashlinsky & Odenwald (2000) and Fig. 11 shows the results for CIB fluctuations from the deep 
2MASS data analysis described in Sec. 5.5. 

Below we itemize the status of CIB measurements in three IR wavelength ranges: NIK, MIR 
and FIR. 

5.6.1. Near-IR 

In the near-IR there are detections of the mean isotropic part of the CIB based on analysis 
of DIRBE and IRTS data (Dwek & Arendt 2000; Matsumoto et al 2000, 2003, Gorjian et al 2000, 
Wright & Reese 2000, Cambresy et al 2001). The measurements agree with each other, although 
the methods of analysis and foreground removal differ substantially. They also agree with the 
measured amplitude of CIB fluctuations. The results seem to indicate fluxes significantly in excess 
of those from observed galaxy populations. Total fluxes from galaxies are shown with filled circles 
in Fig. 9; in near-IR they saturate and are a factor is ~ 2 — 3 below the detected CIB. 

Detections of CIB fluctuations come from three independent experiments and are consis- 
tent with each other. At 0.5° Kashlinsky & Odenwald (2000) find a statistically significant CIB 
anisotropy in the DIRBE first 4 channels. Matsumoto et al (2000, 2003) measured a spectrum of 
CIB anisotropics on degree scales from IRTS data band-averaged to ~ 2 /im ; their results are 
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Fig. 9. — NIR: Mean levels of the near— IR CIB from the IRTS (crosses) from (Matsumoto et al 2000) and open symbols from 
COBE/DIRBE data. Open triangle shows the CIB results of the DIRBE J-band analysis (1.65 /^m ) from (Cambrcsy et al 2001) 
and filled diamonds at 2.2 /im from (Gorjian et al 2001), at 3.5 fj,m from (Dwek & Arendt 1998). FIR,: Thick solid line 
with shade showing the uncertainty shows the results from CORE FIRAS data analysis from Puget et al (1996) and Fixsen 
et al (1998). DIRBE data detection at 140 and 240 nm is shown with crosses from Hauser et al (1998) and open squares 
from Schlegel et al (1998). Filled pentagon shows the lower limit on the CIB from the DIRBE data at 100 fj,m from Dwek et 
al (1998). Open diamonds correspond to results of DIRBE analysis by Finkbcincr et al (2000). Different results at the same 
wavelength are slightly shifted for clearer display. Cumulative fluxes from observed deep galaxy counts are shown with 
the filled circles from Table 5. The cases where the flux from progressively fainter galaxies does not saturate or the saturation 
is not as clear are marked with an upward arrow. 

consistent with Kashlinsky & Odenwald (2000). Because the beam was large in both instruments, 
no galaxy removal was possible in the data and the CIB anisotropics arise from all galaxies, i.e. 
from z = to the earliest times. Based on the amplitude of the present-day galaxy clustering, 
the NIR detections exceed the expectations from galaxies evolving with no or little evolution by a 
factor of ~ 2 — 3 and are consistent with the measured (high) mean CIB levels in the NIR. 

Figure 11 summarizes the 2MASS-based results and compares them with the other measure- 
ments. Analysis of deep 2MASS data (Kashlinsky et al 2002, Odenwald et al 2003) has allowed 
to remove galaxies out to X ~ (18.5 — 19) and measure the spectrum of CIB anisotropics in J, 
H, and K 2MASS bands on sub-arcminute angular scales. These galaxies are typically located 
at cosmological times when the Universe was less than ~ half its present age and when extrapo- 
lated to the present-day the 2MASS based results are consistent with the DIRBE- and IRTS-based 
measurements. This indicates that the near-IR CIB excess, if real, must originate at still earlier 
times. 

All NIR CIB results are mutually consistent, exceed the simple extrapolations from the ob- 
served galaxy populations with little evolution and, if true, probably imply significant energy-release 
activity at early cosmic times. 
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Fig. 10. — Detections of and upper limits on the RMS CIB fluctuation at 0.5° from the DIRBE data analysis by Kashlinsky 
& Odenwald (2000). 

5.6.2. Mid-IR 

Because the zodiacal light is so bright at these wavelengths there is no direct detection of the 
CIB between ~ 10 and ~ 50 fim . There are low upper limits though. Some come from the 
fluctuations analysis of the DIRBE data by Kashlinsky, Mather & Odenwald and use the fact that 
zodiacal light, while very bright, is also very smooth. Whereas light emitted by galaxies must be 
clustered with amplitude of about ~ 10% in flux fluctuations on the DIRBE beam scale. This leads 
to upper limits of about ~ 10 nWm~^sr~^ at these wavelengths. 

Other upper limits come from scattering of 7-ray and CIB photons producing electron positron 
pairs if E^E^^^^ > 2(meC^)^. For 7-rays the cross-section is peaked in the IR wavelength range, 
Apeak — 2.5pj^ /Ltm (Stecker et al 1992). If a drop in the spectrum of an extragalactic 7-ray source 
in TeV range can be observed this should lead to limits on the CIB in the corresponding wavelength 
range (Stecker & dc Jager 1993). Dwek & Slavin (1994) applied this analysis to the 7-ray spectrum 
of Mrk 421. The best upper limits come from analysis of the 1997 7-ray outburst of the blazar Mrk 
501 and are ~ 5 nWm~^sr~^ from 5 to 15 ^m (Stanev & Franceschini 1998, Renault et al 2001). 

The ISO and Spitzcr deep counts provide important lower limits on the CIB at these wave- 
lengths which are very close to the upper limits from 7-ray and fluctuations analyses. It would be 
thus very surprising if the MIR CIB is not inside this 2.5 nWm^^sr"^ ^vl^^^h nWm~^sr~^ range. 



5.6.3. Far-IR 

Far-IR CIB measurements come from the COBE FIRAS (Puget et al 1996, Fixsen et al 1998) 
and COBE DIRBE (Schlegel et al 1998, Hauser et al 1998, Finkbeiner et al 2000) data analysis. 
They are broadly consistent although the COBE DIRBE detections at 140 and 240 /xm give 
larger fluxes (at almost 2-a level) and the Finkbeiner et al (2000) analysis gives an even bigger 
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Fig. 11. The angular spectrum of the CIB fluctuations from the 2MASS analysis by Kashlinsky, Odenwald ct al (2002) 

and Odenwald, Kashlinsky et al (2003) at J, H and K bands. Triangles correspond to the results from the patch from which 
galaxies out to Kg < 18.6 have been removed and crosses to another 512" X 512" patches with galaxies removed out to Kg ^ 19. 
Open diamonds with 92% errors show the fluctuation on larger angular scale from Kashlinsky & Odenwald (2000) analysis of 
COBE DIRBE data. Thick soUd line shows the CIB result from IRTS analysis (Matsumoto et al 2000, 2003). Both the DIRBE 
and IRTS analyses included all galaxies and are consistent with each other and the 2MASS data when account is made of the 
contribution from the removed (brighter) galaxies (Odenwald, Kashlinsky et al 2003). 

discrepancy at 60 /xm if the COBE FIRAS results are simply extrapolated to that wavelength. 
Hauser et al (1998) discuss these differences and point that the differences are probably due to 
different instrumental calibrations. 



where A = (1.3 ± 0.4) x 10-^/3 = 0.64 ± 0.12, T = 18.5 ± 1.2 K, z/q = lOOcm"! = 3000 GHz and 
Bi, is the Planck function. In eq. 24 the uncertainties are highly correlated with correlations +0.98 
for A and P, -0.99 for A,T and -0.95 for /3,T (Fixsen et al 1998). 

The bulk (likely all) of the CIB at these wavelengths comes from dust in high-z and present-day 
galaxies and the SCUBA galaxy counts at 850 fim (Small et al 1997) compare favorably with eq. 



Fixsen et al (1998) suggest a simple fit to the CIB at A> 100 fim : 




(24) 



24. 
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5.6.4- Bolometric CIB flux 

Table 4 shows the integrated IR fluxes over the three ranges of wavelengths where direct CIB 
measurements were possible to meaningfully integrate (hence, e.g. the gap between 4 and 10 /xm ). 
In the NIR out of 37 nWm~^sr~^ approximately 13 nWm^^sr"^ is contributed by the Cambresy et 
al (2001) measurement at the DIRBE's J band; the remaining come from the IRTS data and other 
measurements out to ~ 4 ^um . In order to evaluate the contribution in the MIR region we assumed 
an upper limit between ~ 10 and 50 /v,m of vl^lz 5nWm~^sr~^ from the 7— ray absorption. (An 
upper limit from fluctuations analysis would be a little higher). In the far-IR we integrated the 
approximation eq. 24; including the other DIRBE detections will increase this number. 

The table shows that CIB carries very substantial fluxes testifying to significant energy release 
throughout the evolution of the Universe. For comparison the total flux from the CMB is only 
~ one order of magnitude higher. Table 5 and discussion in Section 6 show that normal stellar 
populations observed in deep galaxy counts surveys contribute about as much CIB flux in the near- 
IR as the total far-IR CIB in Table 4, meaning that approximately half the energy produced by 
stars has been absorbed and re-emitted in the far-IR by dust in (high-z) galaxies. The observed 
NIR galaxies in J,H,K bands contribute only ~ SnWm^^sr"^ . The total fluxes from the deepest 
counts in Spitzer data are also shown although at 5.8 and 8 /xm does not saturate as clearly as 
at 3.6 and 4.8 jim. . If one integrates the flux from galaxies in J, H, K bands and the Spitzer 
3.6 /im one gets 8.6 nWm~^sr~-^ . Extending the range of integration to the other Spitzer IRAC 
channel at 4.5 fim. where the total flux from galaxies saturates (sec Sec. 6) gives 9.7 nWm^^sr^^ 
and integrating over galaxies from 2MASS and all Spitzer IRAC channels gives total flux of ~ 11.4 
nWm~^sr~^ between 1 and 10 jiva. . The observed galaxy populations contribute only a small 
fraction of the observed flux. 

This means a substantial excess of the total CIB measurements over that observed to come 
from total galaxies. It then appears likely that either 1) the total contributions from ordinary 
galaxies are significantly underestimated, or 2) the bulk of the CIB comes from still earlier epochs 
than probed by the systems accessible to current telescopic studies, e.g. from the Population III 
era, or 3) all the the NIR measurements contain systematic errors mistaken for CIB. The first of 
these possibilities is discussed in Sec. 6 and shown not very likely, the second is discussed in Sec. 
7 and the third possibility should always be borne in mind when interpreting the very difficult 
measurements of the CIB. 

Historically the excess CIB flux in the near-IR has been nearly simultaneously and indepen- 
dently detected for the isotropic part of the CIB by Dwek & Arendt (1998) and for fluctuations by 



Table 4. Bolometric/integrated CIB flux (nWm ^sr ^ ) 



NIR: 1 /im < A < 4 /xm 


MIR: 10 /im < A < 50 /im 


FIR: 120 /im <A< 1,000 /im 


Total; 1 /im J; A ^ 1 mm 


37 


< 8 


13 


50 <FtotalS 58 
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Kashlinsky & Odenwald (2000). ^ 



6. 'Ordinary' contributors to CIB 

In this review we use the term "ordinary" galaxy populations to refer to galaxies made of the 
observed normal metal- enriched stellar populations containing Population I and II stars and the 
dust produced by them. Observations from the Hubble Deep (and Ultra Deep) Field suggest that 
metal enrichment has happened fairly early on in the history of the Universe. Similarly, data on 
quasar emission and intrinsic absorption lines indicate that solar metallicities in the surrounding 
medium have already been reached at 2; ^4 (Hammann &; Ferland 1999). This is also suggested 
by the existence of the oldest Population II stars with non-zero metal abundances. Star formation 
processes in the metal enriched medium (metallicity Z^IO^'^ — IO^^Zq) should lead to formation 
of populations with stellar properties not radically different from those observed today. These 
populations must have been preceded by completely metal-free (first) stars, the so-called Population 
III that formed at very high redshifts, 10. Population III stars likely form a completely different 
class of stellar objects, probably have a very different range of masses and also form in a different 
environment. Their case is special, necessarily more speculative, and their possible contribution to 
the CIB is discussed in the following section. 

Before we proceed a simple "back-of-the-envelope" evaluation is in order for the amount of the 
near-IR CIB produced by 'ordinary' stellar populations. If such stars had Salpeter-type IMF most 
of the light would be produced by the high mass end of the stellar mass spectrum, i.e. by stars 
that produced the metal abundances {Z ~ Zq) observed today. If galaxies whose stellar contents 
contribute ~ 3 x 10~^h~^ to the critical density (pcr), then during their main sequence the stars 
would have produced energy density of eQ^pcrC^, where e is the radiative efficiency. For hydrogen 
burning e = 0.007 and this cfRciency is reached for massive stars, MJJIOMq. For present-day 
stellar populations with a Salpeter IMF and a lower cutoff of O.lAf0, the effective efficiency is 
only e ~ 0.001 (Franceschini et al 2001) ^. A fraction /(A > j^fim) of the emitted energy will 
contribute to the near-IR CIB at wavelengths longer than 1 fim. , assuming it was not absorbed or 
re-processed. This would lead to the CIB flux of: 

Ul. ~ ^eZn^PcrC { ) = 8/. ^ ^ 3 X 10-3/.-1 ^ 1 + Z 



*The paper by Kashlinsky & Odenwald (2000) has been received by the Ap.J on April 30, 1998, but due to an 
an editorial mishandling the manuscript was not sent to an independent referee for review until over a year later, in 
May 1999. 

^Solar mass stars burn a core of about ~ 10 — 15 % of their hyrdogen mass during main-sequence stage. During 
their red giant stage the hydrogen burning core is ~ 0.5Mq for M< 2.5Mq (Sweigart, Greggio & Renzini 1990). For 
the IMFs where the bulk of stellar mass is locked in the low end of the stellar mass-spectrum this results in effective 
e considerably smaller than the canonical value of 0.007 (Sweigart, private communication) 
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The stars responsible for producing metals are predominantly hot stars with surface temperature 
T*ii (1 — 2) X 10^ K. In this case /(A > j^^fim) the CIB flux comes from the Rayleigh-Jeans part 
of the spectrum, where / oc T~^(l + z)^ and the rcdshift-dcpcndcnt term on the right-hand-side of 
eq. 25 increases toward high z because a larger part of the emitted spectrum is shifted into the IR 
today. With T* ~ 10^ K and z ~ 2 — 3 one gets contribution to the CIB around a few nWm~^sr~^ 
, but if there were very hot (massive) stars at very early times, their contribution may well exceed 
that of the 'normal' stellar populations. 

Silk (1996) makes a similar argument to eq. 25 differentiating between the metals (iron) 
produced in early and late type galaxies and their respective mass-to-light ratios, and Hauser & 
Dwek (2001) normalize instead to helium production; both get similar numbers. Helium is produced 
by less massive (and longer living) stars, and they could produce an amount of CIB comparable to 
(or less) than eq. 25. 

Comparison with the ncar-IR levels of CIB implies that, if the measurements are true, it is 
difficult to explain these high levels by left-over emissions from galaxy populations with IMF similar 
to that of today. The levels of the near-IR CIB can be increased in one (or all) of the following 
ways over eq. 25: 1) If there exists a population of stars that can emit significant amounts of light 
without producing metals (as happens e.g. in stars more massive than ~ 240Mq, Heger 2003), one 
can have extensive period of energy production that would not violate the observed metallicities. 
2) If emission comes from populations that contribute a larger fraction of rJbaryon — 0.044 than the 
stellar populations of today which have fi* <C i^baryon- Interestingly, most of the baryons today are 
in dark form (Fukugita, Hogan &; Peebles 1998), the dark baryons possibly being in remnants of 
the Population III era. 3) And if there existed at high z a population of very hot stars, the bulk of 
their emission will today be shifted into the NIR bands. We discuss such possibility in the context 
of Population III stars contribution in the following section. 

The remainder of this section discusses the various aspects necessary for estimating the con- 
tributions from ordinary (metal-rich) stellar populations in galaxies to the CIB from NIR to FIR. 

6.1. IMF and star formation history 

SFR indicators include UV and emission-line luminosities indicative of light from short-lived 
massive stars. Conversion to SFR depends on the assumed IMF, but once the IMF assumption is 
made, the relative SFR is less sensitive to IMF. The conversion of the UV luminosities into the 
SFR must also correct for the absorption of UV photons by dust which is generally associated with 
the same young stars. 

In the first calculation of the stellar IMF, Salpeter (1953) approximated the Solar neighborhood 
measurements with IMF having a power law of slope —1.35 and it still is commonly used in modeling 
synthetic stellar populations. Scalo (1986), from measurements of the Solar neighborhood, gives an 
IMF with a mass fraction peak at ~ 0.5 — IMq. While very useful and necessary in constructing 
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synthetic stellar population models, it is not clear whether the assumption of the universal IMF 
is valid in the real Universe. E.g. it cannot reproduce Ha luminosities (Kennicutt et al 1994) or 
match mean galaxy colors (Madau et al 1998) because it has too few stars of M> WMq. 

In galaxies with on-going star formation, the composite UV spectrum is approximately flat. If 
this is combined with the Salpeter IMF between 0.1 and 100 Mq in stellar mass, one can relate the 
rest frame UV luminosity to SFR (Madau et al 1998) for galaxies with continuous star formation 
on timescales longer than >0.1 Gyr: 



Because in starburst galaxies a significant fraction of luminosity, which is dominated by young stars, 
is absorbed and re-emitted in the infrared by dust, the IR luminosity can also be used as tracer 
of SFR. Assuming that the bolometric luminosity is dominated by young stars, Kennicutt (1998) 
obtained in the optically thick limit (when IR luminosity measures the total bolometric luminosity) 
for continuous star bursts lasting at least 10 Myr, solar metal abundance and a Salpeter IMF: 



so luminous IR galaxies with Ljr, > IQ^^Lq form more than 20 Mq/jt in stars. Note that these 
calibrations depend on the IMF in star forming galaxies, since the integrated UV spectrum is 
dominated by massive stars and is assumed to lead to dust (re)emission in the mid- to far-IR. 

Lilly et al (1996) used galaxies from the Canada-France Redshift Survey (CFRS) to construct 
the comoving luminosity, Ci,{z), of the Universe out to z ~ 1 at rest-frame wavelengths of 0.28, 0.44 
and 1 //m . The sample included 730 I-band selected galaxies, with luminosities extrapolated to the 
above rest wavelengths using interpolations from the BVIK photometry with SEDs interpolated 
from synthetic spectra designed to match the observed colors. The present-day I-band corresponds 
to the blue-band for rest frame emission at z ~ 1. They found observational evidence for an increase 
of Ci,{z) out to z ~ 1 in the galaxy rest-frame UV to near-IR (1 ^m ) bands and argued that the 
data requires SFR oc (At)~^-^ where is the time elapsed since the initial star formation. 

With the Lyman-dropout technique (Steidel et al 1996), galaxies are now found to progressively 
higher rcdshifts (Giavalisco 2002). Flux amplification by gravitational Icnsing by known clusters of 
galaxies can lead to galaxy detections at even higher redshifts as amplifications of order ~ 10 — 100 
can be achieved for distant galaxies lying on caustics of the cluster potential. Kneib et al (2004) 
measured optical and IR photometry of one such candidate galaxy at 2; ~ 7 and find that, while 
only < 1 Kpc in extent, it undergoes star formation at the rate of ~ 2.6M0/yr. A claim was made 
recently of a galaxy at 2; ~ 10 among the objects lensed by a known cluster A1835 (Pello et al 
2004, but see Bremer et al 2004). Bouwens et al (2004) have applied the zgso dropout search to 
the Hubble Ultra Deep Field (UDF) data and detected four (and possibly five) potential 2; ~ 7 — 8 
objects that are seen in J and H bands, but have no detection shortward of 8,500 A(one of the 
candidates has a weak signal at 8,500 A). The objects are strongly clustered covering a total area of 
1 arcmin^ and also appear very blue in their rest frame. Using VLA observations of the molecular 



SFR(MQ/2/r) = 1.4 x 10-2^L^(erg/sec/Hz) 



(26) 



SFR(M0/2/r) = 1.71 x 10-^"Lir(8 - lOOO/im)(L0,boi) 



(27) 
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gas in the underlying galaxy of the current record-holder z = 6.4 quasar, Walter et al (2004) detect 
~ 5 X 10^ Mq in molecular gas extended to ~ 2.5 Kpc. The amount of gas is comparable to what is 
found in the nearby ULIRGs and its brightness temperature is about that of the nearby starburst 
centers. Clearly stars must have formed early on in the Universe's history and at vigorous rates. 

Madau et al (1996) have analyzed the Hubble Deep Field (HDF) observations reaching AB 
magnitudes ^28 in four bands from 3000 A to 8140 A. They modeled the HI opacity and used 
stellar population synthesis models to separate low- and high-z galaxies using the Lyman continuum 
drop-out method (Steidel et al 1996). Prom the images they constructed a sample of star-forming 
galaxies from which the comoving UV luminosity density was estimated out to z ~ 4. Assuming 
that the UV luminosities give a complete view of star formation, the SFR density can be evaluated 
from the conversion factors for the Salpeter IMF (Madau ct al 1998). The derived SFR, known 
as the Madau curve, was found to rise rapidly toward early times possibly peaking at z ~ 2. 
Estimating the comoving UV luminosity density from the high redshift galaxies found in the UDF 
data (Bouwens et al 2004) suggests a factor of >3.5 drop in SFR from z = 3.8 to z ^ 7.5 consistent 
with what is found at z ~ 6 (Dickinson et al 2004, Stanway et al 2004). 

Guiderdoni et al (1997) argued that most galaxies are hidden by dust which leads to substantial 
underestimating of the high redshift SFR deduced from optical surveys. They constructed models 

of several galaxy populations to separately match the optical surveys' data and the far-IR CIB 
deduced from the COBE/FIRAS data. The galaxy population responsible for the FIR CIB has 
a high SFR at early times, information on which will be missing from SFR deduced from optical 
bands. Indeed, if most of energy is (re)emitted at IR wavelengths, the SFR may be underestimated 
when deduced from surveys in the UV bands. 

Rowan-Robinson et al (1997) analyzed ISOCAM observations of the HDF at 6.7 and 15 /xm 
(Serjeant et al 1997), where thirteen HDF galaxies have been detected by ISO in the 12.5 hour 
exposure. In eleven of these, there was a substantial mid-IR emission excess, which they interpreted 

as dust emission from a strong star-burst. The SEDs were modeled from 0.3 to 15 /im using the 
available photometry from visible to NIR bands and the Bruzual-Charlot stellar evolution models 
with a Ray leigh- Jeans fall-off beyond 2.5 /im , and the conversion to SFR was done using starburst 
modeling with emission from dust and PAHs. If these star-forming galaxies are typical of the fainter 
HDP galaxies, then Rowan- Robinson et al argued that the true SFR remains flat at z^ 1.5 instead 
of falling off. 

6.2. Normal stellar populations 

Normal stellar populations are defined as those with Population I and II metallicities and the 
(locally) observed IMP, usually taken to be either of Salpeter or Scalo type. 

Constructing realistic SEDs of galaxies is important for proper comparison of contributions of 
galaxies of various types, from various redshifts and at various wavelengths to the observed CIB. 
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In the near-IR such SEDs can be approximated with the spectrum and its evolution derived from 

synthetic stellar population models. This approach requires input of a universal IMF, assumptions 
about the rate of star formation, self-consistent treatment of stellar evolution, metal enrichment 
etc. The early studies have been pioneered by Tinsley (1976) with Bruzual (1983) developing the 
state-of-the art models that can be compared with detailed observations. The current versions of 
the synthetic stellar models e.g. by Bruzual-Charlot (1993) or PEGASE (Fioc &; Rocca-Volmerange 
1997) codes are available on the world- wide- web. The modeling has to be done carefully and many 
priors are required, since the SEDs in turn depend on many parameters: IMF, metallicities, galaxy 
ages, dependence of each of them on other, galaxy luminosity, SFR, etc. Still, such modeling gives 
a fair representation of reality and can certainly provide approximate answers. 

For the purpose of calculations in this section we used the PEGASE ^ code to construct SEDs 
of Early and Late type galaxy stellar populations as follows: Early type galaxy stellar populations 
were assumed to be all described by Salpeter IMF and form at some early time in a single short 
burst of star formation lasting 0.5 Gyr. Late type stellar populations in galaxies were assumed to 
form via an on-going star formation with Scalo IMF starting at some early epoch. The PEGASE 
models were run for a grid of galaxy ages and metallicities. The left panel in Fig. 12 shows 
the SED for the two types of populations for ages between 10 and 14 Gyr and metallicities Z = 
0, 10~^,2 X 10~^,5 X 10~^,10~^. The symbols show the value of the relative flux when averaged 
over SDSS and 2MASS filters to compare with data discussed below. (Plus signs correspond to 
early types and diamonds to late types.) The NIR luminosities are dominated by red giants in the 
old stellar populations and are directly related to the total stellar mass. 

The right panel in Fig. 12 shows how the mass-to-light ratio varies with time for such pop- 
ulations. Basically the differences in the spectra reflect the ratio of young (< 1 Gyr) to old stars 
(Kennicutt 1998). 

Shaded regions show the approximate range of the mass-to-light ratios of the two populations 
(e.g. Fukugita, Hogan & Peebles 1998) plotted for tcosm < 14Gyrs, the age implied by - inter alia 
- WMAP measurements. They agree reasonably well with the plots especially if early populations 
are older than a few Gyrs, although there may not be one 'magical' value to plot for each of the 
galaxies (Faber & Gallagher 1979, Roberts & Haynes 1994). 

In order to compare with GIB measurements, the synthetic SEDs must be averaged over the 
entire ensemble of galaxy populations. This includes averaging over all morphological types and, 
within each type, averaging over the luminosity function. The first is important because early- 
type galaxies contribute substantially to the near-IR EBL (Jimenez & Kashlinsky 1999). The 
second step is also important because, at least for early-type galaxies, the mass-to-light depends 
on galaxy mass or luminosity as suggested by the fundamental plane measurements for elliptical 
galaxies (Faber et al 1989, Djorgovski &: Davis 1987). In the visible, the fundamental plane implies 



'http://www2.iap.fr/uscrs/fioc/PEGASE.html 
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Fig. 12. Runs of synthetic galaxy SED spectra using PEGASE. Green lines correspond to late types, red lines to 

early types stellar populations. Left: black symbols correspond to the SEDs integrated over the SDSS and 2MASS fil- 
ters. The sets of lines are shown for five metallicites oi Z = 0,0.001,0.002,0.005,0.01 and five ages between 10 and 14 
Gyr. Right: Mass-to-{bolometric) light ratios in solar units for early and late stellar populations are plotted vs their age for 
Z = 0,0.001,0.002,0.005,0.01. Cross-shaded areas correspond to typical mass-to-light ratios observed locally at 2 = in the 
early (top cross-shaded area) and late populations with uncertainties from Fukugita, Hogan & Peebles (1999); the observed 
variations in the mass-to-light ratios are larger. 



M/Lb oc L'^ with Kb — 0.25. The slope varies systematically with wavelength from 0.35 i-im 
through the K band (Pahre et al 1998). Also, measurements of cluster galaxies out to z ~ 1, or 
~ ^ look-back time, indicate that the fundamental plane preserves its logarithmic slope while the 
amplitude changes (Treu et al 2002, van Dokkum et al 2001, 2004, van der Wei et al 2004). The 
origin of this wavelength trend may lie in various quantities, like systematic variations in the IMF 
along the luminosity sequence (Pahre et al 1998), or with metallicity (Jimenez &; Kashlinsky 1999). 



6.3. Dust emission from galaxies: mid-IR to sub-mm 

The near-IR CIB data provide information on the brightness and structure of largely un- 
extincted starlight. The mid- and far-IR observations are necessary to complete the picture, by 
providing the corresponding information for starlight that has been absorbed by dust and thermally 
re-radiated. In many environments, even at early epochs, this re-radiated emission is a substantial 
or even dominant fraction of the total luminosity associated with star formation. 

There are excellent reviews on the subject of IRAS galaxies (Soifer, Neugebauer & Houck 
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1987b), ultra luminous infrared galaxies (ULIRGs, Sanders & Mirabel 1996) and sub- mm SCUBA 
galaxies (Blain et al 2002) and the reader is referred to them for more details. This subsection 
summarizes the established properties of these sources as they relate to the CIB. 

The range of galaxy luminosities in the mid- to far-IR greatly exceeds that in the UV, optical 
and NIR. E.g the mid- to far-IR luminosities of IRAS galaxies vary from below a few times 10^^ 
erg/sec to over 5 x 10^^ erg/sec. Many of the more luminous IRAS galaxies are starbursts. The 
bulk of the IRAS detected galaxies are late-type spirals, with only a handful belonging to the early 
(E or SO) type. The shape and slope luminosity function of IRAS galaxies differs from that of the 
visible band galaxies and is best approximated as two different power-laws at the high and low ends 
of the galaxy flux distribution (Soifer & Neugebauer 1991). At the bright end of the bolometric 
luminosity function, the bright mid- to FIR galaxies become the dominant galaxy population in 
the present-day Universe. The ratio of the FIR to visible galaxy luminosity correlates with 60 to 
100 //m color temperature for most galaxies (Soifer et al 1987a,b). Soifer et al (1987), using the 
Draine &; Lee (1984) dust modeling, give a useful approximation for the amount of dust needed to 
produce a given FIR luminosity: 



This shows that even a tiny amount of dust can make the galaxy very bright in the mid- to 
far-IR. Zubko, Dwek k. Arendt (2004) have developed sophisticated ISM dust models by fitting 
simultaneously the far-UV to NIR extinctions, the diffuse IR emission and constraining the dust 
properties from the observed elemental abundances. Assuming that the dust mass-absorption 
coefficient k{\) oc A~"dust^ specific fiux at wavelength A received from a galaxy containing Mfjust 
a distance d away is given by (Dwek 2004): 

/.(A) = 8xl0^ hH(A/l.m)r,.J-lJ Jy 

(29) 

where Ao is some reference wavelength and ndust ~ 1 to 2. 

Most IR emitting galaxies can be classified as normal spirals, starbursts or AGNs. Starburst 
galaxies are especially bright at mid- to far-IR wavelengths and approximately half the galaxies 
brighter than Lfir — 3x 1O^°L0 have star-burst optical spectra (Elston et al 1985). For the ULIRGs 
their measured star formation rate cannot be maintained for longer than (at most) a Gyr, as at 
that rate all their ISM will be consumed. The most luminous IRAS galaxies show features of both 
being a starburst and an AGN. In the less luminous IRAS galaxies their MIR/FIR emission seems 
to be unrelated to their current SFR. For most IRAS galaxies that are not AGNs, the observed 
color-color correlations can be accounted by a two-component dust model, where warmer dust is 
associated with regions of active star formation and the cold component comes from cirrus in the 
galactic disk (Helou 1986). 

One can reconstruct the SED of a typical galaxy in the MIR to FIR range of wavelengths. 
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Fig. 13 shows such synthetic template: it was constructed accounting for continuum dust emission 
assuming T^ust = 35 K and the free-free emission at long wavelengths. FIR emission lines were 
scaled proportional to the flux and adopted from: the Galactic IR lines observed by FIRAS (Wright 
et al 1991), the PAH broad line features at A ^ 12 //m were taken from the Dwek et al (1997) Galaxy 
model of the DIRBE data, and the lines from the ISO Long Wavelength Spectrometer observations 
of nearby bright IR galaxies by Fischer et al (1999). This SED is in good agreement with the IRAS 
galaxies average used by Dwek ct al (1998), the model of Guiderdoni et al (1998) and with the 
template of sub-mm emission from Blain et al (2002), where comparison is given with the data 
from galaxy observations. 

The right panel in Fig. 13 plots the K-correction due to this SED (defined as K^ot = 
-2.51g[I^(i+^)(10pc)//^(10pc)]) for 100 ^um (solid line), 250 /im (dots), 500 fim (dashes), 1000 ^um 
(dashed dotted line) and 2000 /xm (dash-dot-dot). One can see the strongly negative K-correction 
at A ^500 fim and high z: galaxy flux should decrease much slower with increasing redshift and it 
can actually increase with z. This further implies that at sub-mm wavelengths the contribution to 
the GIB should be more heavily weighted toward high redshift sources. 
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Fig. 13. — Left: Typical synthetic sub-mm SED of dusty galaxies including lines from PAHs and ions, atoms and molecules 
as observed by FIRAS and ISO. 

Right: K-correction vs z at various wavelengths: 100 fim (solid line), 250 fim (dots), 500 /im (dashes), 1000 /im (dashed 
dotted line) and 2000 /xm (dash-dot-dot) 

Optical follow up of the IRAS catalog shows that on average about 30 % of the bolometric 
luminosity is emitted in the far-IR via thermal radiation from dust heated by stars (Soifer & 
Neugebauer 1991) and this number rises to up to 95% for ULIRGs (Sanders Sz Mirabel 1996). 
Guiderdoni et al (1997) argue that most galaxies are hidden by dust. The mid- and far-IR number 
counts indicate that the IR-luminous galaxies evolved more rapidly than their optical counterparts 
and make a substantial contribution to the star formation at higher z (Elbaz et al 1999). Locally, 
typical galaxies radiate a little more in the UV to optical bands than in mid-IR to far-IR (Soifer 
Sz Neugebauer 1991). But, while luminous IR galaxies contribute a negligible amount to the local 
rate of star formation, they are major contributors at high z. 
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Spinoglio et al (1995) assembled a sample of ~ 900 IRAS galaxies out to the 12 jj,m limit 
of ~ 0.2 Jy and find that the 12 ^m luminosity is oc Ljr total ^ind suggest that their selection by 
the Li2^m is approximately equivalent to that by the total (bolometric) luminosity . They find 
that the 60 and 25 nm luminosities rise more steeply than linear with the bolometric -LiR,totai 
and the opposite is true for optical bands implying that more luminous disk galaxies have more 
dust shrouded stars. Malkan and Stecker (1998) used these empirical correlations together with 
assumptions of luminosity evolution to reconstruct the CIB across from near- to far-IR. 

Several mid-IR deep surveys have been conducted with ISOCAM, reaching sensitivity levels 
of ~ 50/xJy at 15 /xm (Elbaz et al 1999). Gravitational lensing by known massive clusters was 
exploited by Metcalfe et al (2003) to go deeper reaching the flux limits of 5 and 18 iiJy at 7 and 15 
/xm with ISOCAM (Metcalfe et al 2003). Elbaz et al (2002) demonstrate that MIR luminosities at 
6.75, 12 and 15 /^m are strongly correlated with the total IR luminosity (8 to 1000 /im ). They 
infer the redshift distribution of galaxies from the spectroscopically complete galaxy sample in the 
Hubble Deep Field (North) and find that the correlations hold out to z ~ 1. Chary &; Elbaz (2001) 
show that a wide variety of evolutionary models normalized to these correlations and the observed 
luminosity functions are consistent with the current CIB measurements in the mid- to far-IR. At 
present about 15-20% of the bolometric luminosity is emitted in mid- to far-IR, but it is likely that 
this fraction increases toward early times. Serjeant et al (2001) find good agreement for 90 
galaxies with the local luminosity function from other surveys and evidence for pure luminosity 
evolution at the rate oi S (x (0.98 ± 0.34) lg(l -|- z). 

The FIRBACK ISO survey (Puget et al 1999) was conducted with the ISOPHOT instrument 
at 175 /xm and reached the confusion levels of 45 mJy. It provided a catalog of almost 200 galaxies 
(at 3-sigma level) down to fluxes ~ 180 mJy (Dole et al 2001). The survey covered four square 
degrees in three high Galactic latitude regions. The final catalog of 106 sources (4-sigma) between 
180 mJy and 2.4 Jy was > 85 % complete. The observed galaxy counts at 170 fim require strong 
evolution at flux levels fainter than 500 mJy. Dole et al (2001) estimate that galaxies out to the flux 
limit of 135 mJy account for about 5% of the CIB. Analysis of the first Spitzer MIPS observations 
is consistent with the FIRBACK results (Dole et al 2004). 

Significant progress in understanding the sources contributing to the far-IR CIB has been made 
with the discovery and subsequent studies of the sub-mm galaxies with the SCUBA instrument at 
the James Clerk Maxwell Telescope. ^ SCUBA is an array of 37 bolometers at 850 /xm and 91 
at 450 /im and can observe at both wavelengths simultaneously with angular resolution of about 
15"; the shorter wavelength is more sensitive to atmospheric conditions. Several hundred SCUBA 
sources have by now been detected (Blain et al 2002) and they are mostly high-z galaxies. Redshifts 
of the galaxies have to be determined from optical observations. Because resolution of the SCUBA 
instrument is relatively low and the galaxies are at high z, and therefore faint in optical bands, 
redshift determination is difficult. Hughes et al (1998) identified optical starburst counterparts 



http://www.jach.hawaii.edu/JACpublic/JCMT/Continuum_observing/SCUBA/index.html 
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for the sample of five SCUBA sources brighter than 2 mJy, four of which have z between 4 and 
5. Scott et al (2000) detected a sample of ten FIRBACK galaxies which are bright at 170 ;um 
also at the SCUBA wavelengths with fluxes ^ 10 mJy. Sawicki & Webb (2004) detect another ten 
850 /im SCUBA sources brighter than ~ lOmJy in the general area of the Spitzer First Look 
Survey. Chapman et al (2003) have obtained spectroscopic redshifts of ten representative sub-mm 
galaxies, using the VLA to identify the source positions and obtaining the redshifts with the optical 
spectroscopy from the Keck telescope. They derived a median redshift of 2.4 for the sample and 
the space densities of (3.3 ±2.3) x 10-^ (6.5 ±2.5) x 10"^ (2.4 ± 1.2) x lO'^ Mpc'^ at redshift bins 
z = [0.5 — 1.2], [1.8 — 2.8], [2.8 — 4] respectively. They note that the median redshift range coincides 
with the peak quasar activity suggesting a close relationship between the two. The space density of 
the Chapman et al SCUBA sample increases strongly toward higher z suggesting that they make 
an important component of the SFR at 2; > 2. 

Confusion is a major problem for identifying deep galaxies in counts at sub-mm wavelengths. 
It can be reduced with gravitational lensing by massive clusters which amplifies the flux of galaxies 
and at the same time increases their apparent separation on the sky (Blain 1997). This has been 
successfully applied to the fields of known massive clusters to reach flux limits of below 1 mJy at 
850 /im (Small et al 1997, Blain et al 1999, Cowie et al 2002). The cumulative flux from the galaxy 
counts out to ~ 0.5 — 1 mJy adds up to 0.5 ± 0.2 nWm~^sr~^ accounting for all or most of the 
CIB at that wavelength. Serjeant et al (2004) account for most of the CIB at 450 /xm coming 
from the SCUBA 8 mJy survey maps of Scott et al (2002). These were combined with a Spitzer 
identification and integrating the sub-mm fluxes they find that the galaxies contribute most of the 
measured CIB at 450 /im , but only a small fraction of the CIB at 850 /xm . 

The surface density of SCUBA galaxies is quite high, ~ 300 deg~^ and, in principle, one 

can begin to study their clustering properties. Blain et al (2004) obtained accurate positions for 
a sample of distant SCUBA galaxies. The sample was followed up with Keck spectroscopy to 
determine the redshifts of 73 galaxies, constructing a substantially complete (~ 70 %) redshift 
distribution. They find strong clustering for the sample and, assuming the two-point correlation 
function ^ = {r/ro)~'^ with 7 = 1.8, measure that the SCUBA galaxies at z ~ 2 — 3 have a 
correlation length of rg — (6.9 ± 2.1)/i~^Mpc. For comparison the present day blue galaxies have 
ro = 5.5/i^^Mpc (Maddox et al 1990), and the Lyman break galaxies, that predominantly lie at 
higher z, have clustering length of tq ~ 4/i^^Mpc (Porciani & Giavilsco 2002). This should be 
reflected in larger CIB fluctuations at the sub-mm wavelengths. 

6.4. Contribution from quasars/AGNs 

A different, from stars, source of energy release is accretion onto massive and supermassive 
black holes. The latter manifests itself in the emission from AGNs and quasars, the bulk of whose 
luminosity is produced by accretion onto black holes, not nucleosynthesis. Severgnini (2000) have 
compared an X-ray sample of the [2-10] KeV and 850 /im SCUBA galax;y sources which resolve 
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most of the backgrounds in the two bands. They find the ratio of the 850 //m to [2-10] KeV fluxes 
for the sources much smaller than the value observed for the two backgrounds. The [2-10] KeV 
galaxies brighter than 10^^^ erg/cm^ make up 75 % of the X-ray background in this band, but 
contribute ^7 % to the sub-mm background. Barger et al (2004) have constructed an optical and 
NIR catalog of quasasrs and AGNs from the Chandra Deep Field out to z ~ 3. The catalogue 
currently consists of several hundred galaxies and extends to ttiab ~ 20.5 in K band, but in the 
future could pave a way to direct measurements of the AGN contribution to the NIR CIB levels. 

The current thinking is that the contribution of the active galactic nuclei and quasars to the 
CIB is most likely small. Its precise value depends on bolometric correction, the details of the 
luminosity function and its evolution at early cosmic times. Madau &i Pozzetti (2000) discuss the 
contribution from these sources to the total EBL. They estimate the total mean mass density of 
the quasar remnants today to be pBH — (3 zb 2) x IO^/iMqMpc""^ which should have contributed 
^ ^epBHC^((l + z)~^) by accretion at redshift z to the total bolometric EBL flux. Assuming then 
that these sources radiated with an average efficiency of e ~ 6% corresponding to standard disk 
accretion for a Schwarzschild black hole, they argue that, unless dust-obscured accretion on to 
supermassive black holes results in a much larger efficiency, QSOs peaking at z ~ 2, as suggested 
by observations, are expected to produce no more than 10-20% of the total EBL and a still smaller 
fraction of that will contribute to the total CIB. 

Malkan & Stcckcr (1998) estimated the levels of the CIB between 2 and 300 fim using em- 
pirically constructed SED spectra and the observed correlations between the near-IR and mid-IR 
galaxy luminosities for various galaxy types including Type 1 and 2 Seyfert galaxies. They also 
concluded that the contribution from Seyferts to the CIB is less than 10%. Lagache et al (2003) 
discuss the AGN contribution to the CIB from ~ 10 /xm to sub-mm wavelengths and, assuming 
that the black holes masses which power the AGNs are similar to those measured in the HDF, 
conlcude that the AGN contribution to the CIB is relatively small compared to that of stars. 

Elbaz et al (2002) have summed the observed galaxy counts at 15 /im from the ISOCAM 
measurements and obtain the total flux at 2.4 it 0.5 nWm^^sr^^ or just below the upper CIB limits 
(^4 — 5 nWm~^sr~^ ) discussed above. Of these, they estimate that less than 20 % comes from 
AGNs at 15 /xm and less than 10 % at 140 /xm . Matute et al (2004) estimate the luminosity 
function of AGNs and its evolution from sources in the ELAIS field at 15 /xm and conclude that 
their contribution to the CIB at that wavelengths is small (~ 2 — 3 %). 

6.5. Present-day luminosity density 

Present day luminosity density is an important normalizing factor in determining the CIB 
fluxes via cqs. 2,4. It is determined by the galaxy luminosity functions and morphological types. 
Prom UV to the near-IR it is dominated by stellar emission and at mid- to far-IR by emission from 
dust. For the Schechter (1976) luminosity function, $ = $*(L/L*)~" exp(— L/L*) the luminosity 



-55- 



density, eq. 3, becomes Cu{0) = L^^^T{a + 2). Note that the parameters for the fit are not 
independent and this parametric form, while very convenient, may not give an accurate measure 
of the overall luminosity density and its uncertainties. 

Because the near-IR CIB is produced by stars and the latter emit most of its energy in the vis- 
ible bands, any inter-comparison between the near-IR CIB and the present-day galaxy populations 
must be done over the entire visible to NIR range of wavelengths. Thus in this section we discuss 
the measurements of the present-day luminosity density from the recent galaxy surveys from UV 
to NIR and their consistency across the wavelength range. 

Early measurements of the luminosity function of galaxies (Broadhurst, Ellis &: Shanks 1988, 
Marzke et al 1994, Loveday et al 1992, Gardner et al 1997) have recently been complemented by 
much more extensive (and expensive!) surveys in the B band (2dF), visible bands (SDSS) and the 
near-IR (2MASS) using the latest multi-band and detector technologies. 

In the visible, the new measurements from the SDSS data come in five filters from 0.354 /xm 
to 0.913 fim . Blanton et al (2001) have compiled a substantial catalog of over 10,000 galaxies 
to the depth of ~ 21 — 22 AB magnitudes from SDSS commissioning observations over 140 deg^. 
They computed Petrosian magnitudes over 3" apertures arguing that Petrosian magnitudes best 
reflect the total flux over all the galaxy types. Their flrst findings gave substantially higher total 
luminosity densities than earlier measurements. However, a later re-analysis of the SDSS data 
with a larger catalog (close to 150,000 galaxies) and differently treated evolutionary corrections 
lead to lower values of Cu{0) in agreement with measurements from the 2dF (CoUess et al 2001) 
and Millenium Galaxy Catalog (Liske et al 2003) surveys in b filters using Kron and isophotal 
magnitudes. 

Kochanek et al (2001) selected over 4,000 galaxies with the median redshift of 0.02 from 
the 2MASS catalog and measured the luminosity function using isophotal magnitudes K < 11.25. 
They also subdivided the sample into early and late type galaxies (using the RC3 catalog) providing 
measurements of the luminosity function for the two morphologies. Cole et al (2001) have combined 
the 2MASS and 2dF observations with over 17,000 galaxies with measured redshifts over > 600 
deg^. They computed the luminosity function using (mostly) the Kron magnitudes. The analyzes 
of both Kochanek et al (2001) and Cole et al (2001) agree within their errors. Huang et al (2003) 
analyzed a smaller catalog of 1,056 bright (K < 15) galaxies with median redshift of 0.14 over ~ 8 
deg^ and obtain a significantly higher present-day luminosity density, but claim that the luminosity 
function they measure evolves significantly with time and, hence, its estimation may be sensitive 
to proper evolutionary corrections. 

The mid- to far-IR measurements of the luminosity function come from IRAS galaxies and 
ISO surveys. IRAS galaxies are detected most efficiently at 60 jum in the mid- to far-IR data 

and this wavelength presents the best band for galaxy identification. Soifer & Neugebauer (1991) 
used the 60 fim IRAS Bright Galaxy Sample to derive complete flux-limited samples of galaxies 
at other IRAS wavelengths. From these samples they derive the luminosity functions at of galaxies 
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in the local Universe at 12, 25, 60 and 100 fim . ISO observations of the ELAIS (European Large 
Area ISO Survey) field covered an area of 12 deg2 at 15 and 90 (Oliver et al (2000) and 

resulted in the largest ISO catalog (La Franca et al 2004, Rowan-Robinson et al 2004) enabling the 
measurements of galaxy luminosity function at these wavelengths by Serjeant et al (2001, 2004) 
and Pozzi et al (2004). 

Fig. 14 summarizes the current measurements of the present-day luminosity density from the 
various surveys. Do they indicate that a substantial flux is missing from the near-IR measurements 
when normalized to the measurements in the visible bands? Such possibility has been suggested 
by Wright (2001) using very simplified modeling of galactic SED. Below we discuss this possibility 
using more realistic galaxy modeling and including a simple treatment of galaxy morphology. With 
this modeling we show that the data on the luminosity density at various wavelengths are consistent 
with each other across the entire range from UV to near-IR. 
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Fig. 14. — Open squares show measurements of the total present-day luminosity density at the various wavelengths. The 
UV measurements at 0.28 fim is from Sullivan et al (2000). The next five squares in the direction of increasing wavelength 
correspond to the SDSS measurements from Blanton et al (2001). The near-IR points in K band are taken from Cole et al 
(2001), Kochanek et al (2001) and the highest point from Huang et al (2003). The measurement in J band at 1.2 fim is 
from Cole et al (2001). The five mid- to far-IR points are the IRAS measurements of Soifer & Neugebauer (1991) and the 
90 fxui measurements from the ELAIS field ISO observations (Serjeant et al 2001, 2004). Red sets of lines correspond to the 
contribution from the early type stellar populations alone with metallicities and ages shown in Fig. 12 as described in the main 
text, green lines show the same same for late type populations and blue lines show the simple average of the two contributions. 

Generating realistic SEDs from stellar evolution models is very important in proper interpre- 
tation of the numbers on £-^{0) as Fig. 12 shows. It is interesting to evaluate the total bolometric 



-57- 



fiux contained in the luminosity density measurements. Stellar populations contribute directly to 

emission from UV to near-IR with dust contributions dominating the emissions at longer wave- 
lengths. Integrating the data on the luminosity density from Fig. 14 from UV to 2.2 /xm and 
separately over the IRAS wavelengths we get: 

>Cboi(0.2 - 2/im) = (9.8 ± 1.2) x W^^h erg/sec/Mpc^ 
>Cboi(12 - 100/xm) = (1.5 ± 0.3) x lO^^h erg/sec/Mpc^ (30) 

here we integrated over the near-IR data from Cole et al (2001). Omitting the UV contribu- 
tion which mostly comes from very young stars would give >Cboi(0-32 — 2^m) = (8.5 it 0.7) x 
lO'^^h erg/sec/Mpc^. Taking only UV and visible data gives >Cboi(0.2 - 0.8/xm) ~ (6.4 ± 0.3) x 
10^^/ierg/sec/Mpc^. 

In order to do a proper comparison between the data on JCu{0) at visible and NIR bands we 
divide stellar populations into two types: early (E) and late (L) with the modeling discussed in 
Sec. 6.2. If each type contains I^Ei or f^Lj of the critical density and has stellar populations with 
mass-to-light ratio of M/L = 'jMq/Lq, these populations will produce r2E/L(-^©/-^©)^^7E/Lk7r^ 
leading to: 

_ -^bol /o-|\ 

If we adopt the average values of Qe = 2 x W~^h~^ and JIl = 0.6 x W~^h~^ suggested by 
Fukugita, Hogan & Peebles (1997) with 7e = 5 and 7e = 1.3 from the right panel of Fig. 12, we get 
£boi,stars = 9.1 X W^^h erg/sec/Mpc^ in good agreement with eq. 30. (The old SDSS measurements 
from Blanton et al (2001) would require higher values of Oe/l, but the latter would still be within 
the total baryon density parameter Obaryon)- Thus the measurements of the bolometric luminosity 
density from UV to near-IR suggest that no significant amounts of galaxy fluxes is missing between 
one wavelength measurement to the next. 

Are then the measurements of the luminosity density across the UV to near-IR spectrum 
consistent with each other and with that expected from realistic stellar populations? In order to 
sec the consistency between the near-IR and visible bands data on Ci,{0) we use the SEDs plotted 
in Fig. 12. We further assume that all the parameters (IMF, SED, Z, etc) are independent of L. 
The total luminosity density then becomes: 

J7e _ Cx , . 

JX,E H /A,L — -ir.45, 2 ^ 

7e 7l 10^* /i^ — 

sec Mpc-^ 

Red lines show the least squares fit to eq. 32 for each {Z, t) if only galaxies with the Early type 
stellar populations were present. This leads to Oe^^ ~ (4.5 — 7) x 10~^ for the range of {Z,t). 
Green lines show the least squares fit to eq. 32 for each {Z, t) if only galaxies with the Late type 
stellar populations were present. This leads to Clj^h"^ ~ (1.9 — 2.5) x 10~^. These numbers are in 
broad agreement with the above values on these parameters. It is clear that the luminosity density 
at visible and UV wavelengths reflects mainly emission from late type populations, whereas in the 
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near-IR both early type populations can dominate. In principle, equation 32 can be solved for both 

and r^L) but the answer will be too reflective of the assumptions made in constructing galactic 
SEDs to be useful. Instead we prefer fitting the data in Fig. 14 with a priori values of 0,^ and Ol- 
The blue lines show an example of simple mean of the green and red lines for each {Z,t). 

In conclusion, we note again that the answers one gets in this way are not unique as the re- 
alistic SEDs depends on: 1) galaxy metallicity, Z, and its dependence on other galaxy properties, 
2) the details of the IMF, particularly in the poorly measured stellar mass range and its possible 
dependence on other galaxy properties, 3) galaxy morphology mixes, 4) galaxy luminosity function 
for each galaxy type, 5) galaxy ages for each morphological type, luminosity/mass etc. If a sta- 
tistically significant discrepancy between the measurements of Ci,{0) is found between the various 
bands and surveys, it is probably indicative of the variations in any or all of the above points, but 
we have shown that with realistic (i.e. based on synthetic stellar evolution models) modeling of 
galactic SEDs the measurements can already be fit reasonably well. 

6.6. Deep galaxy counts 

The total flux from galaxies measured in deep count surveys gives another measure of the CIB, 
or more precisely the contribution to the CIB from the known sources. This helps to identify the 
possible CIB excess and in what cosmological populations and at what times it might arise. The 
situation from the current set of measurements from NIR to sub-mm bands is discussed below. 



6.6.1. Near-IR 

In the near-IR the galaxy counts have been measured now to fairly faint limits in J,H,K bands 
with observations coming from ground and Hubble Space Telescope observations. The integrated 
flux of the counts saturates at the levels shown in Table 5 (Gardner 1996, Kashlinsky & Odenwald 
2000, Cambresy et al 2000, Madau & Pozzetti 2000, Fazio et al 2004). Madau & Pozzetti (2000) 
show that the contribution from the HDF galaxies to total EBL from visible to near-IR saturates 
around AB magnitude mAB ~ 20 in all bands and discuss that this is unlikely to result from 
underestimating the abundance of distant sources due to reddening (as their Lyman break shifts 
into progressively longer wavelength bands) or absorption. 

Fig. 15 shows the contributions to the total CIB from galaxies at J,H,K bands. The data come 
from the compilations summarized in Madau & Pozzetti (2000) and Pozzetti et al (1996, 1998) and 
the Subaru Deep Field data (Maihara et al 2001); see caption for the figure. The upper panels 
show dF/dm vs magnitude and the lower panels show the cumulative flux from galaxies brighter 
than magnitude m. The fluxes saturate at mAB ~ 20 — 23. The asymptotic values of the fluxes 
from ordinary galaxies compared to the observed CIB levels are shown in Fig. 9. 
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Fig. 15. Cumulative flux in nWm^^sr^^ contributed by galaxies from a narrow dm magnitude bin is shown vs the Vega 

magnitude from galaxies from deep galaxy surveys in J, H, K. The data are as follows; 

J band: Crosses from Bershady et al (1998), asterisks from Saracco et al (1999), diamonds from Chester et al (1998) and 
triangles from Pozetti et al (1996, 1998) and Madau & Pozzetti (2000). 

H band: Crosses are from Yan et al (1998), asterisks and diamonds from Teplitz et al (1998), triangles from Chester et al 

(1998), and squares from Pozetti ct al (1996, 1998) and Madau & Pozzetti (2000). 

K band: crosses arc from Chester ct al (1998), asterisks from McLcod ct al (1995), small open diamonds from Djorgovski et 
al (1995), small triangles from Mobasher et al (1986), small and large open squares, large triangles and large open diamonds 
from Gardner et al (1993), x-signs from Glazebrook et al (1994), filled diamonds from Soifer et al (1994), filled large circles 
from Pozetti et al (1996, 1998) and Madau & Pozzetti (2000), and filled small circles from Maihara et al (2001). 

ISO galaxy counts at 6.75 //m reach ~ 40//Jy in the HDF (Ohver et al 1997) and Lockman 
hole regions (Taniguchi et al 1997) and Sato et al (2003) reach sources as faint as ~ 10/xJ in another 
region. At these wavelengths the counts have now been complemented by deeper measurements 
with the Spitzer IRAC instrument. 

Fig. 16 shows the contributions to the CIB flux from faint galaxies at 3.6, 4.5, 5.8 and 8 
jum from three recent Spitzer IRAC surveys reported by Fazio et al (2004). The surveys covered 

3 independent fields referred to as QSO field (the deepest with 9.2 hours total exposure), the 
EGS (Extended Groth Strip) field at a higher Galactic latitude and the widest and the shallowest 
Bootes region survey. They separated stars from galaxies and corrected for incompleteness in the 
faintest counts (mostly galaxies). At 3.6 and 4.5 fim. the total flux seems to saturate at the levels 
of shown in Table 5. At 5.8 and 8 /xm the saturation of the cumulative flux from observed galaxies 
is not as clear, although a case can be made that the ceiling of the total flux contribution to the 
CIB has been reached as well. At any rate we interpret the total fluxes at 5.8 and 8 /Ltm as lower 
limits on the contribution from ordinary galaxies to the CIB and plot them with arrows in Fig. 9. 
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The numbers from the Spitzer IRAC surveys shown in Table 5 are direct summation of the fluxes 

from upper panels in Fig. 16. When integrating the counts weighted according to uncertainties, 
Fazio et al (2004) find total fluxes of 5.4, 3.5, 3.6, and 2.6 nWm-^sr"! at 3.6, 4.5, 5.8 and 8 /xm 
respectively. 
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Fig. 16. — Cumulative (bottom set of panels) and differential (top) flxix distribution in nWm^^sr"^ vs the Vega magnitude 
for galaxies in three different surveys from the Spitzer IRAC galaxy surveys by Fazio et al (2004); dots with errors correspond 
to their Bootes field, diamonds to their EOS data and triangles to the QSO1700 field observations. AB magnitudes are shown 
on the upper horizontal axis. 



6.6.2. mid-IR 

Elbaz ct al (2002) used ISOCAM observations of almost 1,000 galaxies at 15 ^m to estimate 
the contribution to the CIB at that wavelength from sources down to 50 /xJy. They obtain that these 
galaxies produce 2.4 it 0.5 nWm~^sr~^ and their contribution may saturate at fluxes $30 — 50;uJy. 

Spitzer/MIPS 24 /im channel has best mid-IR resolution and its confusion limit corresponds 
to fainter fluxes. The MIPS FWHM at 24 /Ltm corresponds to 6". Papovich et al (2004) present 
the number counts of ~ 5 x 10^ sources from the 24 /xm Spitzer deep surveys. They show that the 
counts probe a previously undetected population of very luminous galaxies at high z. The data were 
obtained in five fields from Spitzer characterization and GTO (Guaranteed Time Observations) 
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observations. The largest and shallowest of the fields (Bootes) subtended ca. 9 deg^ with an 

exposure of 87 sees per pixel and the smallest and deepest (ELAIS) was ca. 130 arcmin^ exposed 
for just under one hour in each pixel. The counts out to 60 fim give ~ 2 nWm~^sr~^ and Papovich 
et al (2004) estimate that extrapolating to sources below that limit gives around 3 nWm~^sr~^ in 
total. 
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Fig. 17. — Left: diflferential flux distribution in nWm~^sr~i from the MIPS counts at 24 //m (Papovich et al 2004) 
Right: cumulative flux contribution as function of the source flux from galaxies shown in the left panel. Shaded area corresponds 
to the total flux with its uncertainty including the extrapolation from sources fainter than 60/iJy (Papovich et al 2004). 



6.6.3. Far-IR and sub-mm 

The number counts at mid- and far-IR indicate that the IR-luminous galaxies evolved more 
rapidly than their optical counterparts and make a substantial contribution to the star formation 
at higher z (Elbaz et al 1999). 

A deep survey (FIRBACK) was performed by ISOPHOT at 170 /xm (Lemke et al 1996) out to 
the depth of 120 mJy where galaxy counts can no longer be fitted with a Euclidean slope. This is 
confirmed with the MIPS observations of the Chandra Deep Field (Dole et al 2004) which indicate 
a lower limit on the CIB at 170 fim of 1.4 nWm~^sr~^ corresponding to about 10% of the CIB. 
Dole et al (2004) used early observations with the Spitzer MIPS instrument to measure galaxy 
counts at 70 and 160 /xm down to 15 and 50 mJy respectively. Their counts are consistent with the 
60 urn IRAS counts (Lonsdale et al 1990) and the FIRBACK 170 /xm survey (Dole et al 2001). 
Prom evolutionary modeling they suggest that most of the observed faintest sources lie at z ~ 0.7 



-62- 



with a tail out to z ~ 2. Integrating the flux from the observed sources leads to the total flux of 
0.95 and 1.4 nWm~^sr~^ at 70 and 1 60 fim respectively. As Fig. 9 shows this corresponds to a 
small part of the total CIB levels at these wavelengths, so most of the CIB at 70 and 160 /xm must 
come from still fainter (and more distant) sources. 

The FIR background can be resolved into individual sources with detections by the SCUBA 
instrument at the James Maxwell Telescope (Small et al 1997, Blain et al 2002). The total flux at 
850 /Ltm from the the SCUBA sources is estimated to be 0.5 ± 0.2 nWm^^sr"^ (Small et al 1997). 

Serjeant et al (2004) presented statistical detections of galaxies in the Spitzer Early Release 
Observations through a stacking analysis of their reanalyzed SCUBA 8 mJy survey maps of Scott et 
al (2002) combined with a Spitzer identification of their positions. Integrating the sub-mm fiuxes of 
the Spitzer populations they find that the 5.8 /xm galaxies contribute only 0.12 ± 0.05 nWm^^sr"-*^ 
at 850 ^m , but at the same time contribute 2.4 it 0.7 nWm~^sr~^ at 450 fiui , or almost the 
entire CIB at that wavelength. This contribution is shown with filled circle at 450 /xm in Fig. 9. 

6.7. CIB fluctuations from clustering of ordinary galaxies 

In surveys where the beam is small, galaxies can be resolved and removed down to some limiting 
magnitude allowing, in principle, to probe the CIB from fainter and typically more distant systems. 
At z=l, the angular scale of 1'.' subtends comoving scale 15/i~^Kpc, for the WMAP cosmological 
parameters, so the identified galaxies can be excised almost completely in surveys with arcsec scale 
resolution. In practice, the precise value of the limiting magnitude depends on the instrument noise 
levels, foreground emissions, etc. 

Small beam also means that the shot-noise component of the fiuctuations discussed in Sec. 
3.2.2 may be important. The left panel of Fig. 18 plots the shot noise power spectrum (left vertical 
axis) and the value of (Tg^ (right vertical axis) given by eqs. 16,15 vs the AB magnitude down 
to which galaxies have been removed. The right panel shows dsn vs the wavelength for galaxies 
fainter than ruAB = 20 (open circles), 22 (filled circles), and 25 (filled triangles). The shot noise 
was evaluated from galaxy counts data summarized in Sec. 6.6. In the figure we assumed the beam 
of o^beam = 2.1 X 10~^° sr Or, in case of a spherical beam, a radius of ~ 1.6". We did not attempt to 
calculate the (probably substantial) statistical uncertainties in the shot noise components because 
the various systematic, cosmic variance, etc effects in the different galaxy surveys are difficult to 
quantify (which would probably make such error bars misleading), but the figure should give a 
reliable approximation for the magnitude of the expected shot noise. 

The left panel of Fig. 19 shows the amplitude of the detected CIB fluctuation at = 1 arcsec 
from the analysis of the deep 2MASS data in J, H, K bands (Kashlinsky et al 2002, Odenwald et 
al 2003). It is plotted vs the AB magnitude above which galaxies galaxies were clipped out for the 
three 2MASS bands. Although the present deep galaxy counts data have potentially appreciable 
uncertainties for evaluating the shot noise amplitude, comparison with Fig. 18 shows that the 
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Fig. 18. Left: Shot noise power spectrum (left vertical axis), evaluated according to eq. 16, from galaxies fainter than 

the AB magnitude shown in the horizontal axis. It was computed for galaxy counts summarized in Sec.6.6. Left vertical axis 
shows the corresponding value of <Tsn for a beam of ~ 1.6" in radius. The lines show the shot noise from HDF galaxy counts 

by Madau & Pozzctti (2000): solid, dotted and dashed lines correspond to to J, H and K bands respectively. Symbols show the 
shot noise from galaxies in the IRAC Spitzcr counts data from Fazio et al (2004): plus signs, asterisks, diamonds and triangles 
correspond 3.6, 4.5, 5.8 and 8 ^m respectively. Right: osn vs A evaluated from existing deep galaxy counts data. Open circles 
correspond to galaxies fainter than rriAB = 20, filled circles to rriAB > 22 and filled triangles to rriAB > 25. 

detected CIB fluctuations are substantially higher than the shot-noise fluctuation from galaxies 
remaining in the field. The largest angular scales probed by the analysis were n/q ~ 1.5' which 
at z=l corresponds to ^l^^^Mpc. Hence the data are probing the range of scales where galaxy 
clustering is in the non-linear regime. The lines show the amplitude for passively evolving galaxies 
with synthetic spectra from Jimenez &: Kashlinsky (1999) for J and K bands. They are normalized 
to the observations of the galaxy luminosity densities in Fig. 14, the observed power spectrum of 
galaxy clustering, include the shot-noise component from Fig. 18, and assume no biasing. The high 
numbers for the near-IR fluctuations are consistent with other findings (Kashlinsky &; Odenwald 
2000, Matsumoto et al 2000, 2003) after accounting for the beam difference and the fact that in 
the large-beam DIRBE and IRTS studies no galaxies were removed. 

The figure illustrates that the detected CIB anisotropy levels are large compared to simple 
no-evolution assumptions. They may imply either significant evolution of galaxy populations at 
early cosmic times, or strong biasing of faint (more distant) galaxy population. The latter would 
then have to be wavelength dependent because the discrepancy is different in the various bands. 
Alternatively, Maggliochetti et al (2002) ascribe the difference to contribution from Population III 
stars at z ~ 10 — 20; this contribution is discussed in the next section. 

The slope of the detected power spectrum of the CIB from Kashlinsky et al (2002) depends on 
the clipping magnitude cutoff as shown in the right panel of Fig. 19. If progressively fainter galaxies 
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Fig. 19. — Left: Amplitude of the RMS CIB flux fluctuation at = 1 arcsec from Fig. 11 vs AB clipping magnitude 
beyond which galaxies were excised from the field. Open triangles, open squares and crosses correspond to the 2MASS J,H,K 
bands respectively. Dashed and solid lines show the fluctuation from passively evolving galaxies, including the shot-noise 

component, for J and K bands respectively. Right: Power law index n for the CIB power spectrum P(q) oc q" from the deep 
2MASS data shown in Fig. 11 is plotted vs the limiting clipping magnitude. Same symbol notation as in the left panel. 



are removed the slope of the power spectrum flattens. This is consistent with fainter galaxies being 
at higher z when the clustering pattern was less evolved and the slope reflects evolution of galaxy 
clustering with time. The slope of this dependence is similar in all three bands. The power spectrum 

of the present-day galaxy clustering on non-linear scales has the spectral index n ~ —1.3 (Groth 
&: Peebles 1977, Maddox et al 1990), so it would appear that the deep 2MASS J band patch with 
the lowest rricut probes galaxies out to smaller z. In CDM models the non-linear power spectrum 
of galaxy clustering evolves toward steeper slope (higher |n|) at low z consistent with the trend in 
Fig. 19b. 



6.8. Cumulative flux from galaxy counts vs CIB measurements 

At many IR wavelengths the counts from ordinary galaxies are deep enough to make a 'reason- 
able' guess about their total contribution to the CIB. This would indicate whether other sources 
of radiation existed that produced the observed CIB flux and its fluctuations. Table 5 summarizes 
the total fluxes observed directly by the deepest currently available galaxy populations and one can 
compare them with the CIB measurements summarized in Fig. 9. 

At mid-IR there are no direct CIB measurements to compare, but indirect upper limits are 
close to the lower limits from the counts so probably no surprises are expected. At far-IR most 
of the CIB has been resolved at 450 and 850 /xm with the flux from the observed counts not yet 
saturating. At NIR the counts from the observed 'ordinary' galaxy populations saturate at levels a 
factor of ~ 2 — 3 lower than the claimed levels of the mean CIB measurements and the fluctuations 
analysis. Furthermore, the CIB fluctuations results from the deep 2MASS data suggest that the 



Table 5. Integrated CIB flux from observed NIR galaxy counts 



A 

1.25 lira ^ 
1.65 /im 1 
2.2 /im 1 
3.6 /tm 2 
4.5 /tm 2 
5.8 /jm 2 
8 /im 2 


Flux (nWm-2sr-i ) 
q 71+3.00 
^- ' ^-1.90 

7.92l?;^t 
5.27 ±1.02 
3.95 ± 0.77 

> 2.73 ±0.22 

> 2.46 ± 0.21 


Comments" 

Saturates at m-AB 24 — 25 

Saturates at mAB 24 — 25 

Saturates at rriAB 22 — 23 

Saturates at rriAB jl 22 — 23 

Saturates at mAB ^22 — 23 
Possibly saturates at mAB ^ 22 — 23, but data too uncertain. 
Possibly saturates at mAB 22 — 23, but data too uncertain. 


1.25 /tm < A < 8 /im 
1.25 /tm < A < 4.5 /tm 


11.4 ±4.6 
9.7 ±3.2 


NIR (bolometric) flux 


15 /(m ^ 
24 /tm 4 
70 /tm 5 
170 /im 5 
450 /tm 6 
850 /tm 


2.4 ±0.5 

2 7+1-1 
^•'-0.7 

> 0.95 

> 1.4 
2.4 ±0.7 
0.5 ±0.2 


Saturates at 5< 0.05 mjy 
Saturates at 5 < 0.05 mJy. ^ 
Does not saturate out to 5 > 15 mJy. 
Does not saturate out to 5 > 180 mJy. 
Sub-mm fluxes of Spitzer sources at 24 /tm account for most of the CIB 
Sources out to 5 ~ 0.5 — 1 mJy account for most of the CIB 



Note. — References: i Madau & Pozzetti (2000), 2 Fazio et al (2004), ^ Elbaz et al (2002), * Papovich et al (2004), ^ Dole 
et al (2004), ^ Serjeant et al (2004), Blain et al (1999) 
° Saturation point is defined as the magnitude where ^.^.^^^ < 0.2 

CIB not measured directly at that wavelength 

Counts are not measured to sufficiently faint limits. 
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excess arises in galaxy populations with mAB~21 and likely comes from early times. 

Fig. 20 plots the NIR CIB excess left over after galaxy counts contributions from Table 5 are 
subtracted from the observed CIB levels shown in Fig. 9. The last point is plotted at 4 fim , 
where IRTS data are still available and galaxy counts were extrapolated to this wavelengths using 
the numbers in Table 5. Integrating over the points in the figure leads to the total bolometric flux 
for the near-IR CIB excess between 1 and 4 jim of 29.4 ± 13.0 nWm~^sr~^ . Ignoring the last 
point (open circle) in the figures leads to 26.4 ± 12.2 nWm^^sr"-*^ . 

Zodiacal light at these wavelengths would come from reflected Solar emission and observations 
suggest a slope oc A^^ (Kelsall et al 1998 and Fig. 4). Assuming that zodiacal light errors are 
proportional to the zodiacal flux, the excess may be attributable to the residual zodi emission 
if the points at 1.25 and 4 ij,m can be discarded. Comparison with the zodiacal light surface 
brightness shown in Fig. 4 suggests that in this case, the relative errors in the zodiacal emission 
modeling would have to have amplitude of ^ 40-50 % (sec Kelsall et al 1998 and Fig. 4) at the 
near-IR bands. They would also have to be the same in all the different CIB measurements with 
different instruments and different zodiacal models used. Galactic cirrus emission would at these 
wavelengths contribute emission with a much smaller and a roughly flat SED (Leinert et al 1998). 

CIB excess: Fboi(l-4 yu.m) = 29±13 nW/rnVsr 
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Fig. 20. CIB excess at near-IR: filled circles correspond to CIB measurements minus the cumulative fiux from the observed 

galaxy populations. As discussed in the text, the various measurements are consistent, but for the purposes of this plot, the 
CIB mean levels were taken from Cambresy et al (2001) at J band, and from the IRTS measurements at 1.63 to 4 fira . Filled 
circles show wavelengths where galaxy counts can be summed directly. Open circle corresponds to extrapolation of the observed 
fluxes from galaxy counts given in Table 5 to A = 4 jum , the maximal wavelengths of the IRTS measurements. 
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7. Population III 

The epoch of the first stars, when Population III stars formed, is now emerging as the next 
cosmological frontier. It is not clear what these stars' properties were, when they formed or how 
long their era lasted before leading to the stars and galaxies we see today. On the other hand, if the 
excess ncar-IR CIB is real, as many independent measurements suggest, it may well be attributable 
to emissions by these stars and would make the CIB and its structure a unique probe directly into 
the epoch and efficiency of Population III formation and evolution. This would provide a powerful 
application of the CIB-related science to one of the most outstanding remaining questions in the 
standard cosmological picture. 

Recent measurements of quasar spectra reveal the Gunn-Peterson trough at Arest < 1216A 
at z>6 — 7 indicating the location of the reionization epoch (Becker et al 2001, Djorgovski et al 
2001). These quasars already show the presence of metals suggesting that the first metal free stars 
had to form at still earlier epochs and enrich the IGM by z ~ 7. At epochs corresponding to 
z^7 the Universe had to contain significant amounts of neutral hydrogen leading to absorption of 
all radiation at wavelengths shorter than the Lyman limit. Interestingly, the J band filter probes 
wavelengths longer than the Lyman limit at zUlA and the K band at 20, so much of the near-IR 
CIB can probe the Population III era even if reionization has not yet occurred. 

7.1. What were the first stars? 

Population III stars are the first stars to form in the metal-free Universe and so far remained an 
entirely hypothetical (and theoretical) class of objects since introduced by Rees in the late seventies 
(Rees 1980 and references cited therein). Much of what can be attributed to consequences of that 
early and possibly brief era depends on the the mass Population III stars had and we start this 
section with a brief discussion of the theoretical spectrum of possibilities and the current ideas. 

Should the first stars have been different from modern day stellar objects? We know from 
observations that star formation takes places in gas clouds much more massive than Mchandra = 
(^^)^/^mp ~ 1.44M0, a mass-scale of a typical self-gravitating object supported by nuclear fusion, 
and must be caused by gravitational collapse (stellar material is much more dense than the clouds 
in which they form) and fragmentation (stars are much smaller in mass than the parental clouds). 
Thus the gravitational collapse processes must lead to final fragments that are just in the right 
mass range to support nuclear burning. 

Hoyle (1953) provided a plausible and elegant answer to this that became known as "opacity 
limited fragmentation theory" . It goes as follows: the cooling time is much shorter in galactic mass 
gas clouds than the dynamical collapse time and the ratio of the two decreases with increasing 
gas density (Rees &: Ostriker 1977). Gravitational collapse thus proceeds isothcrmally leading to a 
decrease in the Jeans mass inside the cloud, Mjeans T'^/^p"^/^. The cloud becomes susceptible 
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to fragmentation into masses of ^ Mjeans, which themselves collapse and fragment into smaller 
pieces. This process continues until the density increases enough so that the optical depth across 
the fragment is sufficiently high to trap the radiation released by the collapse. This requires 
r oc Kp^/^M^/^ of order unity (for absorption; for scattering this happens at higher r). For clouds 
of solar metallicity the collapse proceeds at T ~ lOK and fragmentation stops at ~ O.OIMq (Low & 
Lynden-Bell 1977, Silk 1977). For metal free gas clouds the collapse proceeds at higher T (> lO^K 
in the absence of hydrogen molecules, Dalgarno &: McCray 1972). Thus the early thinking in the 
pre-CDM era was that the first stars would have to be much more massive. On the other hand, 
Rees (1976) has shown, from very general arguments of maximal cooling efficiency, that the minimal 
mass is roughly independent of T and is Mmin Mc]^{kBT /rripC^)^^* . I.e. in metal free gas the 
cooling is very inefficient so T is high, but at the same time - because opacity is generated by the 
same coolants - r is smaller and the fragmentation stops at higher densities and smaller Mjeans- 

However, a consensus based on recent numerical investigations within the framework of the 
ACDM models is now emerging that fragmentation of the first collapsing clouds at redshifts ~ 
10 — 30 was very inefficient and that the first metal-free stars after all were very massive objects 
with mass IOOMq ((Abel et al 2002; Bromm et al 1999), see recent review by Bromm & Larson 
2004). Such stars would form in small mini-halos and live only a few million years, much less 
than the age of the Universe (~ 2 x 10^ years at z = 20), making their direct detection still more 
difficult. On the other hand the net radiation produced by these massive stars could give substantial 
contributions to the total diffuse background light ((Rees 1978)) and since their light is red-shifted 
much of that contribution will be today in the infrared bands ((Bond et al 1986)). This could make 
them significant contributors to the near-IR CIB. 

Such massive, metal-free stars will be dominated by radiation pressure and would radiate 
close to the Eddington limit: L ~ ^Edd = IZI^^isiiTV/ ~ 1.3 x 10'^^M/Mq erg/sec, where ax is the 
cross section due to electron (Thompson) scattering. The energy spectrum for emission from these 
stars will be quite featureless and close to that of a black body at T ~ 10^'^~^K (Schaerer 2003, 
Tumlinson et al 2003). Unlike their metal rich counterparts Population III stars are not expected to 
have significant mass-loss during their lifetime (Baraffe et al 2001). If sufficiently massive (^ 24OM0, 
Bond ct al 1984, Heger et al 2003) such stars would also avoid SN explosions and collapse directly 
to black holes in which case their numbers can significantly exceed that required to produce the 
metallicities observed in Population II stars. The lifetime of these stars will be independent of 
mass: i/, ~ eMc^/L ~ 3 x 10^ years, where e = 0.007 is the efficiency of the hydrogen burning. 
These numbers are in good agreement with detailed computations (e.g. Schaerer 2002). 

From the WMAP large-scale polarization measurements (Kogut et al 2003) it is known that the 
Universe's optical depth to recombination is r ~ 0.2 requiring the re-ionization to occur by ~ 20. 
Numerical calculations suggest that when a Population III star forms it is surrounded by a gaseous 
nebula from the gas of the host halo which was not incorporated into the stars (Bromm et al 1999). 
The nebula and the IGM would remain neutral in the absence of Population III. In general, the 
resultant SED would thus be made up of three components: 1) direct stellar emission not absorbed 
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by Lyman continuum; 2) Lyman emission of the absorbed part of stellar SED, 3) free-free emission 
from the IGM. Consequently, Santos et al (2002) consider two extremes for the reprocessing of the 
ionizing radiation from Population III: 1) Each Population III star is surrounded by a dense nebula 
with all the reprocessing of ionizing radiation taking place there, and 2) the nebula is optically 
thin to the ionizing photons and all the reprocessing takes place in the IGM. Fig. 21 shows the 
SED from Santos et al (2002) produced by a 1000 Mq metal-free star for the two extreme cases of 
photon propagation. Superimposed are also the J and K filters redshifted to z = 20 and 10. 
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Fig. 21. — An example of the Population III SED produced by a 1000 Mq metals-free star from Santos et al (2002). Light 
and dark shaded area show the range probed respectively by the J and K band DIRBE filter at z = 20 and 10. Dashed-dotted 

line shows the spectrum of a 300 Mq metal free star; thin solid line shows the part of that SED loiigward of the Lyman limit. 
Dotted lines show the free-free emission from the surrounding medium and dashed lines show the nebulae emission; thick and 
thin lines of each type correspond to the two extremes photon propagation considered by Santos et al (2002). Thick solid lines 
correspond to the total SED produced by the Population III emission in the two extreme cases. 

There are several intuitive reasons why Population III produce significant both CIB levels and 
fluctuations: 

• Each unit mass of Population III (if made of massive stars) emits > 10^ more light than 
normal stars. 

• A higher fraction of the total luminosity would be redshifted into the NIR bands from hot 
objects at these redshifts than from z ~ 2 — 3. E.g. the K band redshifted to z ~ 20 corresponds 
to rest wavelength of 1000 A, roughly the Lyman limit wavelength. 

• Massive stars (assuming Population III were indeed massive stars) radiate with a higher 
mean radiative efficiency, e = 0.007, than the present day stellar populations. 

• Because their era was presumably brief. Population III epochs contain less projected volume 
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than the ordinary galaxy populations spanning the epochs of Population I and I stars. Hence, the 
larger relative fluctuations. 

• Biasing is higher for Population III because they form out of rarer regions which leads to the 
amplified correlations. 



7.2. Isotropic component of CIB 

It was noted by several authors that the near-IR CIB signals such as e.g. detected in the 
2MASS long integration data (Kashlinsky et al 2003) may come from the Population III stars 
(Magliocchetti, Salvaterra &: Ferrara 2003, Santos, Bromm k. Kamionkowski 2002, Salvaterra & 
Ferrara 2003, Kashlinsky et al 1999, Cooray et al 2004, Kashlinsky et al 2004). The argument 
can be shown to be (almost) model-independent provided Population III stars were indeed very 
massive. 

Our discussion in Sections 7.2 through 7.4 follows Kashlinsky et al (2004): 

For a flat Universe the co-moving volume per unit solid angle contained in the cosmic time 
interval dt is dV = c(l -|- z)~^d\dt, where dL is the luminosity distance. Each Population III star 
will produce flux Lhi,i{l+z)lAiTdPi^, where is the fraction of the total energy spectrum emitted per 
unit frequency and v' = v{l+z) is the rest-frame frequency. (The Population III SED is normalized 
so that J hydu = 1). The co-moving mass density in these stars is J Mn{M,t)dM = ribaryongi^/*) 
where is the fraction of the total baryonic mass in the Universe locked in Population III stars 
at time t. The net flux per unit frequency from a population of such stars with mass function 
n{M)dM is given by: 



d^ jLniM,t)dM . ^ . , dV c.,..L.. 

J^I. = (1 + Z) {K') = ^(V)(^)/*Pbaryonic 



(33) 



Here {b,^) = J Ln{M,t)b,ydM/ J Ln{M,t)dM denotes the mean Population III SED averaged over 
their initial mass function and {-^) = / Ln{M,t)dM/ J Mn{M,t)dM. For Gaussian density field 
/* ~ 5x 10~^— 3x 10~^ if on average Population III formed in 2-3 sigma regions (see later). Provided 
that L/M = constant, this result does not depend on the details of the initial mass function of 
Population III stars (cf. Fig. 5 of Salvaterra &: Ferrara 2003). (Note that for the present day stellar 
populations their mass-to-light ratio depends strongly on stellar mass and {^) is much smaller 
than the Population III value of AirGmpc/aT — 3.3 x 10^Lq/Mq leading to substantially smaller 
net fluxes). Assuming no significant mass loss (Baraffe et al 2001) during their lifetime tz for 
Population III stars, these stars will produce CIB of amplitude: 

3 1 ^ 1 f r , n . dt 



Note that in eq. 34, c^/G is the maximal luminosity that can be achieved by gravitational 
processes. It enters there because ultimately the nuclear burning of stars is caused by gravity 
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as they evolve in gravitational equilibrium with the (radiation) pressure. Eq. 34, has a simple 
meaning illustrating its relative model-independence: the cumulative bolometric flux produced 
is Lmax = c^/G ~ IO^^Lq distributed over the surface of the Hubble radius, AttR^, times the 
model dependent, but more-or-less understood, dimensionless factors. The spectral distribution of 
the produced diffuse background will be determined by (6,^) = / + z)~^hydt/ J f^,{l + z)~^dt 
i.e. the mean SED averaged over the Population III era redshifts. The value of -Z^max/47ri?|^ is 
~ 3 X 10^ nW/m^/sr so even with small values of J^baryon, e, /*, the net flux from Population III 
stars could be indeed substantial. 

The first stars had to form in the rare regions of the density field that reached the turn-around 
over-density while the bulk of the matter was still in linear regime (density contrast < 1). They 
would then turn-around and, if certain conditions are met, collapse and form the first generation of 
stars. For Gaussian density field, as is the case with the CDM models, this would specify the value 
of /* at each time. It is perhaps a bit unexpected that the WMAP polarization measurements 
indicated that first stars started forming at = 20l'lg°, but the most straightforward explanation 
is that the earliest stars formed out of very high peaks of the density field leading to smaller 

Once they turn around and begin to collapse the primordial clouds will quickly heat up to their 
virial temperature, Tyir, and only efficient cooling will enable isothermal collapse when gas pressure 
is smaller than gravity (or T < T^n) ■ The details of what determines collapse and formation of 
Population III stars are not yet clear, but molecular cooling seems critical in the initial stages of 
metal- free collapse. In general H2 cooling is effective at T^AQO K (Santos et al 2002), but some 
simulations suggest rotational H2 cooling can effectively dissipate binding energy of the cloud only 
if T > 2000 K (Miralda-Escude 2003). The numbers that follow have thus been evaluated for Tyir = 
400K and 2000K. The relation between comoving scale r, the mass contained by that radius and 
the virial temperature at the onset of the collapse is: 

In the spherical linear approximation, the mass that at some early epoch, Zi, had density 
contrast {Scoi = 1-68 times the growth factor between Zi and z) will collapse at redshift z. For 
Gaussian density field, the regions that collapse correspond to 77 = <5coi/[(<^m)]^ > 1 standard 
deviations the value of ^ crfc(?7/\/2). The value of the RMS density contrast, at 
redshift z can be evaluated given the power spectrum of the assumed model (ACDM) normalized 
to large-scale CMB data and can be estimated from Fig. 3 and cqs. 35. Population III stars had 
to form out of ~ 2 — 3 sigma rare fluctuations which for Gaussian distribution correspond to 
varying between ~ 5 x 10~^ and ~ 3 x 10~^; departures from spherical symmetry would accelerate 
the collapse thereby decreasing 77 and increasing Because 77 is a decreasing function of decreasing 
z, and ft, is a (rapidly) increasing function of decreasing rj, the average /* will be dominated by 
the late times of Population III evolution. 

Fig. 22 shows the total bolometric flux evaluated according to eq. 34 evaluated for ACDM 
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Fig. 22. — Bolometric flux produced from the Population III era lasting until the redshift Zf shown along the horizontal axis. 

Solid line assumes Population III stars form in collapsing regions with Tvir = 400 K and dashed line corresponds to Tyi^ = 400 
K. Shaded area shows the bolometric flux of the CIB excess at near-IR evaluated from Fig. 20 with its uncertainty. 

model for Tyir = 400 K (solid line) and Tyir = 2000 K plotted vs the redshift marking the end of 
the putative Population III era, Zf. The lines in the figure assume that efficiency of Population III 
formation is regulated only by parameter The prediction for the total bolometric flux produced 
by Population III is fairly robust. On the other hand, the way it is distributed among the various 
bands or its energy spectrum depends on the details of how Lyman-a ionizing photons are absorbed 
and re-emitted along the line-of-sight. Santos et al (2002) model the redistribution of the total flux 
produced from the Population III era and produce good flts to the observed near-IR CIB at J and 
K bands. Clearly, in order to account for the near-IR CIB excess measured at 1-4 fiva , the diffuse 
bolometric flux from Population III stars must exceed the observed excess; any additional flux will 
then have to be redistributed to different wavelengths. This can be achieved if the Population III 
stars were massive and their era lasted at least until Zf ~ 15. Interestingly, the existence of the 
CIB excess at J band suggests that if this excess originated from the Population III era, the latter 
should have extended to z< 14. 

This discussion suggests that Population III era could have produced CIB at levels comparable 
to those in Fig. 20, but it could also have produced substantially lower (or higher) fluxes, while 
the measurements of the CIB mean levels can be significantly affected by the systematics and be 
mistaken for the various residual errors. On the other hand, as is shown below, Population III 
stars, whose emission arises at epochs when the spatial spectrum of galaxy clustering is not yet 
evolved, should have produced a unique and measurable signature via the CIB fluctuations. It is 
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that signature, both its spectrum in the angular and energy frequency domain, that could provide 
the ultimate insight into the Population III epochs (Cooray et al 2004, Kashlinsky et al 2004). 

7.3. Contribution to anisotropies in CIB 

In order to evaluate the amplitude and spectrum of the CIB from Population III era we start 
witli the density field of the ACDM model. For reference, a 1 arcminute angular scale subtends 
comoving scale ~ 1.5/i~^Mpc between z = 6 and z = 15 for WMAP cosmological parameters and 
equation of state with w = —1. As Fig. 3 shows, on sub-arcminute scales the density field is in 
quasi-linear to non-linear regime (density contrast ^0.2) for the ACDM model. There the spectrum 
due to clustering evolution was modified significantly over that of the linear ACDM spectrum, and 
the fluctuations amplitude has increased from its primordial value especially since the effective 
spectral index on these scales for ACDM model was — 2. 

At 2; = 20 the Universe is ~ 2 x 10^ years old which is much larger than the age of the individual 
Population III stars, ~ 3 x 10^ years. If Population III era lasted for only one generation forming 
at 2;*, the CIB fluctuations will be ^cib ~ V^'^ildA^'', ^)\z*i where A(/c) is given by eq. 12. In what 
follows we define U as the time-length of the period in which Population III stars were the dominant 
luminosity sources in the Universe. 

Several points are worth noting about eq. 11 when applied to Population III era: 1) For given 
amplitude of the mean CIB levels from Population III, the value of A (A;) is inversely proportional 
to ^/t^ and thus measures the duration of the Population III era. The density perturbations grow 
with decreasing z, thus most of the contribution to the integral in eq. 10 comes from the low end 
of z over some range of b,^. At some wavelengths the overall dependence on wins out at some 
(shortest NIR) wavelengths and the longer the Population III phase, the smaller arc the relative 
fluctuations of the CIB from them. 2) In the Harrison-Zeldovich regime of the power spectrum, 
P3 oc k, one would have Sqib oc q^-^. 3) The transition to the Harrison-Zeldovich regime occurs in 
the linear regime at all relevant redshifts and happens at the co-moving scale equal to the horizon 
scale at the matter-radiation equality. All this would allow probing of the Population III era, its 
duration, and the primordial power spectrum at high redshifts on scales that are currently not 
probed well. Interestingly, short duration of the Population III will lead to smaller CIB flux, but 
larger relative fluctuations and vice versa. 

In addition to the small angular scale increase due to non-linear gravitational evolution, the 
fluctuations in the clustering distribution of Population III systems will be amplified because, within 
the framework of the ACDM model these systems had to form out of rare peaks of the primordial 
density field (Kaiser 1984, Jensen & Szalay 1986, Politzer & Wise 1996, Kashlinsky 1987, 1991, 
1998). This leads to biased (enhanced) 2-point correlation function of the Population III regions, ^5 
over that of the underlying density field, ^(r) = ^ f PACBM{k)k'^jo{kr)dk. In the ACDM model 
the Population III stars at z ~ 10 — 20 will form in regions having 77 > 1. In this limit the biasing 
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factor given by eq. 29 of Kashlinsky (1998) reduces to check: 




In the Umit of ^ — > this reduces to the more famihar form of ^^(r) ~ ':^^^if)- We adopt eq. 36 
in evaluating tlie numbers below. 

There will also be shot-noise fluctuations due to individual Population III systems entering 

— 1/2 

the beam. The relative magnitude of these fluctuations will be N^^^^^, where A'^bcam is the number 
of the Population III systems within the beam. This component may be important at very small 
angular scales, where iVbeam ~ !> and will contribute to the power spectrum: PsN = where 
n2 = cj n*(i|^(l + z)~^dt is the projected angular number density of Population III systems. The 
detection of the shot-noise component in the CIB power spectrum at small angular scales will 
give a direct measure of both the duration of the Population III era and constrain the makeup 
and masses of the Population III systems. Calculations show that unless the Population III era 
was very short, the shot noise correction to the CIB would be small on scales greater than a few 
acrseconds (Kashlinsky et al 2004). 

Because of the Lyman absorption by the surrounding matter prior to reionization, essentially 
all the flux emitted by Population III stars at rest wavelength less than 912 A will be absorbed 
by the Lyman continuum absorption in the surrounding medium (Gunn & Peterson 1965, Yoshii 
& Peterson 1994, Haiman 2002, Santos 2004) and the contribution to the CIB will be cut-off at 
A < 912(1 -|-zj) A. At longer wavelengths the flux, including emission from the nebulae around 
each star, will propagate without signifleant attenuation. Population III stars emit as black bodies 
at logT ~ 4.8 — 5 (Schaerer 2002) and the Lyman limit at z ~ 20 is shifted to ~ 2//m at z = 0. The 
2MASS J band filter contains emission out to 1.4 /xm and represents the shortest wavelength where 
excess in the CIB over that from ordinary galaxies has been measured (Kashlinsky Sz Odenwald 
2000, Cambresy et al 2001, Kashlinsky et al 2002, Odenwald et al 2003). If the measured excesses 
in CIB flTictuations and isotropic component at J band are indeed attributable to Population III, 
their era must have lasted until 14. 

Fig. 23 shows the resultant CIB fluctuations from Population III stars from 1.25 to 8 ^um 
. The fluctuations are normalized to reproduce the CIB excess at 2.2 /im of 10 nWm^^sr^^ . 
All Population III stars were assumed to start forming at z = 20, but difl^erently colored lines 
correspond to different values of Zf, the redshift of the end of Population III era. Because of the 
Lyman continuum absorption there would be no appreciable CIB fluctuations at J band if Zf ^14 
(no blue lines) and the contribution at J band will come from z ~ 13, rather than z = 20. Taken 
at face value, the presence of excess CIB fluctuations in J band indicates that the Population III 
stars were possibly the dominant sources of luminosity until Zf^lA. 



7.4. Can CIB anisotropies from Population III be measured? 



Both DIRBE and 2MASS data indicate CIB anisotropies at amplitudes larger than the con- 
tribution from ordinary galaxy populations (Kashlinsky & Odenwald 2000, Kashlinsky et al 2002, 
Odenwald et al 2003) and are consistent with significant contributions due to Population III. How- 
ever, because the measured signal contains contributions from remaining galaxies (all galaxies for 
DIRBE and K > 19 galaxies for 2MASS), it is difficult to isolate the contribution from the Popu- 
lation III stars. 




Fig. 23. CIB fluctuations, 'J cp-P^ jl-K , from Population III stars forming at 2* = 20 are shown with colored lines. Blue 

and green lines correspond to Population III era ending at = 15 and 10 respectively. Thick and thin colored lines correspond 
to the two extremes of Population III SED modeling after Santos et al (2002) and shown in Fig. 21. Solid and dashed colored 
lines correspond Population III forming in haloes with Tvir > 400K and Tvir > 2000K. Filled squares at 1.25 and 2.2 ^lra 
show residual atmospheric fluctuations from ground based 2MASS measurements after 1 hour of exposure from Odenwald et al 
(2003). Thick long dashes are adopted from Kashlinsky et al (2004) and denote the upper limit on zodiacal light fluctuations 
from (Abraham et al 1997) scaled to the corresponding band. Thick dashed-triple-dotted line denotes cirrus fluctuations: these 
arc upper limits at 2.2 micron and an estimate from (Kiss et al 2003) at 3.5 micron. Dot-dashed lines correspond to shot 
noise from galaxies fainter than AB magnitude = 22 (thickest), and rriAB = 20 (thin). Diamonds with error bars show the 
CIB fluctuations at ~ 0.7° from the COBE DIRBE fluctuations analysis (Kashlinsky & Odenwald 2000). Note that because 
of the large DIRBE beam, these results include contribution from all galaxies as well as other sources such as Population III. 
The K band CIB fluctuation from deep integration 2MASS data is shown at the largest scale accessible to those data with X ; 
the 2MASS data shown in the figure were taken for the the patches for which galaxies were removed brighter than ii" ~ 19 
(Kashlinsky, Odenwald et al 2002, Odenwald, Kashlinsky et al 2003). The cosmic variance 1-sigma uncertainty is shown with 
shaded regions for 1.25 /^m . The darkest shade corresponds to a total of 1,000 deg^ observed, and the lighter shade corresponds 
to a total of 100 dcg^ . 
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Population III stars, if massive, should have left a unique and measurable signature in the near- 

IR CIB anisotropies over angular scales from ~ 1 arcminute to several degrees as Fig. 23 shows. 
Detection of these fluctuations depends on the identification and removal of various foreground 
emission (and noise) contributions: atmosphere (for ground-based measurements), zodiacal light 
from the Solar system. Galactic cirrus emission, and instrument noise. 

Kashlinsky et al (2004) discuss observability of the Population III CIB fluctuations vs the 
confusion arising from the various foregrounds. The confusion may arise from 1) Galactic cirrus 
emission, 2) galaxies with ordinary stellar populations not removed from the data, 3) zodiacal 
light fluctuation, and 4) atmospheric fluctuations from ground measurements. They estimated the 
contributions from the various foregrounds to the fluctuations from 1 to 3.5 fim . The numbers 
from their discussion are plotted in the top panels of Fig. 23. 

Below is brief discussion of the various contaminants when searching for the CIB fluctuations 
from the Population III era: 

• Atmosphere: Filled squares in Fig. 23 show the residual atmospheric fluctuations at 1.25 and 
2.2 fim on sub-arcminute scales after one hour of integration from one of the deep 2MASS fields 
measurements (Odcnwald, Kashlinsky ct al 2003). Atmospheric gradients become important on 
larger angular scales (;^ 0.05 — 0.1°) and their efi^ects can be highly variable on a wide range of time 
scales (e.g. Adams Sz Skrutskie 1995, http://pegasus.phast.umass.edu/adams/airglowpage.html) 
making detection of arcminute and degree scale CIB fluctuations difficult in ground based obser- 
vations. Yet these problems can be completely avoided with space-based experiments. 

• Ordinary galaxies must be eliminated to sufficiently faint levels so that the remaining fluctu- 
ations in their cumulative emission are sufficiently small. The dashed-dotted lines in Fig. 23 show 
the shot-noise fiuctuations estimated from the observed deep counts from 1 to 8 fim . Thin lines 
correspond to galaxies brighter than ttiab = 20 identified and removed and thick lines when the 
same is possible out to uiab = 22. One can see that ideally one would have to identify galaxies 
out to mAB ^23 — 24 in order to measure the CIB fluctuations from Population III era across the 
entire range of scales. 

• Zodiacal light emission from interplanetary dust (IPD) is the brightest foreground at most 
IR wavelengths over most of the sky. There are some structures in this emission associated with 
particular asteroid families, comets, and an earth-resonant ring, but these structures tend to be 
confined to low ecliptic latitudes or otherwise localized. The main IPD cloud is generally modeled 
with a smooth density distribution. Observationally, intensity fluctuations of the main IPD cloud 
have been limited to < 0.2% at 25 //m (Abraham et al 1997). Extrapolating this limit to other 
DIRBE wavelengths using the observed mean high-latitude zodiacal light spectrum from COBE 
DIRBE (Kelsall et al 1998), Kashlinsky, Arendt et al (2004) arrive at the limits shown in Figure 
23. Because the Earth is moving with respect to (orbiting within) the IPD cloud, the zodiacal 
light varies over time. Likewise, any zodiacal light fluctuations will not remain fixed in celestial 
coordinates. Therefore repeated observations of a field on timescales of weeks to months should be 
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able to distinguish and reject any zodiacal light fluctuations from the invariant Galactic and CIB 

fluctuations. The upper limit on zodiacal light fluctuations can be evaluated also at the bottom 
panels of Fig. 23 using extrapolations of the DIRBE- measured zodiacal light spectrum from Fig. 
4. The zodiacal light is expected to be very smooth at all the near-IR bands. 

• Galactic cirrus. Intensity fluctuations of the Galactic foregrounds are perhaps the most 
difficult to distinguish from those of the CIB. Stellar emission may exhibit structure from binaries, 
clusters and associations, and from large scale tidal streams ripped from past and present dwarf 
galaxy satellites of the Milky Way. At long IR wavelengths, stellar emission is minimized by virtue 
of being far out on the Rayleigh- Jeans tail of the stellar spectrum (apart from certain rare classes 
of dusty stars). At near-IR wavelengths stellar emission is important, but with sufficient sensitivity 
and angular resolution most Galactic stellar emission, and related structure, can be resolved and 
removed. IR emission from the ISM (cirrus) is intrinsically diffuse and cannot be resolved. Cirrus 
emission is known to extend to wavelengths as short as 3 /xm (Arendt et al 1998). Statistically, 
the structure of the cirrus emission can be modeled with power-law distributions. Using the mean 
cirrus spectrum, measurements made for the cirrus fluctuations in the far-IR with ISO (Kiss et 
al 2003) were scaled all the way to 3.5 /im by Kashlinsky, Arendt et al (2004), providing the 
estimated fluctuation contribution from cirrus that is shown in Figure 23. The extrapolation to 
shorter wavelengths is highly uncertain, because cirrus (diffuse ISM) emission has not been detected 
at these wavelengths, and the effects of extinction may become more significant than those of 
emission, but cirrus contribution is generally expected to be several orders of magnitude lower than 
the CIB fluctuations (Arendt et al 1998). We do not show the expected cirrus fluctuations at the 
wavelengths corresponding to the bottom panels of Fig. 23 as the simple extrapolations of the cirrus 
SED from Fig. 4 may be inadequate due to the presence of PAH emission at these wavelengths. 

Figure 23 shows that CIB fluctuations from Population III would be the dominant source of 
diffuse light fluctuations on arcminute and degree scales even if Population III stars epoch was 
briefer, and their diffuse flux smaller, than the current CIB numbers suggest. Their angular power 
spectrum should be very different from other sources of diffuse emission and its measurement thus 
presents a way to actually discover Population III and measure the duration of their era and their 
spatial distribution. The latter would provide direct information on primordial power spectrum 
on scales and at epochs that are not easily attainable by conventional surveys. This measurement 
would be imperative to make and is feasible with the present day space technology. 

The previous discussion suggests that CIB fluctuations from the Population III era are ob- 
servable, but because of the large-scale atmospheric gradients, space observations are required for 
detection. The following is needed: 

1) In order to make certain that the signal is not contaminated by distant ordinary galaxies 
with normal stellar populations, one would need to conduct a deep enough survey in order to 
identify and eliminate ordinary galaxies from the field. In practice, as Fig. 9 indicates, going to 
ruAB — 24 would be sufficient. 
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Fig. 24. Black rectangular box shows the regions of the parameters required for a space mission to measure 

the CIB anisotropics from Population III. Cross-hatched region in the right panel displays the magnitude range where 
Population III stars are expected to dominate CIB anisotropics. Other colors correspond to the currently existing 
datasets or planned NASA missions which can be utilized for this purpose: red is the Hubble Deep Field from HST 
(http://www.stsci.edu/ftp/science/hdf/hdf.html), dark blue is for the deep 2MASS data (15; 20), green is the NASA's Spitzer 
GOODS project (http://www.stsci.edu/ftp/science/goods/), light blue corresponds to the planned NASA WISE MIDEX mis- 
sion (http://www.astro.ucla.edu/ wright/WISE), yellow corresponds to the field-of-view of the NASA James Webb Space 
Telescope planned for the end of the decade (http:/ /jwst. gsfc.nasa.gov) and orange color to the three surveys planned by SNAP 
(http: / /snap. lbl.gov). 



2) A direct signature of Population III signal is that there should be a Gunn-Peterson Lyman- 
break feature in the spectral energy distribution of the CIB fluctuations, i.e. at wavelengths < 
912(1 + Zf)~^ A. That spectral drop would provide an indication of the epoch corresponding to the 
end of the Population III era. At longer wavelengths the spectral energy distribution of the CIB 
fluctuations would probe the history of energy emission and its re-distribution by the concurrent 
IGM. At wavelengths 1 10;um zodiacal light fluctuations may become dominant. Hence one would 
want 0.9/im^A^5 — 10/xm. 

3) On angular scales from a few arcminutes to ~ 10° Population III would produce CIB 
anisotropics with a distinct and measurable angular spectrum, but its measurement will be limited 
by the cosmic or sampling variance. Calculations show that for reliable results on scales up to 9 
one would need to covers area a few times larger (Kashlinsky et al 2004). Shaded regions in the 
upper left panel of Fig. 23 show that in order to get reliable results on scales up to 0ii 1° one would 
need to observe areas of Z 10° across. 

Fig. 24 shows the parameters a space-based survey should have in order to probe Population 
III CIB anisotropics and compares it with the currently planned space missions. The left panel 
shows the projection on the (angular scale)-wavelength plane; the right panel shows the projection 
onto the (angular scale)-sensitivity plane. Clearly, such a Population III experiment needs to have 
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a wide field-of-view (FOV) and at the same time have good angular resolution (preferably up to 

a few arcseconds). The sensitivity requirements were taken to reach the ability to identify and 
eliminate from the data ordinary galaxies that at all the near-IR wavelengths contribute ^ 1 
nWm~^sr~^ ; Fig. 18 shows that this is reachable at mAB~24 which with modern technologies and 
a ~ 0.5 m mirror is achievable from space in ~ 1 hr integrations. None of the currently planned 
surveys cover perfectly the required range of parameters, although a combination of various missions 
instruments may do a part of the job. Cooray ct al (2004) propose a rocket experiment similar to 
an earlier experiment by Xu ct al (2002) with an FOV of a few degrees and ~ arcmin resolution, 
whereas Kashlinsky et al (2004) argue for a dedicated space mission based on the existing detector 
technologies. 

8. Snapshot of the future 

It is important to resolve as much of the CIB as possible into individual contributors and 
isolate the contributions to the total mean flux and anisotropies by the cosmic epoch. It is also 
important to isolate the epochs from which luminous sources contribute to the part of the CIB 
which remains unresolvable by contemporary instruments and techniques. Significant progress in 
these directions is expected to happen with the current (Spitzer) and planned space missions which 
should cover larger scales, have better angular resolution and go deeper across the entire range of 
the IR spectrum. In this section we give a brief overview of the IR missions and their potential for 
CIB studies. 

Spitzer 

The Spitzer telescope (initially known as SIRTF - Space InfraRed Space Facility) has been 
launched aboard a Boeing Delta II rocket in August 2003. The telescope carries three instruments: 
InfraRed Array Camera (IRAC), Multiband Imaging Photometer (MIPS) and InfraRed Spectro- 
graph (IRS). IRAC is a four-channel camera that provides simultaneous 5.12' x 5.12' images. The 
pixel scale is 1.2" in all IRAC bands. The Multiband Imaging Photometer for SIRTF (MIPS) is 
designed to provide very deep imaging and mapping at 24, 70, and 160 jum . In integrations of 2000 
seconds, it reaches 5— cr detection limits at these wavelengths of 0.2, 0.5, and 8 mJy, respectively. 
The pixel scale is 2.4" at 24 fim , 4.5" at 70 fim and 15" at 160 fim . 

On the observational side the following surveys have already been approved by the Spitzer 
Legacy ^ and Director's Discretionary ^ programs: 



'http://ssc.spitzer.caltech.edu/legacy/ 
' http : // ssc. spit zer . caltech.edu /fls/ 
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• First Look Survey: The extragalactic component of the First Look Survey (FLS) is intended to 

reach the 1-sigma sensitivities of ~ 10/iJ with IRAC and 1 to 27 mJy with MIPS over a 4 sq. deg. 
and a 1 sq. deg. field. A smaller 0.25 sq. deg. region will be covered with 10 times the nominal 
integration time. Although the survey is very shallow, the 160 /im data are already expected to 
be confusion limited even for this shallow survey. 

• SWIRE: The SWIRE project (Lonsdale et al 2004) will cover 7 fields between 5 and 15 sq. deg. 
in size with sensitivities that are several times better than the FLS. A stated goal of the SWIRE 
project is to investigate galaxy formation in the 0.5 < z < 2.5 range. 

• GOODS: The GOODS project will observe 300 square arcmin divided into two fields: HDF 
North and Chandra Deep Field South. Mean exposure time per position is 25 hours per band 
with IRAC. At 24 fim , MIPS will obtain 10 hour exposures, if on orbit experience shows that 
this will be a distinct advantage over planned Guaranteed Time Observations (GTO). MIPS GTO 
proposals for these regions at 70 and 160 are sufficiently deep, that the GOODS project will 
not duplicate these observations. If it will prove useful, the GOODS observations will also include 
an ultra-deep field in the HDF-N region with up to 100 hours of integration. The GOODS team 
intends to use the resolved galaxy counts, which should include L* galaxies to redshifts to z = 5, 
to establish the best lower limits on the 3.6 - 24 /xm mean GIB. 

On the data analysis side of the long term programs, the National Science Foundation has 
approved: 

• LIBRA, or Looking for Infrared Background Radiation Anisotropics, is an NSF-supported 5 
year project to measure GIB anisotropics from early epochs. The integrated light of all galaxies in 

the deepest counts to date fails to match the observed mean level of the GIB, indicating a significant 
high-redshift contribution to the GIB. The project will use the Spitzer and deep 2MASS data at 1 
- 160 fim to measure the spatial fiuctuations of this residual high-z portion of the GIB. Analysis 
of these spatial fluctuations will provide information on the luminosity of the early universe, and 
on the developing structure of galaxies and clusters of galaxies at early epochs. 



Astro-F 

Japan's Astro-F or the Infrared Imaging Surveyor (IRIS) is an infrared astronomy mission 
scheduled for launch in 2005 into a sun-synchronous polar orbit. IRIS employs a 70 cm telescope 
cooled to 6 K using liquid helium. ASTRO-F is designed for advanced surveys with two observation 
modes. The first is a general survey, one rotation in one orbital period, which is used for the all 
sky survey. The second is a pointing mode for imaging and spectroscopic observations of a limited 



http: / /www. stsci.edu /science / goods 

http: / /haiti. gsfc.nasa.gov/kashlinsky /LIBRA 
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region of the sky. Two focal-plane instruments are installed. One is Far-Infrared Surveyor (FIS) 
and the other is Infrared Camera (IRC). The FIS is a photometer optimized for all-sky survey 
with far-infrared arrays, and is expected to produce catalogs of infrared sources with much better 
sensitivity and higher angular resolution than the IRAS. The FIS can be operated as an imager or 
a Fourier-transform spectrometer in the pointing mode. The IRC is a three-channel camera that 
covers the wavelength bands from 2 fj,m. to 25 /xm , and has the capability to perform low-resolution 
spectroscopy with prisms/grisms on filter wheels. The field of view of the IRC is 10 arcmin and the 
spatial resolution is approximately 2 arcsec. Large format arrays are used to attain the deep survey 
with wide field and high angular resolution. IRC observations are carried out only in pointing 
mode. 

The detection limits will be 1 - 100 fiJy in the near-mid infrared and 10-100 mJy in the far 
infrared. It is planned to conduct a large area (~4 square degree) survey at the K and L bands for 
CIB studies. 

WISE 

WISE (Wide-field Infrared Survey Explorer) has been selected as NASA's next Medium- 
Class Explorer (MIDEX) mission with a tentative launch date in 2008. It will have four bands at 
3.5, 4.6, 12 and 23 /xm with c± 3" pixels. During its 6 month mission WISE satellite will provide 
an all-sky coverage about 1,000 more sensitive than IRAS. It should detect ULIRGs out to z ~ 3 
and bright IR galaxies out to z ^ 1. 

Herschel and Planck 

The European Space Agency's Herschel Space Observatory ^'^ (formerly called Far Infrared 
and Submillimetre Telescope, or FIRST), with an anticipated launch in 2007, will be the first 
space observatory covering the full far infrared and sub-millimeter waveband, and its passively 
cooled telescope will have a 3.5 meter diameter mirror. It will be located at L2 and will provide 
photometry and spectroscopy in the 57 to 670 ^m range (Pilbratt 2004). The Observatory will 
have 3 instruments: The Photodetector Array Camera and Spectrometer (PACS), the Spectral 
and Photometric Imaging REceiver (SPIRE) instrument, and the Heterodyne Instrument for the 
Far Infrared (HIFI) instrument. SPIRE is made of two sub-instruments: a three-band imaging 
photometer operating at 250, 360 and 520 /xm , and an imaging Fourier Transform Spectrometer 
(FTS) covering 200-670 /xm . The field of view of the photometer will be 4' x 8', observed simul- 
taneously in the three spectral bands. An internal beam steering mirror allows spatial modulation 



http://wise.ssl.berkeley.edu/ 
littp://astro.esa.int /herschel/ 
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of the telescope beam and will be used to jiggle the field of view in order to produce fully-sampled 

images. Observations can also be made by scanning the telescope without chopping. In addition to 
pointed imaging and spectroscopy of distant galaxies, it is planned that the telescope will conduct 
deep confusion limited surveys of > 100 deg2 of the sky. 

Planck satellite is designed mainly for CMB measurements and will be launched together 
with Herschel. Its High Frequency Instrument (HFI) has bands at 350, 550, 850 and 1380 /um with 
resolution of 10 to 5 arcminutes and will be useful for CIB studies. Its all sky survey will probably 
be too shallow to constrain galaxy evolution, but will give good measurements of the bright end 
of galaxy counts. With its scanning strategy Planck will survey some high Galactic latitude areas 
much deeper (a factor 3 in flux) and they could be useful for studies of the far-IR to sub-mm 
CIB fluctuations at arcminute to degree angular scales. 



JWST 

Successor to the Hubble Space Telescope, the James Webb Space Telescope is a large, 
infrared-optimized space telescope scheduled for launch in August, 2011 (Sabelhaus &: Decker 2004). 
JWST will have a large mirror, 6.5 meters in diameter. It will reside in an L2 Lissajous orbit. It 
will have two instruments of direct relevance to CIB science: the NIR camera (NIRCAM) and 
the MIR Instrument (MIRI). The NIRCAM will cover wavelengths from 0.6 to 5 /xm with 0.03 
arcsec/pixcl resolution from 0.6 to 2.3 /im and 0.064 arcsec/pixel at 2.3 to 5 fim . The MIRI will 
provide the JWST with imaging and spectroscopy at wavelengths from 5 through 27 microns with 
^ 0.2 arcsec/pixel. The FOV will be ~ 1.5' x 1.5'. The telescope will be able to detect Lyman 
break galaxies out to 2;^ 10 and supernovae out to 2; ^20 (Gardner et al 2004). Its science goals in 
cosmology will be to identify the end of dark ages and the assembly of galaxies at observer bands 
from visible to mid-IR. 



SNAP 

The Super Novae Acceleration Probe (SNAP), or the Joint Dark Energy Mission, is planned 
as a combined DOE and NASA space mission designed to probe the equation-of-statc of the Universe 
with high z supernovae out to z 2 (Aldcring et al 2004). It is scheduled tentatively for launch 
early in the next decade. SNAP will cover wavelengths from visible to 1.7 //m and will be useful 
for NIR CIB studies. Panoramic, wide and deep surveys are currently planned to cover up to 
~ 10^, 1000 and 15 deg^ and will go to AB magnitudes of ~ 26.5, 28 and 30 respectively. Fig. 24 



http://www.rssd.esa.int/Planck 

http://www.jwst.nasa.gov 

http://snap.lbl.gov 
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shows that these SNAP surveys at the longest SNAP wavebands should be very useful in uncovering 
the potential high- 2; contribution to the CIB fluctuations at <1.5 fim including from Population 
III. 



SCUBA-2, SOFIA, ALMA 

SCUBA-2 1^ will replace the current SCUBA on the James Clerk Maxwell Telescope in 2006. 
It will have a total of ~ 10, 000 bolometers and offer simultaneous imaging of a 50 arcmin^ FOV at 
450 and 850 //m mapping the sky up to 1,000 times faster than the current SCUBA array (Audley 
et al 2004). 

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is under development by 
NASA, and is expected to have its initial science operations in early 2007. This observatory, with 
its 2.5 m telescope, operates at 41,000 to 45,000 feet aboard a Boeing 747 aircraft, providing access 
to much of the far-IR spectrum. It will provide complementary capabilities to the Spitzcr Space 
Telescope, permitting higher spatial resolution, and hence lower confusion limits at far infrared 
wavelengths. The High Angular Resolution Widefield Camera (HAWC) provides imaging capability 
at 50, 90, 160, and 215 /xm . Surveys are planned to detect high redshift galaxies (D. A. Harper, 
private communication). The instrument will reach the confusion limit in about 40 hours, so such 
observations will be carried out over very limited areas. The instrument will be useful to follow up 
observations of high redshift sources detected by Spitzer or other surveys. 

The Atacama Large Millimeter Array (ALMA) in Chile is planned as a synthesis radio 
telescope built in international collaboration between North America and Eiirope which will operate 
at millimeter and sub-mm wavelengths. It will reach resolution of ^ 1" over fields of several square 
arcmin. 



SAFIR, SPIRIT, SPECS 

Planned for the end of the next decade these prospective NASA missions will achieve both 
high angular resolution (arcsecond to sub-arcsecond) and high sensitivity (~ AtJy) necessary to 
resolve and characterize most of the sources comprising the cosmic infrared background at far-IR 
and sub-mm wavelengths. 



http://www.jach.hawaii.edu/ J ACpublic/JCMT/scuba/scuba2/ 
^*http://sofia.arc.nasa.gov/Sofia/ 
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Recommended by the National Academy of Sciences Decadal Review SAFIR stands for a 
Single Aperture Far-Infrared observatory (Lester et al 2004). It is a large cryogenic space telescope 
scheduled for launch around 2015 to 2020. The SAFIR telescope will operate between 20 fim and 1 
mm and will be cooled to about 5 K. The combination of large mirror size and cold temperature will 
make SAFIR significantly more sensitive than the Spitzer and Herschel instruments with sensitivity 
limited only by the irreducible noise of photons in the astrophysical backgrounds. 

In the Community Plan for Far-IR/Submillimiter Space Astronomy it is envisaged that 
SAFIR will be followed by a kilometer baseline FIR interferometer designed to provide both wide 
field-of-view imaging and spectroscopy. This interferometer is widely known as the Sub-millimeter 
Probe of the Evolution of Cosmic Structure (SPECS). SAFIR and SPECS were selected for study 
under the ROSS/Vision Mission study program, and the Space Infrared Interferometer Telescope 
(SPIRIT), a science pathfinder for SPECS, was selected for study under the NASA ROSS Origins 
Science mission concept study program. These will be the first instruments that would achieve 
sub-arcsecond resolution at FIR/sub-mm wavelengths and provide the /xJy- level sensitivity over a 
field-of-view of 1 arcmin^ (Leisawitz 2004). All three observatories, in addition to high sensitivity 
and resolution, will have spectral resolution of the order of ~ 1000. 



9. Concluding remarks 

Cosmic Infrared Background presents a complementary, and sometimes the only, way to de- 
tect cumulative emission from galaxies at all cosmic times, including from objects inaccessible to 
current or future telescopic studies. Thanks to better detector sensitivities and new measurement 
techniques it is now possible to identify and remove foreground emissions at IR wavelengths and 
begin to measure the CIB. Because of the difficulties in identifying the mean level of the CIB 
in the presence of strong foreground emissions, it is often useful to attempt to measure the CIB 
fluctuations spectrum. The latter also allows to isolate the contributions to the CIB from progres- 
sively fainter galaxies and, on average, earlier cosmic epochs from surveys with ~ arcsec angular 
resolution. 

We have reviewed how mean levels of the CIB and, both the amplitude and the angular 
spectrum of, its anisotropics relate to the properties of the underlying galaxy populations. The 
various foreground contributions and the latest measurements of the mean CIB and its fluctuations 
were summarized from the near-IR to sub-mm. There are mutually consistent detections of the 



Astronomy and Astrophysics in the New MiUennium 2001, National Research Council 
http://www.nap.edu/books/0309070317/html/ 

■^^ http://safir.jpl.nasa.gov/ 

^^in Proceedings "New Concepts for Far-Infrared/Submillimiter Space Astronomy", eds. D.J. Benford and D. 
Leisawitz (Washington, DC: NASA), NASA CP-2003-212233, pp. XV-XXV (2003) 
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mean levels and fluctuations of the CIB in the near-IR, low upper limits in the mid-IR and detections 
of the mean CIB in the far-IR longward of 100 /xm . 

There are two broad classes of contributors to the CIB: from 'ordinary' (metal rich) stellar 
populations and from objects from the Population III era. The data on the present-day galaxy 
luminosity density present an important normalization point for interpretation of the CIB mea- 
surements. We showed that by constructing realistic SEDs for the stellar component one can get 
a fair agreements between the different luminosity density measurements at various wavelengths. 
This implies that it is unlikely that significant flux is unaccounted for in galaxy measurements at 
different bands and with different instruments and methods, and the results are consistent with the 
independently derived measurements of the present-day stellar density parameter, Q^,. The same 
holds true at longer wavelengths where emission from dust is the dominant contributor. 

Total flux from observed ordinary galaxy populations gives a lower limit on their contribu- 
tion to the CIB and can be estimated from a variety of galaxy counts data from the ncar-IR to 
sub-mm bands. We summarized such estimates based on various most recent measurements and 
conclude that 1) at far-IR wavelengths the observed background is most likely accounted for by 
the observed galaxy populations, which are probably located at z> 1 — 2; 2) at mid-IR wavelengths 
the total fluxes saturate at levels which are just below the best current upper limits on the CIB 
and probably account for at least a large fraction of the CIB; 3) at near-IR wavelengths the total 
contributions from the observed galaxy population saturate at levels significantly lower than the 
claimed detections of the CIB mean levels and anisotropics. The latter detections, if taken at 
face value, would indicate substantial contributions from much earlier epochs than probed by the 
observed faintest galaxy populations. 

A plausible candidate to explain the excess in the near-IR CIB would be emissions from 
Population III era objects. Their expected contributions should have left a unique and observable 
signature in the spectrum of the CIB anisotropics if Population III were massive stars as is currently 

expected. The contributions from the various foreground emissions are such that this signature can, 
and should be, observed with future space-based instruments and properly tuned observations. 

Finally, the current (Spitzer) and planned space missions and new instruments in the IR and 
sub-mm are expected to bring a wealth of high quality observational data. This data can and, no 
doubt, will be used for further progress in the CIB studies and in identifying the history of light 
emission from collapsed objects in the early Universe. 
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